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(57) Abstract: A method and system of physically solving the charge, mass, and current density functions of polyatomic molecules, 
polyatomic molecular ions, diatomic molecules, molecular radicals, molecular ions, or any portion of these species using Maxwell's 
equations and computing and rendering the physical nature of the chemical bond using the solutions. The results can be displayed 
on visual or graphical media. The display can be static or dynamic such that electron motion and specie's vibrational, rotational, and 
translational motion can be displayed in an embodiment. The displayed information is useful to anticipate reactivity and physical 
properties. The insight into the nature of the chemical bond of at least one specie can permit the solution and display of those of 
other species to provide utility to anticipate their reactivity and physical properties. 
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SYSTEM AND METHOD OF COMPUTING AND RENDERING THE NATURE OF 
POLYATOMIC MOLECULES AND POLYATOMIC MOLECULAR IONS 

This application claims priority to U.S. Application Nos.; 60/73 0,882, filed October 28, 
2005; 60/732,154, filed November 2, 2005; 60/737,744, filed November 18, 2005; 
60/758,528, filed January 13, 2006; 60/780,518, filed March 9, 2006; 60/788,694, filed April 
4, 2006; 60/812,590, filed June 12, 2006; and 60/815,253, June 21, 2006, the complete 
disclosures of which are incorporated herein by reference. 

Field of the Invention: 

This invention relates to a system and method of physically solving the charge, mass, 
and current density functions of polyatomic molecules, polyatomic molecular ions^ diatomic 
molecules, molecular radicals, molecular ions, or any portion of these species, and computing 
and rendering the nature of these species using the solutions. The results can be displayed on 
visual or graphical media. The displayed information provides insight into the nature of these 
species and is useful to anticipate their reactivity, physical properties, and spectral absorption 
and emission, and permits the solution and display of other species. 

Rather than using postulated unverifiable theories that treat atomic particles as if they 
were not real, physical laws are now applied to atoms and ions. In an attempt to provide 
some physical insight into atomic problems and starting with the same essential physics as 
Bohr of the e~ moving in the Coulombic field of the proton with a true wave equation, as 
opposed to the diffusion equation of Schrodinger, a classical approach is explored which 
yields a model that is remarkably accurate and provides insight into physics on the atomic 
level. The proverbial view deeply seated in the wave-particle duality notion that there is no 
large-scale physical counterpart to the nature of the electron is shown not to be correct. 
Physical laws and intuition may be restored when dealing with the wave equation and 
quantum atomic problems. 

Specifically, a theory of classical quantum mechanics (CQM) was derived from first 
principles as reported previously [reference Nos. 1-8] that successfully applies physical laws 
to the solution of atomic problems that has its basis in a breakthrough in the understanding of 
the stability of the bound electron to radiation. Rather than using the postulated Schrodinger 
boundary condition: Ml F -> 0 as r -» oo which leads to a purely mathematical model of the 
electron, the constraint is based on experimental observation. Using Maxwell's equations, the 
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classical wave equation is solved with the constraint that the bound n = l -state electron 
cannot radiate energy. Although it is well known that an accelerated point particle radiates, 
an extended distribution modeled as a superposition of accelerating charges does not have to 
radiate. A simple invariant physical model arises naturally wherein the predicted results are 
extremely straightforward and internally consistent requiring minimal math, as in the case of 
the most famous equations of Newton, Maxwell, Einstein, de Broglie, and Planck on which 
the model is based. No new physics is needed; only the known physical laws based on direct 
observation are used. 

Applicant's previously filed WO2005/067678 discloses a method and system of 
physically solving the charge, mass, and current density functions of atoms and atomic ions 
and computing and rendering the nature of these species using the solutions. The complete 
disclosure of this published PCT application is incorporated herein by reference. 

Applicant's previously filed WO2005/1 16630 discloses a method and system of 
physically solving the charge, mass, and current density functions of excited states of atoms 
and atomic ions and computing and rendering the nature of these species using the solutions. 
The complete disclosure of this published PCT application is incorporated herein by 
reference. 

Applicant's previously filed U.S. Published Patent Application No. 20050209788 Al, 
relates to a method and system of physically solving the charge, mass, and current density 
functions of hydrogen-type molecules and molecular ions and computing and rendering the 
nature of the chemical bond using the solutions. The complete disclosure of this published 
application is incorporated herein by reference. 

Background of the Invention 

The old view that the electron is a zero or one-dimensional point in an all-space 
probability wave function ^(x) is not taken for granted. The theory of classical quantum 
mechanics (CQM), derived from first principles, must successfully and consistently apply 
physical laws on all scales [1-8]. Stability to radiation was ignored by all past atomic models. 
Historically, the point at which QM broke with classical laws can be traced to the issue of 
nonradiation of the one electron atom. Bohr just postulated orbits stable to radiation with the 
further postulate that the bound electron of the hydrogen atom does not obey Maxwell's 
equations — rather it obeys different physics [1-12]. Later physics was replaced by "pure 
mathematics" based on the notion of the inexplicable wave-particle duality nature of electrons 
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which lead to the Schrcidinger equation wherein the consequences of radiation predicted by 
Maxwell's equations were ignored. Ironically, Bohr, Schrodinger, and Dirac used the 
Coulomb potential, and Dirac used the vector potential of Maxwell's equations. But, all 
ignored electrodynamics and the corresponding radiative consequences. Dirac originally 
attempted to solve the bound electron physically with stability with respect to radiation 
according to Maxwell's equations with the further constraints that it was relativistically 
invariant and gave rise to electron spin [13]. He and many founders of QM such as 
Sommerfeld, Bohm, and Weinstein wrongly pursued a planetary model, were unsuccessful, 
and resorted to the current mathematical-probability-wave model that has many problems [9- 
1 6], Consequently, Feynman for example, attempted to use first principles including 
Maxwell's equations to discover new physics to replace quantum mechanics [17]. 

Physical laws may indeed be the root of the observations thought to be "purely 
quantum mechanical 1 ' 5 and it was a mistake to make the assumption that Maxwell's 
electrodynamic equations must be rejected at the atomic level. Thus, in the present approach, 
the classical wave equation is solved with the constraint that a bound n = 1 -state electron 
cannot radiate energy. 

Herein, derivations consider the electrodynamic effects of moving charges as well as 
the Coulomb potential, and the search is for a solution representative of the electron wherein 
there is acceleration of charge motion without radiation. The mathematical formulation for 
zero radiation based on Maxwell's equations follows from a derivation by Haus [18]. The 
function that describes the motion of the electron must not possess spacetime Fourier 
components that are synchronous with waves traveling at the speed of light. Similarly, 
nonradiation is demonstrated based on the electron's electromagnetic fields and the Poynting 
power vector. 

It was shown previously [1-8] that CQM gives closed form solutions for the atom 
including the stability of the n = 1 state and the instability of the excited states, the equation 
of the photon and electron in excited states, and the equation of the free electron and photon, 
which predict the wave particle duality behavior of particles and light. The current and 
charge density functions of the electron may be directly physically interpreted. For example, 
spin angular momentum results from the motion of negatively charged mass moving 
systematically, and the equation for angular momentum, r x p, can be applied directly to the 
wave function (a current density function) that describes the electron. The magnetic moment 
of a Bohr magneton, Stern Gerlach experiment, g factor, Lamb shift, resonant line width and 
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shape, selection rules, correspondence principle, wave particle duality, excited states, reduced 
mass, rotational energies, and momenta, orbital and spin splitting, spin-orbital coupling, 
Knight shift, and spin-nuclear coupling, and elastic electron scattering from helium atoms, are 
derived in closed-form equations based on Maxwell's equations. The calculations agree with 
experimental observations. 

The Schrodinger equation gives a vague and fluid model of the electron. Schrodinger 
interpreted * (x) x P(x) as the charge-density or the amount of charge between x and 
x + dx Q¥ * is the complex conjugate of ). Presumably, then, he pictured the electron to 
be spread over large regions of space. After Schrodinger' s interpretation, Max Bom, who 
was working with scattering theory, found that this interpretation led to inconsistencies, and 
he replaced the Schrodinger interpretation with the probability of finding the electron 
between x and x + dx as 



Born's interpretation is generally accepted. Nonetheless, interpretation of the wave function 
is a never-ending source of confusion and conflict. Many scientists have solved this problem 
by conveniently adopting the Schrodinger interpretation for some problems and the Born 
interpretation for others. This duality allows the electron to be everywhere at one time: — yet 
have no volume. Alternatively, the electron can be viewed as a discrete particle that moves 
here and there (from r = 0 to r = oo ) 5 and X P K P * gives the time average of this motion. 

In contrast to the failure of the Bohr theory and the nonphysical, adjustable-parameter 
approach of quantum mechanics, multielectron atoms [1,4] and the nature of the chemical 
bond [1, 5] are given by exact closed-form solutions containing fundamental constants only. 
Using the nonradiative wave equation solutions that describe the bound electron having 
conserved momentum and energy, the radii are determined from the force balance of the 
electric, magnetic, and centrifugal forces that corresponds to the minimum of energy of the 
system. The ionization energies are then given by the electric and magnetic energies at these 
radii. The spreadsheets to calculate the energies from exact solutions of one through twenty- 
electron atoms are given in f 06 Mills GUT [1] and are available from the internet [19]. For 
400 atoms and ions, as well as hundreds of molecules, the agreement between the predicted 
and experimental results is remarkable. 

The background theory of classical quantum mechanics (CQM) for the physical 
solutions of atoms and atomic ions is disclosed in R. Mills, The Grand Unified Theory of 
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Classical Quantum Mechanics, January 2000 Edition, BlackLight Power, Inc., Cranbury, 
New Jersey, (" '00 Mills GUT"), provided by BlackLight Power, Inc., 493 Old Trenton Road, 
Cranbury, NJ, 08512; R. Mills, The Grand Unified Theory of Classical Quantum Mechanics, 
September 2001 Edition, BlackLight Power, Inc., Cranbury, New Jersey, Distributed by 
Axnazon.com (" '01 Mills GUT"), provided by BlackLight Power, Inc., 493 Old Trenton 
Road, Cranbury, NJ, 08512; R. Mills, The Grand Unified Theory of Classical Quantum 
Mechanics, July 2004 Edition, BlackLight Power, Inc., Cranbury, New Jersey, (" '04 Mills 
GUT"), provided by BlackLight Power, Inc., 493 Old Trenton Road, Cranbury, NJ, 08512; R. 
Mills, The Grand Unified Theory of Classical Quantum Mechanics, January 2005 Edition, 
BlackLight Power, Inc., Cranbury, New Jersey, (" '05 Mills GUT"), provided by BlackLight 
Power, Inc., 493 Old Trenton Road, Cranbury, NJ, 08512 (posted at 

www.blacklightT3ower.com) : R. L. Mills, "The Grand Unified Theory of Classical Quantum 
Mechanics", June 2006 Edition, Cadmus Professional Communications-Science Press 
Division, Ephrata, PA, ISBN 0963517171, Library of Congress Control Number 
2005936834, (" '06 Mills GUT"), provided by BlackLight Power, Inc., 493 Old Trenton 
Road, Cranbury, NJ, 08512 (posted at 

http.7/www.blacklightpower.com/bookdownload.shtml) : in prior PCT applications 
PCT/US02/35872; PCT/US02/06945; PCT/US02/06955; PCT/US01/09055; PCT/US01/ 
25954; PCT/US00/20820; PCT/US00/20819; PCT/US00/09055; PCT/US99/17171; 
PCT/US99/17129; PCT/US 98/22822; PCT/US98/14029; PCT/US96/07949; 
PCT/US94/02219; PCT/US9 1/08496; PCT/US90/01998; andPCT/US89/05037 and U.S. 
Patent No. 6,024,935; the entire disclosures of which are all incorporated herein by reference 
(hereinafter "Mills Prior Publications"). 

The following list of references, which are also incorporated herein by reference in 
their entirety, are referred to in the above sections using [brackets]: 

1. R. L. Mills, "The Grand Unified Theory of Classical Quantum Mechanics", June 2006 
Edition, Cadmus Professional Communications-Science Press Division, Ephrata, PA, 
ISBN 0963517171, Library of Congress Control Number 2005936834; posted at 
http://www.blacklightDQwer.coin/bookdownload.shtml . 

2. R. L. Mills, "Classical Quantum Mechanics Physics Essays, Vol. 16, No. 4, December, 
(2003), pp. 433-498; posted with spreadsheets at 
www.blacklightpower.com/techpapers.shtml. 
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3. R. Mills, "Physical Solutions of the Nature of the Atom, Photon, and Their Interactions to 
Form Excited and Predicted Hydrino States" , submitted. 

4. R. L. Mills, "Exact Classical Quantum Mechanical Solutions for One- Through Twenty- 
Electron Atoms", in press, posted with spreadsheets at 
http://www.blacklightpower.com/techpapers,shtml. 

5. R. L. Mills, "The Nature of the Chemical Bond Revisited and an Alternative Maxwellian 
Approach", Physics Essays, Vol. 17, (2004), pp. 342-389, posted with spreadsheets at 
http:y/www.blacklightpower.com/techpapers.shtmL 

6. R, L. Mills, "Maxwell's Equations and QED: Which is Fact and Which is Fiction", in 
press, posted with spreadsheets at http:/ywww.blacklightpower.com/techpapers.shtml. 

7. R. L. Mills, "Exact Classical Quantum Mechanical Solution for Atomic Helium Which 
Predicts Conjugate Parameters from a Unique Solution for the First Time", submitted, 
posted with spreadsheets at http://www.blacldightpowerxom/theory/theory.shtml.. 

8. R. Mills, "The Grand Unified Theory of Classical Quantum Mechanics", Int. J. Hydrogen 
Energy, Vol. 27, No. 5, (2002), pp. 565-590. 

9. R. L. Mills, "The Fallacy of Feynman's Argument on the Stability of the Hydrogen Atom 
According to Quantum Mechanics", Annales de la Fondation Louis de Broglie, Vol. 30, 
No. 2, (2005), pp. 129-151, posted at http://www.blacklightpower.com/techpapers.shtml. 

10. R. Mills, The Nature of Free Electrons in Superfluid Helium — a Test of Quantum 
Mechanics and a Basis to Review its Foundations and Make a Comparison to Classical 
Theory, Int. J. Hydrogen Energy, Vol. 26, No. 10, (2001), pp. 1059-1096. 

11. R. Mills, "The Hydrogen Atom Revisited", Int. J. of Hydrogen Energy, Vol. 25, Issue 12, 
December, (2000), pp. 1171-1183. 

12. F. Laloe, Do we really understand quantum mechanics? Strange correlations, paradoxes, 
and theorems, Am. J. Phys. 69 (6), June 2001, 655-701. 

13. P. Pearle, Foundations of Physics, "Absence of radiationless motions of relativistically 
rigid classical electron", Vol. 7, Nos. 11/12, (1977), pp. 931-945. 

14. V. F. Weisskopf, Reviews of Modern Physics, Vol. 21, No. 2, (1949), pp. 305-315. 

15. H. Wergeland, "The Klein Paradox Revisited", Old and New Questions in Physics, 
Cosmology, Philosophy, and Theoretical Biology, A. van der Merwe, Editor, Plenum 
Press, New York, (1983), pp. 503-515. 

16. A. Einstein, B. Podolsky, N. Rosen, Phys. Rev., Vol. 47, (1935), p. 777. 

17. F. Dyson, "Feynman's proof of Maxwell equations", Am. J. Phys., Vol. 58, (1990), pp. 



6 



WO 2007/053486 



PCT/US2006/042140 



209-211. 

18. Hans, H. A., "On the radiation from point charges", American Journal of Physics, 54, 
(1986), pp. 1126-1129. 

19. http://www.blackligh^)ower.com/new.shtml. 

SUMMARY OF THE INVENTION 

The present invention, an exemplary embodiment of which is also referred to as 
Millsian software, stems from a new fundamental insight into the nature of the atom. 
Applicant's new theory of Classical Quantum Mechanics (CQM) reveals the nature of atoms 
and molecules using classical physical laws for the first time. As discussed above, traditional 
quantum mechanics can solve neither multi-electron atoms nor molecules exactly. By 
contrast, CQM produces exact, closed-form solutions containing physical constants only for 
even the most complex atoms and molecules. 

The present invention is the first and only molecular modeling program ever built on 
the CQM framework. All the major functional groups that make up most organic molecules 
have been solved exactly in closed-form solutions with CQM. By using these functional 
groups as building blocks, or independent units, a potentially infinite number of organic 
molecules can be solved. As a result, the present invention can be used to visualize the exact 
3D structure and calculate the heat of formation of almost any organic molecule. 

For the first time, the significant building-block molecules of chemistry have been 
successfully solved using classical physical laws in exact closed-form equations having 
fundamental constants only. The major functional groups have been solved from which 
molecules of infinite length can be solved almost instantly with a computer program. The 
predictions are accurate within experimental error for over 375 exemplary molecules. 

Applicant's CQM is the theory that physical laws (Maxwell's Equations, Newton's 
Laws, Special and General Relativity) must hold on all scales. The theory is based on an often 
overlooked result of Maxwell's Equations, that an extended distribution of charge may, under 
certain conditions, accelerate without radiating. This "condition of no radiation" is invoked to 
solve the physical structure of subatomic particles, atoms, and molecules. 

In exact closed-form equations with physical constants only, solutions to thousands of 
known experimental values arise that were beyond the reach of previous outdated theories. 
These include the electron spin, g-factor, multi-electron atoms, excited states, polyatomic 
molecules, wave-particle duality and the nature of the photon, the masses and families of 
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fundamental particles, and the relationships between fundamental laws of the universe that 
reveal why the universe is accelerating as it expands. CQM is successful to over 85 orders of 
magnitude, from the level of quarks to the cosmos. Applicant now has over 65 peer-reviewed 
journal articles and also books discussing the CQM and supporting experimental evidence. 

The molecular modeling market was estimated to be a two-billion-dollar per year 
industry in 2002, with hundreds of millions of government and industry dollars invested in 
computer algorithms and supercomputer centers. This makes it the largest effort of 
computational chemistry and physics. 

The present invention's advantages over other models includes: Rendering true 
molecular structures; Providing precisely all characteristics, spatial and temporal charge 
distributions and energies of every electron in every bond, and of every bonding atom; 
Facilitating the identification of biologically active sites in drugs; and Facilitating drug 
design. 

An objective of the present invention is to solve the charge (mass) and current-density 
functions of polyatomic molecules, polyatomic molecular ions, diatomic molecules, 
molecular radicals, molecular ions, or any portion of these species from first principles. In an 
embodiment, the solution for the polyatomic molecules, polyatomic molecular ions s diatomic 
molecules, molecular radicals, molecular ions, or any portion of these species is derived from 
Maxwell's equations invoking the constraint that the bound electron before excitation does 
not radiate even though it undergoes acceleration. 

Another objective of the present invention is to generate a readout, display, or image 
of the solutions so that the nature of polyatomic molecules, polyatomic molecular ions, 
diatomic molecules, molecular radicals, molecular ions, or any portion of these species be 
better understood and potentially applied to predict reactivity and physical and optical 
properties. 

Another objective of the present invention is to apply the methods and systems of 
solving the nature of polyatomic molecules, polyatomic molecular ions, diatomic molecules, 
molecular radicals, molecular ions, or any portion of these species and their rendering to 
numerical or graphical form to all atoms and atomic ions. 

These objectives and other objectives are obtained by a system of computing and 
rendering the nature of at least one specie selected from a group of diatomic molecules 
having at least one atom is other than hydrogen, polyatomic molecules, molecular ions, 
polyatomic molecular ions, or molecular radicals, or any functional group therein, comprising 
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physical, Maxwellian solutions of charge, mass, and current density functions of said specie, 
said system comprising processing means for processing physical, Maxwellian equations 
representing charge, mass, and current density functions of said specie; and an output device 
in communication with the processing means for displaying said physical, Maxwellian 
solutions of charge, mass, and current density functions of said specie. 

Also provided is a composition of matter comprising a plurality of atoms, the 
improvement comprising a novel property or use discovered by calculation of at least one of 
a bond distance between two of the atoms, a bond angle between three of the atoms, and 
a bond energy between two of the atoms, orbital intercept distances and angles,charge-density 
functions of atomic, hybridized, and molecular orbitals, the bond distance, bond angle, and 
bond energy being calculated from physical solutions of the charge, mass, and current density 
functions of atoms and atomic ions, which solutions are derived from Maxwell's equations 
using a constraint that a bound electron(s) does not radiate under acceleration. 

The presented exact physical solutions for known species of the group of polyatomic 
molecules, polyatomic molecular ions, diatomic molecules, molecular radicals, molecular 
ions, or any functional group therein, can be applied to other species. These solutions can be 
used to predict the properties of other species and engineer compositions of matter in a 
manner which is not possible using past quantum mechanical techniques. The molecular 
solutions can be used to design synthetic pathways and predict product yields based on 
equilibrium constants calculated from the heats of formation. Not only can new stable 
compositions of matter be predicted, but now the structures of combinatorial chemistry 
reactions can be predicted. 

Pharmaceutical applications include the ability to graphically or computationally 
render the structures of drugs that permit the identification of the biologically active parts of 
the specie to be identified from the common spatial charge-density functions of a series of 
active species. Novel drugs can now be designed according to geometrical parameters and 
bonding interactions with the data of the structure of the active site of the drug. 

The system can be used to calculate conformations, folding, and physical properties, 
and the exact solutions of the charge distributions in any given specie are used to calculate the 
fields. From the fields, the interactions between groups of the same specie or between groups 
on different species are calculated wherein the interactions are distance and relative 
orientation dependent. The fields and interactions can be determined using a finite-element- 
analysis approach of Maxwell's equations. 
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Embodiments of the system for performing computing and rendering of the nature of 
the polyatomic molecules, polyatomic molecular ions, diatomic molecules, molecular 
radicals, molecular ions, or any portion of these species using the physical solutions may 
comprise a general purpose computer. Such a general purpose computer may have any 
number of basic configurations. For example, such a general purpose computer may 
comprise a central processing unit (CPU), one or more specialized processors, system 
memory, a mass storage device such as a magnetic disk, an optical disk, or other storage 
device, an input means, such as a keyboard or mouse, a display device, and a printer or other 
output device. A system implementing the present invention can also comprise a special 
purpose computer or other hardware system and all should be included within its scope. 

Although not preferred, any of the calculated and measured values and constants 
recited in the equations herein can be adjusted, for example, up to + 10%, if desired. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 illustrates a block diagram of an exemplary software program; and 
Figs. 2 and 3 illustrate pictures of an exemplary software program. 
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DETAILED DESCRIPTION OF INVETTOlsT 

Organic molecules comprising an arbitrary number of atoms can be solved using the 
similar principles and procedures as those used to solve alkanes of arbitrary length. Alkanes 
can be considered to be comprised of the functional groups of CH 3 , CH 2 , and C - C . These 
groups with the corresponding geometrical parameters and energies can be added as a linear 
sum to give the solution of any straight chain alkane as shown in the Continuous-Chain 
Alkanes section. Similarly, the geometrical parameters and energies of all functional groups 
such as alkanes, branched alkanes, alkenes, branched alkenes, alkynes, alkyl fluorides, alkyl 
chlorides, alkyl bromides, alkyl iodides, alkene halides, primary alcohols, secondary alcohols, 
tertiary alcohols, ethers ? primary amines, secondary amines, tertiary amines, aldehydes, 
ketones, carboxylic acids, carboxylic esters, amides, N-alkyl amides, N,N-dialkyl amides, 
urea, acid halides, acid anhydrides, nitriles, thiols, sulfides, disulfides, sulfoxides, sulfones, 
sulfites, sulfates, nitro alkanes, nitrites, nitrates, conjugated polyenes, aromatics, heterocyclic 
aromatics, substituted aromatics, and others can be solved. The functional-group solutions 
can be made into a linear superposition and sum, respectively, to give the solution of any 
organic molecule. The solutions of the functional groups can be conveniently obtained by 
using generalized forms of the geometrical and energy equations. The equations and sections 
that are referenced by not contained in this text refer to those sections and equations of the 
book by R. L. Mills entitled, "The Grand Unified Theory of Classical Quantum Mechanics", 
June 2006 Edition, Cadmus Professional Communications-Science Press Division, Ephrata, 
PA, ISBN 0963517171, Library of Congress Control Number 2005936834; posted at 
http://ww w.blacklightpower.com/bookdownload.shtml which is incorporated in its 
entirety by reference. ) 

Consider the case wherein at least two atomic orbital hybridize as a linear 
combination of electrons at the same energy in order to achieve a bond at an energy 
minimum, and the sharing of electrons between two or more such orbitals to form a MO 
permits the participating hybridized orbitals to decrease in energy through a decrease in the 
radius of one or more of the participating orbitals. The force generalized constant k* of a 
H 2 -type ellipsoidal MO due to the equivalent of two point charges of at the foci is given by: 
C x C 2 2e 2 

* • r«.~ <•*•» 
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where C x is the fraction of the H 2 -type ellipsoidal MO basis function of a chemical bond of 
the molecule or molecular ion which is 0.75 (Eq. (13.59)) in the case of H bonding to a 
central atom and 0.5 (Eq. (14.152)) otherwise, and C 2 is the factor that results in an 
equipotential energy match of the participating at least two molecular or atomic orbitals of the 
chemical bond. From Eqs. (13.58-13.63), the distance from the origin of the MO to each 
focus c' is given by: 



The internuclear distance is 



(15.2) 



2c' = 2. h^r ( 15 - 3 > 
\2Cf 2 

The length of the semiminor axis of the prolate spheroidal MO b = c is given by 



And, the eccentricity, e , is 



a 



(15.5) 



From Eqs. (1 1.207-1 1.212), the potential energy of the two electrons in the central field of the 
nuclei at the foci is 



-2e 2 a + 4a 2 -b 2 n » 

V e = n.c,c 2 , In , (15.6) 

' 1 l2 %7teja 2 -b 2 a-4a 2 -b 1 
The potential energy of the two nuclei is 

v fl , (15.7) 

" • %7te 4 'a 2 -b 2 

o 

The kinetic energy of the electrons is 

h 2 , a + \la 2 -b 2 , K(A 

T = n ^ ; — rr^ 111 — rr^ (15 * 8) 

2m a 2 -b 2 a-4<r -b 2 
And, the energy, V m , of the magnetic force between the electrons is 



F « = w ^ — rr^ ln — tt=^ (15,9) 

Ama^a 2 -b 2 a-y]a ~b 2 
The total energy of the if-type prolate spheroidal MO, E T (h^o), is given by the sum of the 
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energy terms: 

E T ( ff2 MoyV+T + V m + V p (15.10) 



JT { \ ^ 

L {h 2 mo )— — - 



C l C 2 



^ a) a-^Ja 2 -b 2 



S7te Q c* 



^12— *|ln ;-l 

\ a J a- 



(15.11) 



c 

where n x is the number of equivalent bonds of the MO and applies in the case of functional 
groups. In the case of independent MOs not in contact with the bonding atoms, the terms 
based on charge are multiplied by c BO , the bond-order factor. It is 1 for a single bond, 4 for 

an independent double bond as in the case of the C0 2 and N0 2 molecules, and 9 for an 
independent triplet bond. Then, the kinetic energy term is multiplied by c % BO which is 1 for a 
single bond, 2 for a double bond, and 9/2 for a triple bond. c x is the fraction of the H 2 -type 
ellipsoidal MO basis function of an MO which is 0.75 (Eqs. (13.67-13.73)) in the case of H 
bonding to an unhybridized central atom and 1 otherwise, and c 2 is the factor that results in 
an equipotential energy match of the participating the MO and the at least two atomic orbitals 
of the chemical bond. Specifically, to meet the equipotential condition and energy matching 
conditions for the union of the H 2 -type-ellipsoidal-MO and the HOs or AOs of the bonding 

atoms, the factor c 2 of a i7 2 4ype ellipsoidal MO may given by (i) one, (ii) the ratio of the 

Coulombic or valence energy of the AO or HO of at least one atom of the bond and 
13.605804 eV , the Coulombic energy between the electron and proton of H , (iii) the ratio of 
the valence energy of the AO or HO of one atom and the Coulombic energy of another, (iv) 
the ratio of the valence energies of the AOs or HOs of two atoms, (v) the ratio of two c 2 

factors corresponding to any of cases (ii)-(iv), and (vi) the product of two different c 2 factors 

corresponding to any of the cases (i)-(v). Specific examples of the factor c 2 of a H 2 -type 

ellipsoidal MO given in previous sections are 

0.936127 , the ratio of the ionization energy of N 14.53414 eV and 13.605804 eV , the 
Coulombic energy between the electron and proton of H ; 

0.91771 , the ratio of 14.82575 eV , -E Couhmb (c,2,sp 3 ) 5 and 13.605804 eV ; 
0.87495 , the ratio of 15.55033 eV , ~E Coulomb (c /W ,2^ 3 ), and 13.605804 eV ; 
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0.85252 , the ratio of 15.95955 eV , -E Coyhmb (c ,^ we ,2^ 3 ), and 13.605804 eV ; 
0.85252 , the ratio of 15.95955 eV , -E Couhmb (C bemme ,2$p 3 y and 13.605804 eF, and 

0.86359 , the ratio of 15.55033 eV , -E Couhmb (C a/ ^ e ,2^? 3 ), and 13.605804 eV . 

In the generalization of the hybridization of at least two atomic-orbital shells to form a 
shell of hybrid orbitals, the hybridized shell comprises a linear combination of the electrons 
of the atomic-orbital shells. The radius of the hybridized shell is calculated from the total 
Coulombic energy equation by considering that the central field decreases by an integer for 
each successive electron of the shell and that the total energy of the shell is equal to the total 

Coulombic energy of the initial AO electrons. The total energy E T (atom, msp 3 ^ (m is the 

integer of the valence shell) of the AO electrons and the hybridized shell is given by the sum 
of energies of successive ions of the atom over the n electrons comprising total electrons of 
the at least one AO shell. 

E T (atom,msp 3 }= ~Y, Ip m (15.12) 



where IP m is the mth ionization energy (positive) of the atom. The radius r^ 3 of the 

<Z-q)e 



msp- 

hybridized shell is given by: 

Z-l fry N 2 



* 3=Z ~ / ' 3 x (15.13) 

msp q =z-n %tus q E t [atom, msp J 

Then, the Coulombic energy E CoulQmb ^tom^msp^ of the outer electron of the atom msp" 
shell is given by 

/ \ -e 2 

E coui omb {atom,msp 3 y g ^ (15.14) 

o mS p^> 

In the case that during hybridization at least one of the spin-paired AO electrons is unpaired 
in the hybridized orbital (HO) 3 the energy change for the promotion to the unpaired state is 
the magnetic energy E(magnetic) at the initial radius r of the AO electron: 

Eirnagnetic) = = ^ (15.15) 

m r r 

e 

Then, the energy E (atom, msp 3 ^ of the outer electron of the atom msp 3 shell is given by the 
sum of E Gotllomb (atom,msp 3 ^ and E(magnetic): 
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E(atom,msp 3 )= "** + lKf2 ff (15.16) 
v 7 S7te n r 3 mrr 5 

« msp* e 

Consider next that the at least two atomic orbitals hybridize as a linear combination of 
electrons at the same energy in order to achieve a bond at an energy minimum with another 
atomic orbital or hybridized orbital. As a further generalization of the basis of the stability of 
the MO. the sharing of electrons between two or more such hybridized orbitals to form a MO 
permits the participating hybridized orbitals to decrease in energy through a decrease in the 
radius of one or more of the participating orbitals. In this case, the total energy of the 

hybridized orbitals is given by the sum of E (atom, msp 3 ) and the next energies of successive 

ions of the atom over the n electrons comprising the total electrons of the at least two initial 

AO shells. Here, E(atom,msp 3 ^ is the sum of the first ionization energy of the atom and the 

hybridization energy. An example of E (atom, msp 2 ^ for E(c,2sp 3> )is given in Eq. (14.503) 
where the sum of the negative of the first ionization energy of C, -11.27671 eV , plus the 
hybridization energy to form the C2sp 3 shell given by Eq. (14.146) is 

E(c,2sp 3> )= -14.63489 eV . 

Thus, the sharing of electrons between two atom msp 3 HOs to form an atom-atom- 
bond MO permits each participating hybridized orbital to decrease in radius and energy. In 
order to further satisfy the potential, kinetic, and orbital energy relationships, each atom msp 3 
HO donates an excess of 25% per bond of its electron density to the atom-atom-bond MO to 
form an energy minimum wherein the atom-atom bond comprises one of a single, double, or 

> triple bond. In each case, the radius of the hybridized shell is calculated from the Coulombic 
energy equation by considering that the central field decreases by an integer for each 
successive electron of the shell and the total energy of the shell is equal to the total 
Coulombic energy of the initial AO electrons plus the hybridization energy. The total energy 

E T (nioLatom, msp 3 ) (m is the integer of the valence shell) of the HO electrons is given by 

> the sum of energies of successive ions of the atom over the n electrons comprising total 
electrons of the at least one initial AO shell and the hybridization energy: 

E T lmol.atom,msp 1= E(atom,msp 3 y^JP m (15.17) 

m=2 
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where IP m is the mth ionization energy (positive) of the atom and the sum of ~IP } plus the 
hybridization energy is E(atom,msp 3 ^. Thus, the radius r , of the hybridized shell is given 
by: 

Vg=z-/; J i$7V8 0 E T \moLatom 9 msp \ 

where j = 1,2,3 for a single, double, and triple bond, respectively. The Coulombic energy 

E couhmb (mol.atom,msp 3 ^ of the outer electron of the atom msp 3 shell is given by 

2 

E co U i omb (^ ol * tom > ms P 3 )= -~~ r — (15.19) 

0 msp* 

In the case that during hybridization at least one of the spin-paired AO electrons is unpaired 
in the hybridized orbital (HO), the energy change for the promotion to the unpaired state is 
the magnetic energy E(magnetic) at the initial radius rof the AO electron given by Eq. 

(15.15). Then, the energy E(moLatom,msp 3 ^ of the outer electron of the atom msp 3 shell is 
given by the sum of E CoulQmb (moLalom, msp 3 ^ and E(magnetic) : 

(\ —e 2 2.71U e 2 fi 2 
moLatom, msp 3 )= + ' ° (1 5 .20) 

u mspr e 

E T {atom- atom, msp 3 y the energy change of each atom msp 3 shell with the formation of the 
atom-atom-bond MO is given by the difference between E (nioLatom, msp 3 ^ and 
E(atom,msp 3 ^: 

E T (atom - atom, msp 3 ^= E (moLatom, msp 3 ^- E (atom, msp 3 ^ (15.21) 
As examples from prior sections, E Couhmb (moLatom, msp 3 ) is one of: 

E Coulomb (p ethylene > 2S p 3 )> ^Coulomb (f ethane > 2s P^ )> E Coulomb acetylene > 2iS p 3 )' 3Jld 

E Coulomb (fpafkzme ? 2s P 

. . E cou !om6 (?tom,m S p 3 ) is one of E^ ulomi (casp 3 )and E Cou]omi (ci,3sp 3 }, 
^ E(mol.atom,msp 2 ) is one of E(c ethylene ,2sp 3 ), E(C elhane ,2sp 3 ), 
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E (atom,, msp 3 ^ is one of and E (c, 2sp 3 ^ and E (ci, 3sp 3 ^; 
E T (atom - atom, msp 3 ) is one of e(c~- C,2$p 3 e(c = C, 2^ 3 ), and. e(c = C, 2,sp 3 

arrow? »25p 3 is one of C2sp 3 , Cl3sp 3 
E T (atom - atom (s x ^ msp 3 ^ is ^ - C,2sp 3 ^ and E T (atom - atom (s 2 ^ msp 3 ^ is 

E T (c = C,2sp*y and 

r 3 is one of r - , r , „ 3 , r , , 3 ? r , , r , 3 9 and r . . 

mspr Clsp* ethane2sp 3 ethylemlsp* acetylene2sp J p alkamlsp* Cllsp* 

In the case of the C2sp 3 HQ, the initial parameters (Eqs. (14.142-14.146)) are 

r ,=V -r x = 7 >: = 0.91771a rt (15.22) 

^ $7te Q (^148.25751 eV) 8^ 0 (el48.25751 eV) 0 

2 2 

£ CD1I , . (C,2?p 3> )= ~ 6 = — = -14.82575 eV (15.23) 

Coulomb^ V J g 8TO.0.91771O, V ' 

E(magnetic) = ^° = ^* =0.19086 (15.24) 

m 2 (r 3 J (o.84317«J 

2 2 W 2 h 2 



, =-14.82575 eF + 0.19086 eV (15.25) 

= -14.63489 eV 

In Eq. (15.18), 

(Z-<?) = 10 (15.26) 

?=Z-n 

Eqs. (14.147) and (15.17) give 

E T (mol.atom,msp 3> )= E T (C e , Aawe ,2sp 3 )= -151.61569 eV (15.27) 

Using Eqs. (15.18-15.28), the final values of r Q ^ , E Cmlomb (C2sp 3 ), and £(C2sp 3 ), and the 
resulting E r ^C-c,C2sp*j of the MO due to charge donation from the HO to the MO where 
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BO 



C - C refers to the bond order of the carbon-carbon bond for different values of the 
parameter s are given in Table 15,1. 

Table 15.1. The final values of r^, E Coulomb (cisp 3 ^ and E(C2sp 3 ^ and the resulting 

f BO \ B ° 

E T \^c B -C,C2sp 3 j °f the MO due to charge donation from the HO to the MO where C - C 
refers to the bond order of the carbon-carbon bond. 



MO 
Bond 
Order 
(BO) 


s 1 


s 2 


Final 


(eV) 
Final 


E(C2sp^ 

(eV) 
Final 


( BO \ 

E T \C - C.dsp 3 J 

(eV) 


I 


1 


0 


0.87495 


-15.55033 


-15.35946 


-0.72457 


n 


2 


0 


0.85252 


-15.95955 


-15.76868 


-1.13379 


in 


3 


0 


0.83008 


-16.39089 


-16.20002 


-1.56513 


IV 


4 


0 


0.80765 


-16.84619 


-16.65532 


-2.02043 



x E (mol.atom (s. \ wisp 3 V E (mol.atom fa \ msp 3 ) 
atom,msp 3 y-^— . Vl/ \ V W J ' (15.28) 



In another generalized case of the basis of forming a minimum-energy bond with the 
constraint that it must meet the energy matching condition for all MOs at all HOs or AOs, the 

energy E^moLatom, msp 3 ^ of the outer electron of the atom msp 3 shell of each bonding atom 
must be the average of E(niol.atom, msp 1 ^ for two different values of s : 

Elmol.L , 

^ / 2 

In this case, E T (atom - atom, msp 2 ^ 5 the energy change of each atom msp 3 shell with the 

formation of each atom-atom-bond MO, is average for two different values of s : 

s v E T (atom - atom (s, \ msp 3 Yh E rj , (atom - atom (s 0 \ msp 3 ) 

E T {atom - atom, msp 3 )= — y r , v ^ (1 5.29) 

2 

Consider an aromatic molecule such as benzene given in the Benzene Molecule 
section. Each C = C double bond comprises a linear combination of a factor of 0.75 of four 
paired electrons (three electrons) from two sets of two C2sp 3 HOs of the participating carbon 
atoms. Each C - H bond of CH having two spin-paired electrons, one from an initially 
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unpaired electron of the carbon atom and the other from the hydrogen atom, comprises the 
linear combination of 75% H 2 -type ellipsoidal MO and 25% C2sp 3 HO as given by Eq. 
(13.439). However, E T (atom - atom, msp" ) of the C-i^ -bond MO is given by 
0.5£ r (c = C,2ap 3 ) (Eq, (14.247)) corresponding to one half of a double bond that matches 

the condition for a single-bond order for C - H that is lowered in energy due to the aromatic 
character of the bond. 

A further general possibility is that a minimum-energy bond is achieved with 
satisfaction of the potential, kinetic, and orbital energy relationships by the formation of an 
MO comprising an allowed multiple of a linear combination of # 2 -type ellipsoidal MOs and 
corresponding HOs or AOs that contribute a corresponding allowed multiple (e.g. 0.5, 0.75, 
1) of the bond order given in Table 15.1. For example, the alkane MO given in the 
Continuous-Chain Alkanes section comprises a linear combination of factors of 0.5 of a 
single bond and 0.5 of a double bond. 

; Consider a first MO and its HOs comprising a linear combination of bond orders and 
a second MO that shares a HO with the first. In addition to the mutual HO, the second MO 
comprises another AO or HO having a single bond order or a mixed bond order. Then, in 
order for the two MOs to be energy matched, the bond order of the second MO and its HOs or 
its HO and AO is a linear combination of the terms corresponding to the bond order of the 
mutual HO and the bond order of the independent HO or AO. Then, in general, 
E T (atom~ atom.msp 3 ^ the energy change of each atom msp 3 shell with the formation of 

each atom-atom-bond MO, is a weighted linear sum for different values of s that matches the 
energy of the bonded MOs, HOs, and AOs : 

E T {atom- atom,msp 3 y=^c s E T (atom - atom (s^ msp 3 } (15.30) 
where c is the multiple of the BO of s„ . The radius r , of the atom msp 3 shell of each 
bonding atom is given by the Coulombic energy using the initial energy E Coalomi (atom,msp^ 

and E T (atom - atom,msp 3 ^ the energy change of each atom msp* shell with the formation 
of each atom-atom-bond MO: 

-e 2 

r ,= , 5 : : C15.31) 



msp* 



8xs 0 a 0 (E Couhnb (atom,msp i yt E T (atom - atom,msp 3 ^ 
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where E Couhmb (C2^ 3 )= -14.825751 eV . The Coulombic energy E CoulQmb (wtoLatom,msp 3 ^ 

of the outer electron of the atom msp 3 shell is given by Eq. (15.19). In the case that during 
hybridization, at least one of the spin-paired AO electrons is unpaired in the hybridized 
orbital (HO) 9 the energy change for the promotion to the unpaired state is the magnetic energy 
E{magnetic) (Eq. (15.15)) at the initial radius rof the AO electron. Then, the energy 

E(mol.atom,msp 3 ^ of the outer electron of the atom msp 1 shell is given by the sum of 

"^Coulomb (^^ • 

atom.msp 3 ^ and Eimagnetic) (Eq. (15.20)). E T (atom — atom, msp 3 ^ the 
energy change of each atom msp 3 shell with the formation of the atom-atom-bond MO is 
given by the difference between E (mohatom, msp 3 ) and E (atom, msp 3 ^ given by Eq. 

(15.21). Using Eq. (15.23) for E Coulomb (ca$p 2 )in Eq. (15.31), the single bond order energies 

given by Eqs. (15.18-15.27) and shown in Table 15. 1 ? and the linear combination energies 
(Eqs. (15.28-15.30)), the parameters of linear combinations of bond orders and linear 
combinations of mixed bond orders are given in Table 1 5.2. 

Table 15.2. The final values of (c2*/? 3 ), and E(c2sp 3 } and the resulting 

E r ^C~ C,C2sp*j of the MO comprising a linear combination of H 2 -type ellipsoidal MOs and 

corresponding HOs of single or mixed bond order where c Sn is the multiple of the bond order 
parameter E T (atom - atom (s n } msp 3 ) given in Table 15.1. 
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MO 
Bond Order 
(BO) 


s 1 


r* 


s 2 




s 3 


1/21 


1 


0.5 


0 


0 


0 


1/2II 


2 


0.5 


0 


0 


0 


I + 1/2II 


1 


0.5 


2 


0.25 


0 


1/2II + (I + II) 


2 


0.25 


1 


0.25 


2 


3/4II 


2 


0.75 


0 


0 


0 


I + II 


1 


0.5 


2 


0.5 


0 


i + m 


1 


0.5 


3 


0.5 


0 


I + IV 


1 


0.5 


4 


0.5 


0 


11 + in 


2 


0.5 


3 


0.5 


0 


n + rv 


2 


0.5 


4 


0.5 


0 


in+iv 


3 


0.5 


4 


0.5 


0 


rv + iv 


4 


0.5 


4 


0.5 


0 
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MO 
Bond Order 
(BO) 


c 


Final 


(eV) 
Final 


£(C2«y) 

(eV) 
Final 


(BO \ 

E T \ C-CClsp" 
T \ J 

(eV) 


1/21 


0 


0.89582 


-15.18804 


-14.99717 


-0.36228 


1/2II 


0 


0.88392 


-15.39265 


-15.20178 


-0.56689 


I + 1/2II 


0 


0.87941 


-15.47149 


-15.28062 


-0.64573 


i/2ii+(i + n) 


0.25 


0.87363 


-15.57379 


-15.38293 


-0.74804 


3/4II 


0 


0.86793 


-15.67610 


-15.48523 


-0.85034 


I + II 


0 


0.86359 


-15.75493 


-15.56407 


-0.92918 


i + m 


0 


0.85193 


-15.97060 


-15.77974 


-1.14485 


i+rv 


0 


0.83995 


-16.19826 


-16.00739 


-1.37250 


n + in 


0 


0.84115 


-16.17521 


-15.98435 


-1.34946 


H + rv 


0 


0.82948 


-16.40286 


-16.21200 


-1.57711 


III + IV 


0 


0.81871 


-16.61853 


-16.42767 


-1.79278 


IV + IV 


0 


0.80765 


-16.84619 


-16.65532 


-2.02043 



Consider next the radius of the AO or HO due to the contribution of charge to more 
than one bond. The energy contribution due to the charge donation at each atom such as 
carbon superimposes linearly. In general, the radius r 3 of the C2sp 3 HO of a carbon 

atom of a given molecule is calculated using Eq> (14.514) by considering y\E (MO,2sp 3 X 

the total energy donation to each bond with which it participates in bonding. The general 
equation for the radius is, given by 

= 

W "8« 0 (^(c.2^>25 r f (jffl,V)) 

e 2 " (15.32) 

8ot 0 ^14.825751 eV + Y\ E r mol (MO,2sp 3 

The Coulombic energy ' E Coulomb (mol.atom.msp 3 ^ of the outer electron of the atom msp 3 shell 
is given by Eq. (15.19). In the case that during hybridization, at least one of the spin-paired 
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AO electrons is unpaired in the hybridized orbital (HO), the energy change for the promotion 
to the unpaired state is the magnetic energy Eimagnetic) (Eq. (15.15)) at the initial radius 

rof the AO electron. Then, the energy E (moLatom, msp 3 ) of the outer electron of the 

atom msp 3 shell is given by the sum of E Coulomb (moLatom,msp^ and E (magnetic) (Eq. 
(15.20)). 

For example, the C2$p s HO of each methyl group of an alkane contributes 

-0.92918 eV (Eq. (14.513)) to the corresponding single C-C bond; thus, the 

corresponding C2sp 3 HO radius is given by Eq. (14.514). The C2sp 3 HO of each methylene 

group of C n H 2n+2 contributes -0.92918 eV to each of the two corresponding C-C bond 

MOs. Thus, the radius (Eq. (15.32)), the Coulombic energy (Eq. (15.19)), and the energy 
(Eq. (15.20)) of each alkane methylene group are 



-e 2 

y 



. 8 ^o (Be*** ( C 'V)^E E t„^ (methylene C - C,2sp 3 )) 



e 2 



8>r£ 0 (el4.825751 eF + e0.92918 eF + e0.92918 eV^ 
= 0.81549« 0 



(15.33) 



^Coulomb 



2 

( C n,e t Me n e 2s P 3> h — "T^ -v = -16.68412 (15,34) 

V methylene JfJ. ( 0 g ! 549^ ) ^ ■> 

^/O, ' ,\ -e 1 2n;u n e 2 h 2 

E \C melhytsne 2sp 3 )^- , Q . + - = -16.49325 eV (15.35) 

V 7 8^ 0 (0.81549« 0 ) OTe 2 (o.84317aJ 

In the determination of the parameters of functional groups, heteroatoms bonding to 
C2sp 3 HOs to form MOs are energy matched to the C2sp 3 HOs. Thus, the radius and the 
energy parameters of a bonding heteroatom are given by the same equations as those for 
C2sp 3 HOs. Using Eqs. (15.15), (15.19-15.20), (15.24), and (15.32) in a generalized fashion, 
the final values of the radius of the HO or AO, r AtomHOAO , E Coulomb (mol.atom,msp 3 ^, and 
E (C mol 2sp 3 ) are calculated using V E T (mo, 2sp 3 \ the total energy donation to each bond 

^ S ' group \ S 

with which an atom participates in bonding corresponding to the values of e t \C - C 9 C2sp^j °f 
the MO due to charge donation from the AO or HO to the MO given in Tables 15.1 and 15.2. 
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Table 15.3.A. The final values of r AtomJfOAO , E Couhmb (moLatom 9 msp 3 y and E(C mal 2sp 3 ) 
calculated using the values of eIc- C,C2sp*\ given in Tables 15.1 and 15.2. 
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Table 15.3.B. The final values of r AtomHOAO , E Couhmb (moLatom y msp 3 y and E[C mol 2sp 3 ) 
calculated for heterocyclic groups using the values of eJ\C B - ' C 9 C2sp*j gi yen m Tables 15.1 
and 15.2. 
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The energy of the MO is matched to each of the participating outermost atomic or 
hybridized orbitals of the bonding atoms wherein the energy match includes the energy 
contribution due to the AO or HO's donation of charge to the MO. The force constant k' 
(Eq. (15.1)) is used to determine the ellipsoidal parameter c 1 (Eq. (15.2)) of the each H 2 - 
type-ellipsoidal-MO in terms of the central force of the foci. Then, c' is substituted into the 
energy equation (from Eq. (15.11))) which is set equal to n x times the total energy of H 2 
where n l is the number of equivalent bonds of the MO and the energy of H 2 , 
-31.63536831 eV , Eq. (11.212) is the minimum energy possible for a prolate spheroidal 
MO, From the energy equation and the relationship between the axes, the dimensions of the 
MO are solved. The energy equation has the semimajor axis a as it only parameter. The 
solution of the semimajor axis a then allows for the solution of the other axes of each prolate 
spheroid and eccentricity of each MO (Eqs. (15.3-15.5)). The parameter solutions then allow 
for the component and total energies of the MO to be determined. 

The total energy, E T (p 2 Mo\ is given by the sum of the energy terms (Eqs. (15.6- 

15.11)) plus E T (AO/HO): 

E T (H 2 Mo)=V e +T + V m +V p +E T (AO/ HO) (15.36) 



E T (fi 2 Mo)= - 



n i e 



C l C 2 



^ a ) a-^a 2 -b 2 



+ e t (aoi ho) 



~„ a + c „ 
c,cj 2 — s. |in 1 



(15.37) 



+ 



e t (aoi HO) 



where «, is the number of equivalent bonds of the MO, c x is the fraction of the H 2 -type 
ellipsoidal MO basis function of a chemical bond of the group, c 2 is the factor that results in 
an equipotential energy match of the participating at least two atomic orbitals of each 
chemical bond, and E T (AO I HO) is the total energy comprising the difference of the energy 

e(ao/ho) of at least one atomic or hybrid orbital to which the MO is energy matched and 
any energy component AE ff ^ MO (AO / HO) due to the AO or HO's charge donation to the 
MO. 

E T (AO/HO)^ e(A0/H0)- AE ff2MO (AO/Ho) (15.38) 
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As specific examples given in previous sections, E T (AO I HO>) is one from the group of 



e t (ao/ho 

E t (A0/ HO 

e t (ao/ho 
e t (ao/ho 
e t (ao/ho 
e t (aoi ho 
e t (ao/ho 

E T (AO I HO 

e t (ao/ HO 
e t (aoi ho 
e t (ao/ho 
e t (ao/ho 

E T (AOI HO 



= E(02p shell)= -E(jonization; O) = -13.6181 eV ; 
i= E(N2p shell)= -Eijonization; N) = -14.53414 eV ; 

E(c,2sp 3S )= -14.63489 eV ; 
= E CoM p>3sp>)= -14.60295 eV ; 

E(ionization\ C) + E (ionization; C + ) ; 
^E(C elhane ,2sp 3 )= -15.35946 eV ; 
\=+E (C elh}ilenii , 2sp 3 )— E ({C etfvlene , 2sp 3 ^; 

= E(c,2sp 3 y-2E r (C = C,2jp 3 )= -14.63489 eV - (-2.26758 eV); 
= E(C acelylene ,2sp 3 y E(C a ce(ylsne ,2sp 3 }- E(C acetylene ,2sp 3 )= 16.20002 eV : 
i= E(c,2sp 3 y 2E T (c = C,2sp 3 )= -14.63489 eV - (-3.13026 eV); 

= E(c,2sp 3 )- E t (c = C,2sp 3 )= -14.63489 eV - (-1.13379 eV), and 
= E(C alkane ,2 S p 3 )=~ 15.56407 eV . 



To solve the bond parameters and energies, c } = a. 



0 _ 



m e 2 2QC 0 a 

e 12 



' 2C 1 C 2 



(Eq. 



(15.2)) is substituted into E T (h 2 mo) to give 



E T (h 2 mo)= — 



n,e 



a o1i a + yja 2 -b 2 



c,cj 2 2-1 In " ' 

^ <*) a- 4 a 2 - 



"12 



-1 



+ e t (ao/ho) 



n x e 



\ a J a- 



+ e t (aoi ho) 



(15.39) 



C a 0 V V 2CC 2 
^[2-^jln V^- 1 



\2C X C 2 



+ e t (aoiho) 



The total energy is set equal to E(basis energies) which in the most general case is given by 
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the sum of a first integer n x times the total energy of H 2 minus a second integer n 2 times the 
total energy of H 9 minus a third integer n 3 times the valence energy of E(AO^ (e.g. 

e(n^= -14.53414 eV) where the first integer can be 1.2,3... ? and each of the second and 
third integers can be 0,1,2, 3... . 

E(basis energies) = n x (-31.63536831 eV)- n 2 (-13.605804 eV)- n 2 E(AO) (15.40) 
In the case that the MO bonds two atoms other than hydrogen, E(basis energies) is n x times 
the total energy of H 2 where n x is the number of equivalent bonds of the MO and the energy 
of H 2 , -31.63536831 eV , Eq. (11.212) is the minimum energy possible for a prolate 
spheroidal MO: 

E(basis energies) = n x (-31.63536831 eV^ (15.41) 

E T (h 2 mo), is set equal to E(basis energies) , and the semimajor axis a is solved. 
Thus, the semimajor axis a is solved from the equation of the form: 



n x e 2 



2>7V€ t ' °~ 



aa Q 

a+ 1 - 



„ a 0 V V 2C i C 2 



pCf 2 



+ E T (AO/ HO)= E(basis energies) (15.42) 



The distance from the origin of the H 2 -type-ellipsoidal-MO to each focus c' , the internuclear 
distance 2c 1 , and the length of the semiminor axis of the prolate spheroidal H 2 -type MO 
b = c are solved from the semimajor axis a using Eqs. (15.2-15.4). Then, the component 
energies are given by Eqs. (15.6-15.9) and (15.39). 

The total energy of the MO of the functional group, E T (k>), is the sum of the total 
energy of the components comprising the energy contribution of the MO formed between the 
participating atoms and E T (atom - atom,msp 3 .AO}, the change in the energy of the AOs or 
HOs upon forming the bond. From Eqs. (15.39-15.40), E T (mo) is 

E (mo)= E(b asis energies) + E T (atom — atom, msp 3 .AO^ (15.43) 
During bond formation, the electrons undergo a reentrant oscillatory orbit with 
vibration of the nuclei, and the corresponding energy E osc is the sum of the Doppler 5 E , and 
average vibrational kinetic energies, E K b : 
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(15.44) 



\mc" z \ ju 

where n x is the number of equivalent bonds of the MO , k is the spring constant of the 
equivalent harmonic oscillator, and fx is the reduced mass. The angular frequency of the 
reentrant oscillation in the transition state corresponding to E D is determined by the force 
between the central field and the electrons in the transition state. The force and its derivative 
are given by 



BO 



C, CL e 2 
4xe 0 R 3 



and 



r(a)=2c B 



BO 



(15.45) 



(15.46) 



such that the angular frequency of the oscillation in the transition state is given by 



i 



m 



C. C, e l 

4xe 0 R 3 



m 



(15.47) 



where R is the semimajor axis a or the semiminor axis b depending on the eccentricity of 
the bond that is most representative of the oscillation in the transition state, c B0 is the bond- 
order factor which is 1 for a single bond and when the MO comprises n x equivalent single 
bonds as in the case of functional groups. c B0 is 4 for an independent double bond as in the 
case of the C0 2 and N0 2 molecules and 9 for an independent triplet bond. C lo is the 
fraction of the H 2 -type ellipsoidal MO basis function of the oscillatory transition state of a 
chemical bond of the group, and C 2o is the factor that results in an equipotential energy match 
of the participating at least two atomic orbitals of the transition state of the chemical bond. 
Typically, C Xo = C x and C 2o = C 2 . The kinetic energy, E K 9 corresponding to E D is given by 
Planck's equation for functional groups: 



E K = hco = , 



lo 2o 



1 



4xe Q R 3 



(15.48) 



m 



The Doppler energy of the electrons of the reentrant orbit is 
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\2E K 



= E ^ 



2h^ 


C C e 2 


4ns 0 R 3 




m c 2 



E osc given by the sum of E D and E Kvlb is 



( 



hv 



2h^ 


C, C, e 2 


4tvs 0 R 3 




„„ 2 



(15.49) 



+ E. 



vib 



(15.50) 



E of a group having bonds is given by E T (mo)/ n x such that 



r 



2K \l k 



E T+osc (group) is given by the sum of E T (a/o) (Eq. (15.42)) and E^ (Eq. (15.51)): 



(( 



aa n 



Cl c 2 [2-^ln- 



2££l_i 



\2Cf 2 



^+E T (AO I ifO)+ ^ (atom - atom,msp\AO^ 



1+ | 



2^ 




Arcs R 3 

o 




mc 2 



1. * 

'2 V 



energies 



) + E r(c 



atom - a/<wz, /wsp* 



(15.51) 
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(15.52) 

The total energy of the functional group E T Qnmp) is the sum of the total energy of the 
components comprising the energy contribution of the MO formed between the participating 
atoms, E(fiasis energies) , the change in the energy of the AOs or HOs upon forming the 

bond ( E T (atom - atom, msp 3 .AO^), the energy of oscillation in the transition state, and the 

change in magnetic energy with bond formation, E mag . From Eq. (15.52), the total energy of 

the group E T is 



( 



E T (group)— 



(E(basis energies) + E T (atom - atom.msp 3 ,AO^ 1 + ^ 



2ft 



C C e 



1 



Ans R 3 

O 



m c 



\ 1 Kvib mag 



(15.53) 



The change in magnetic energy E mag which arises due to the formation of unpaired electrons 
in the corresponding fragments relative to the bonded group is given by 

(15.54) 



2 3 

m r 



where r 3 is the radius of the atom that reacts to form the bond and c 3 is the number of 
electron pairs. ; 



( 



E T (group) - 



(E(basis energies) + E T (atom - atom.msp 3 AO^ 



1 + 



2ft^ 



AxsR 3 



m 



m c 



(15.55) 



The total bond energy of the group E D (group) is the negative difference of the total energy of 
the group (Eq. (15.55)) and the total energy of the starting species given by the sum of 

C A^ initial 

(c 4 AO I HO) and c 5 E lnWa , (c 5 AO I HO): 



35 



WO 2007/053486 



PCT/US2006/042140 



E D <£roup)= - 



a C 9 e 2 



Aks R 3 



m 



(e (basis energies) + E T (atom - atom.msp 3 ,AO^\ 1 + 1 



(15.56) 



In the case of organic molecules, the atoms of the functional groups are energy matched to the 
C2sp 3 HO such that 

E(AO/ HO)= -14.63489 eV (15.57) 
For examples of E from previous sections: 



p. = c, , S ""'^ ^ = £,0.14803 eV 



(0.91771oJ 

T n 



E ; (N2p)=c, - :" =c. 



8 w4 



mag 



^3 3 

r 



— { - s ^ — = c,0.14185 eV 



(15.58) 

(15.59) 
(15.60) 



(0.93084a 0 J 

In the general case of the solution of an organic functional group, the geometric bond 
parameters are solved from the semimajor axis and the relationships between the parameters 
by first using Eq. (15.42) to arrive at a. Then, the remaining parameters are determined 
using Eqs. (15.1-15.5). Next, the energies are given by Eqs. (15.52-15.59). To meet the 
equipotential condition for the union of the H 2 -type-ellipsoidal-MO and the HO or AO of the 
atom of a functional group, the factor c 2 of a i? 2 -type ellipsoidal MO in principal Eqs. 
(15.42) and (15.52) may given by 



(i) one: 

c 2 = l (15.61) 

' ■ 1 (ii) the ratio that is less than one of 13.605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1.243), and the magnitude of the 
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Coulombic energy of the participating AO or HO of the atom, E Couhmb (MO.atom,msp 3 ^ given 
by Eqs. (1 5 . 1 9) and (1 5 .3 1 - 1 5.32). For | E Coulomb (MO.atom, msp 3 ]> 13.605804 eV : 



e 2 



8^ 0 a 0 13.605804 eV 

— =17 tttt: n ( 15 - 62 ) 



_f! \ E couiomb (MO.atom, msp 3 ) 

% KS 0 r A-B AorBsp 3 

For \E Coulomb (MO. atom, tnsp 3 ^< 1 3.605804 eV : 



e 2 



_ ^^A-B AorBsp 3 \ E Coulomb (MO .CltOm,mSp 3 } 



13.605804 eV (15 ' 63) 



(iii) the ratio that is less than one of 13.605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1.243), and the magnitude of the 
valence energy, Eiyalence), of the participating AO or HO of the atom where Eiyalence) is 

the ionization energy or E (MO.atom, msp 3 ) given by Eqs. (15.20) and (15.31-15.32). For 

\E(yalence)\ >13.605804 eV : 



e 2 



8x£ 0 a 0 13.605804 eV 

e 2 \Eiyalence)\ (1564) 



0 A-B AorBsp 3 



For \E(yalence)\ < 13.605804 eV : 



e 2 



r %Tts Q r A -BAorBtf \E(yalence)\ 

2 ~ e 2 ~ 13.605804 eV (15,65) 

8x& Q a Q 



(iv) the ratio that is less than one of the magnitude of the Coulombic energy of the 
participating AO or HO of a first atom, E Coulomb (MO.atom 9 msp^ given by Eqs. (15.19) and 
(15.31-15.32), and the magnitude of the valence energy, E(yalence) , of the participating AO 
or HO of a second atom to which the first is energy matched where Eiyalence) is the 
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ionization energy or E^MO.atom ? msp 3 ^ given by Eqs. (15.20) and (15.31-15.32). For 
E couiomb fao. atom, msp 3 j > E{yalence) : 

c 2 = -. if*^ 1 3 , (15.66) 

\Ecouiomb\MO.atom,msp J| 



For 



Ecouiomb (MO. atom, msp 3 ^ < Eiyalence) : 



\Ecouiomb (MO. atom, msp ) 

c 2 = ! r-^ ; ^ (15.67) 

2 \E(yalence)\ v J 



(v) the ratio that is less than one of the magnitude of the valence-level energies, 
E n (yalence), of the AO or HO of the nth participating atom of two that are energy matched 

where Eiyalence) is the ionization energy or E(MO.atom,msp 1 '^ given by Eqs. (15.20) and 

(15.31-15.32): 

^ E^alence) 

E 2 (yalence) 

(vi) the factor that is the ratio of the hybridization factor c 2 (l) of the valence AO or 

HO of a first atom and the hybridization factor c 2 (2^ of the valence AO or HO of a second 

atom to which the first is energy matched where c 2 (n) is given by Eqs. (15.62-15.68); 

alternatively c 2 is the hybridization factor c 2 (ij of the valence AOs or HOs a first pair of 

atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 
pair to which the first two are energy matched: 



a 



(2) 



(15.69) 



(vii) the factor that is the product of the hybridization factor c 2 (l) of the valence AO 
or HO of a first atom and the hybridization factor c 2 (2) of the valence AO or HO of a second 
atom to which the first is energy matched where <? 2 (V) is given by Eqs. (15.62-15.69); 
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alternatively c 2 is the hybridization factor c 2 (l) of the valence AOs or HOs a first pair of 

atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 
pair to which the first two are energy matched: 

c 2 =c 2 (l)c 2 (2) (15.70) 
The hybridization factor c 2 corresponds to the force constant k (Eqs. (11.65) and (13.58)). 
In the case that the valence or Coulombic energy of the AO or HO is less than 13.605804 eV , 
the magnitude of the Coulombic energy between the electron and proton of H given by Eq. 
(1.243), then C 2 corresponding to (Eq. (15.1)) is given by Eqs. (15.62-15.70). 

Specific examples of the factors c 2 and C 2 of a H 2 -type ellipsoidal MO of Eq. 
(15.51) given in following sections are 

c t (C2s P >HO ,0 F y ^Zgpc, :%%g*> *1771)- 0.77087 , 
c, ,0 ay <c* sp >Hoy Zl||^ ( 0. 91771 ) = „. 8131 , , 

C, ( C 2^HO ,0 Br y 1 ^±y 2 (02^0)- =!i||I^(„. 91771 )= 0 . 74081 , 

^c^Bo^y^^ici^^y^^^rny^, 

o)= i |^ Cj (« O )=zl||i^ (0 . 9I771)=0 . 85 3 95; 

fn. ,oaA -14.53414 eV ' 

(H t0 1 E^)= -15.35946 eV = ^ '> 

c 2 (C2s^HO t o N)= ^0^c 2 (Gl^HOy > 9177l)- 0.91140 ; 

c 2 {H to -jf^L = -".33414 «F - 0.93383 ; 

V ' E(c,2sp 3 ) -15.56407 eV ' 

^ (m * tj\ E (S,3p) -10.36001 eV nrr ^ AA 
C 2 [S3p to H) = ^, \' = = 0.76144 ; 

V ' E(H) -13.60580 eV 

C 2 (C2sp>HO to S)= , E @> . c 2 (C2sp 3 HO)= - 10 - 36001 eF ( 0 .91771> 0.64965 ; 
K J E(C,2sp 3 ) V 7 -14.63489 eV K J 
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m 

-13.61806 eV . 



, (p to S3sp 3 to C2sp 3 HO)= -j^c 2 (C2sp 3 HO) 

: (0.91771); 



-10.36001 eV 
= 1.20632 



( 5)= -11.57099 *r . S5045 

v ' B(H) -13.60580 eV 



C, (Clsp'HO ,o SV> fr* > EiliUS^- 85045 )' 0 66951 » 

C 2 (S3sp 3 to O to C2sp 3 HO)= ~^^^ c 2 (C2sp i HO) 

- 11 - 5 ' 2126 * F (0.91771) 5 



-13.61806 eF 
= 0.77641 



c 2 (O to N2p to C2sp 3 HO)= (C2sp 3 HO) 



13.61806,K (091771) . 



-14.53414 eF 
= 0.85987 



x c 2 (C2sp 3 HO to N) 0.91140 

c 2 (N2p to 02p)= ^ — < = = 1.06727 ; 

c 2 (C25^ 3 iTOto O) 0.85395 

C 2 (benzeneC2sp 3 HO)= c 2 (benzeneC2sp 3 HO)= 13 - 60S8Q4 eV = 0 .85252 ; 
v y v 7 15.95955 eF 

c 2 (arylC2sp 2 HO to o)= 3^ 2 (arylC2sp 3 HO) 

^ 3 - 6I806eF (0.85252) ; 



-14.63489 eV 
0.79329 



/rr -E(A0 -14.53414 <?F 

c 2 (# to anline JV)= / v = 0.92171 ; 

2V ' E{C s 2sp 3 ) -15.76868 eV 

c 2 (arylC2sp 3 HO to n)= f^ N \\ c 2 (arylC2sp 3 HO) 



-14.53414 eF . n __x 

= (0.85252) , and 

-14.63489 eV K ' 

= 0.84665 
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n (q>x + 14, ~n 3^ E(S,3p) -10.36001 eV 

C 2 [S3p to aiyl-type C2sp 3 HO )=—r — = = 0.65700 . 

v J E\C,2sp ) -15.76868 eV 



MO INTERCEPT ANGLES AND DISTANCES 

Consider the general case of Eqs. (13.84-13.95) wherein the nucleus of a B atom and the 
nucleus of a A atom comprise the foci of each H 2 -type ellipsoidal MO of an A- B bond. 
The parameters of the point of intersection of each H 2 -type ellipsoidal MO and the A -atom 
AO are determined from the polar equation of the ellipse: 
1 + e 



r = j\ 



(15.71) 



°l + ecos<9' 

The radius of the A shell is r A 9 and the polar radial coordinate of the ellipse and the radius of 
the A shell are equal at the point of intersection such that 



1 + 



(15.72) 



1 + — cos<9' 
a 



The polar angle 9' at the intersection point is given by 
( ( . c< v 



0' = cos" 1 



1+- 

(a-c')— 3— I 



(15.73) 



Then, the angle 0 A AO theradial vector of the A AO makes with the internuclear axis is 

0aao = ™O°-0' (15.74) 
The distance from the point of intersection of the orbitals to the internuclear axis must be the 
same for both component 1 orbitals such that the angle (ot-0 u between the internuclear 

axis and the point of intersection of each H 2 -type ellipsoidal MO with the A radial vector 
obeys the following relationship: 

r A Sin0 AAO = b ^ 0 H 2 MO (15.75) 

such that 

. ^sinfi* 

— din -1 _2 A AO 



b 



(15.76) 
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The distance d H MO along the internuclear axis from the origin of H 2 -type ellipsoidal MO to 



the point of intersection of the orbitals is given by 

The distance d AAO along the internuclear axis from the origin of the A atom to the point of 
intersection of the orbitals is given by 

d =c x -d (15.78) 

Cl A AO U H, MO V 



BOND ANGLES 

Further consider an ACB MO comprising a linear combination ofC-A -bond and C- B- 
bond MOs where C is the general central atom. A bond is also possible between the A and 
B atoms of the C- A and C-B bonds. Such A- B bonding would decrease the C - A 
and C-B bond strengths since electron density would be shifted from the latter bonds to the 
former bond. Thus, the ZACB bond angle is determined by the condition that the total 
energy of the # 2 -type ellipsoidal MO between the terminal A and B atoms is zero. The 
force constant k' of a H 2 -type ellipsoidal MO due to the equivalent of two point charges of 
at the foci is given by: 

7 -,_ C i C 2 2g2 (15.79) 
4xe o 

where is the fraction of the H 2 -type ellipsoidal MO basis function of a chemical bond of 
the molecule which is 0.75 (Eq. (13.59)) for a terminal A - H ( A is H or other atom) and 1 
otherwise and C 2 is the factor that results in an equipotential energy match of the 
participating at least two atomic orbitals of the chemical bond and is equal to the 
corresponding factor of Eqs. (15.42) and (15.52). The distance from the origin of the MO to 
each focus c y of the A- B ellipsoidal MO is given by: 



a h247rs o = aa » (15.80) 
\m/2C x C 2 a \2Cf 2 



The internuclear distance is 



2c < = 2 l-^ (15-81) 



2C l C 2 
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The length of the semiminor axis of the prolate spheroidal A- B MO b = c is given by Eq. 
(15.4). 

The component energies and the total energy, E T (k 2 mo), of the A- B bond are given 
by the energy equations (Eqs. (11.207-11.212), (11.213-11.217), and (11.239)) of H 2 except 
that the terms based on charge are multiplied by c BO , the bond-order factor which is 1 for a 
single bond and when the MO comprises n x equivalent single bonds as in the case of 
functional groups. c B0 is 4 for an independent double bond as in the case of the C0 2 and 
N0 2 molecules. The kinetic energy term is multiplied by c x BO which is 1 for a single bond, 2 
for a double bond, and 9/2 for a triple bond. The electron energy terms are multiplied by c x , 
the fraction of the H 2 -type ellipsoidal MO basis function of a terminal chemical bond which 
is 0.75 (Eq. (13.233)) for a terminal A - H {A is H or other atom) and 1 otherwise. The 
electron energy terms are further multiplied by c f 2 , the hybridization or energy-matching 
factor that results in an equipotential energy match of the participating at least two atomic 
orbitals of each terminal bond. Furthermore, when A - B comprises atoms other than H, 

E T (atom - atom.msp 3 AO^, the energy component due to the AO or HO's charge donation to 
the terminal MO, is added to the other energy terms to give E T (h 2 mo): 



E (h 2 mo)= ; 



,+ E T (atom - atom 9 msp 3 .AO^ 



(15.82) 

The radiation reaction force in the case of the vibration of A - B in the transition state 
corresponds to the Doppler energy, E D , given by Eq. (11.181) that is dependent on the 
motion of the electrons and the nuclei. The total energy that includes the radiation reaction of 
the ^4- B MO is given by the sum of E t (h 2 mo) (Eq. (15.82)) and E osc given Eqs. (11.213- 

11.220), (11.231-11.236), and (11.239-11.240). Thus, the total energy E T (A- B) of the 
A - B MO , including the Doppler term is 
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tf n \ a-hc' 



c,c'\ 2c„-c' -9- 1 In 




—J In - 1 -f E T [atom - msp 3 .^OJJ 



(15.83) 

where C lo is the fraction of the i/ 2 -type ellipsoidal MO basis function of the oscillatory 
transition state of the A- B bond which is 0.75 (Eq. (13.233)) in the case of H bonding to a 
central atom and 1 otherwise, C 2o is the factor that results in an equipotential energy match of 
the participating at least two atomic orbitals of the transition state of the chemical bond, and 



is the reduced mass of the nuclei given by Eq. (11.1 54). To match the boundary 



condition that the total energy of the A- B ellipsoidal MO is zero, E T (a - given by Eq. 
(15.83) is set equal to zero. Substitution of Eq. (15.81) into Eq. (15.83) gives 



0 = 



—e 






aa 0 


2Cf 2 




\ a % 


2Cf 2 



-1 



+ 



E T (atom ^ atom, msp* .AO^ 



C BO e 



aa n 



. ' (15.84) 

The vibrational energy-term of Eq. (15.84) is determined by the forces between the central 
field and the electrons and those between the nuclei (Eqs. (11.141-11.145)). The electron- 
central-field force and its derivative are given by 



/(«>-< 



Cl c' 2 e 2 
BO 4xs 0 a 3 



(15.85) 



and 
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The nuclear repulsion force and its derivative are given by 

7. 

/(a + <?')=■ 



8fte o (a + c'^ 



(15.86) 



(15.87) 



and 



/'(a + <?')= - 



47re o (a + c'^ 

such that the angular frequency of the oscillation is given by 



(15.88) 



-3 
a 



1 



1 4fts 0 a 3 



(15.89) 



Since both terms of E osc = E D + E Kvib are small due to the large values of a and c\ to very 
good approximation, a convenient form of Eq. (15.84) which is evaluated to determine the 
bond angles of functional groups is given by 
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8^ h 
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(15,90) 



From the energy relationship given by Eq. (15.90) and the relationship between the axes 
given by Eqs. (15.2-15.5), the dimensions of the A-B MO can be solved. The most 
convenient way to solve Eq. (1 5.90) is by the reiterative technique using a computer. 

A factor c 2 of a given atom in the determination of c' 2 for calculating the zero of the 
total A - B bond energy is typically given by Eqs. (15.62-15.65). In the case of a H- H 
terminal bond of an alkyl or alkenyl group, c r 2 is typically the ratio of c 2 of Eq. (15.62) for 
the H - H bond which is one and c 2 of the carbon of the corresponding C-H bond: 
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c 9 1 ^ 13.605804 eV 

° 2 c 2 (C2sp') E Coulomb (C-HC2sp 3 ) ( " } 

In the case of the determination of the bond angle of the ACH MO comprising a linear 
combination of C - A -bond and C-H -bond MOs where A and C are general, C is the 
central atom, and c 2 for an atom is given by Eqs. (15.62-15.70), c f 2 of the A~ H terminal 
bond is typically the ratio of c 2 of the A atom for the A- H terminal bond and c 2 of the C 
atom of the corresponding C - H bond: 

, _ c 2 (A(A-H)msp 3 ) 
C >-c 2 (C(C~H)(msp>) (15 - 92) 

In the case of the determination of the bond angle of the COH MO of an alcohol comprising 
a linear combination of C - O-bond and O-H -bond MOs where C , O, and H are carbon, 
oxygen, and hydrogen, respectively, c' 2 of the C-H terminal bond is typically 0.91771 
since the oxygen and hydrogen atoms are at the Coulomb potential of a proton and an electron 
(Eqs. (L236) and (10.162), respectively) that is energy matched to the C2sp 3 HO, 

In the determination of the hybridization factor c' 2 of Eq. (15.90) from Eqs. (15.62- 
15.70), the Coulombic energy, E Coulomb (AdO.atom, msp 2 ^ or the energy, E (MO.atom, msp 3 
the radius r R A w 3 of the A or B AO or HO of the heteroatom of the A- B terminal bond 

A-B AorBsp 

MO such as the C2sp 3 HO of a terminal C-C bond is calculated using Eq. (15.32) by 

considering y\E T {MO,2sp 3 \ the total energy donation to each bond with which it 

participates in bonding jas it forms the terminal bond. The Coulombic energy 

E couhmb (MO. atom.msp 3 ^ of the outer electron of the atom wisp* shell is given by Eq. (15.19). 

In the case that during hybridization, at least one of the spin-paired AO electrons is unpaired 
in the hybridized orbital (HO), the energy change for the promotion to the unpaired state is 
the magnetic energy E(magnetic) (Eq. (15.15)) at the initial radius r of the AO electron. 

Then, the energy E (hdO.atom, msp 3 ) of the outer electron of the atom msp 3 shell is given by 

the sum of E Coufomb (MO.atom 9 msp 3 ) and E(jnagnetic) (Eq. (15.20)). 

In the specific case of the terminal bonding of two carbon atoms, the c 2 factor of each 
carbon given by Eq. (15.62) is determined using the Coulombic energy 
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E (C-C C2sp 3 } of the outer electron of the C2sp 3 shell given by Eq. (15.19) with the 

Coulomb \ J 

radius r , of each C2sp 3 HO of the terminal C-C bond calculated using Eq. (15.32) 

C— C Clsp 

by considering V£ T (MO,2sp 3 \ the total energy donation to each bond with which it 

' mot \ / 

participates in bonding as it forms the terminal bond including the contribution of the 
methylene energy, 0.92918 eV (Eq. (14.513)), corresponding to the terminal C-C bond. 
The corresponding E T (atom - atom, msp 3 .AO^m Eq. (15.90) is 

E T (c~C C2sp 3 )= -1.85836 eV . 

In the case that the terminal atoms are carbon or other heteroatoms, the terminal bond 
comprises a linear combination of the HOs or AOs; thus, c' 2 is the average of the 
hybridization factors of the participating atoms corresponding to the normalized linear sum: 

c' 2 =- (c' 2 (atom l)+ c' 2 (atom 2)) (1 5.93) 

In the exemplary cases of C-C, O-O, and N-N where C is carbon: 

( e 2 e 2 



2 2 



-A 
~ 2 



8^g 0 a 0 t 8xs 0 a 0 



e 2 e 2 



\^ nG ^ r A-A\AOIHO ■ ^^A-AA^OIHO ) 

\ 13.605804 eV 13.605804 eV 

I E Cmlomb (A - A.A.AO I Hdf E Couhmb (a - A.A.AO I HO)) 



(15.94) 



(15.95) 



In the exemplary cases of C-N , C-O, and C-S , 
J 13.605804 eV (c „ A 

where C is carbon and c 2 (c to B^ is the hybridization factor of Eqs. (15.52) and (15.84) that 
matches the energy of the atom B to that of the atom C in the group. For these cases, the 
corresponding 1 E T (atom- atom,msp\AO) term in Eq. (15.90) depends on the hybridization 
and bond order of the terminal atoms in the molecule, but typical values matching those used 
in the determination of the bond energies (Eq. (15.56)) are 

E T (C - OC2sp i D2 p)= -1.44915 eV ; E T (c-0 C2sp 3 .02p)= -1.65376 eV ; 
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E T {C-N C2sp\N2 P y -1.44915 eV ; E T (c-S C2sp\S2p)= -0.72457 eV ; 

E T (0-0 02p.02p)= -1.44915 eV ; E r(P~° 02p.02p)= -1.65376 eV ; 

E T (N-N N2 p.N2p*)= -1 .4491 5 eV ; E t (n~ON2 pD2 p)= -1 .449 1 5 eV ; 

E T (F-F F2p.F2p)=~lA49l5 eV ; £ r (c/- C/ Cl3p.Cl3p)= -0.92918 ; 

£ r (Br - Br Br4p.Br4p)= -0.92918 eV ; E T (l - I I5p.I5p)= -0.36229 eV ; 

E T (C-F C2sp\F2p)= -1.85836 eV; E T (c~Cl C2sp\Cl3p)= -0.92918 eV ; 

E T (c~Br C2sp\Br4p)= -0.72457 eV ; E T (C- / C2,s/? 3 ./5 J p)= -0.36228 , and 

E T (p-Cl G2p.Cl3p)= -0.92918 eV . 

In the case that the terminal bond is X-X where X is a halogen atom, c 1 is one, 
and c f 2 is the average (Eq. (15.93)) of the hybridization factors of the participating halogen 
atoms given by Eqs. (15.62-15.63) where E^^QkfO.atom^msp 3 ^ is determined using Eq. 
(15.32) and E Couhmb (MO.atojn,msp 3 ^= 13.605804 for X=I. The factor C x of Eq. 

(15.90) is one for all halogen atoms. The factor C 2 of fluorine is one since it is the only 

halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude of 
the Coulombic energy between the electron and proton of H given by Eq. (1.243). For each 
of the other halogens, C7, Br , and I, C 2 is the hybridization factor of Eq. (15.52) given by 

Eq. (15.70) with c 2 (l) being that of the halogen given by Eq. (15.68) that matches the 

valence energy of X (E l (valence)}) to that of the C2sp 3 HO ( E 2 (valence) = -14.63489 eV 9 

Eq. (15.25)) and to the hybridization of C2sp 3 HO (c 2 (2)== 0.91771, Eq. (13.430)). 

E T (atom- atom^msp 3 \AO^ of Eq. (15.90) is the maximum for the participating atoms which 

is -1.44915 eF, -0.92918eF, -0.92918 eV , and -0.33582eF for F, CI, £r,and /, 
respectively. 

Consider the case that the terminal bond is C - X where C is a carbon atom and X 
is a halogen atom. The factors c x and C x of Eq. (15.90) are one for all halogen atoms. For 

X = F 9 c r 2 is the average (Eq. (15.95)) of the hybridization factors of the participating carbon 

and F atoms where c 2 for carbon is given by Eq. (15.62) and c 2 for fluorine matched to 
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carbon is given by Eq. (15.70) with c 2 (l) for the fluorine atom given by Eq. (15.68) that 
matches the valence energy of F (E^alence)^ -17.42282 eV) to that of the C2sp 3 HO 
(E 2 (yalence)= -14.63489 eV , Eq. (15.25)) and to the hybridization of C2sp 3 HO 
(c 2 (2)= 0.91771 , Eq. (13.430)). The factor C 2 of fluorine is one since it is the only halogen 
wherein the ionization energy is greater than that 13.605804 eV , the magnitude of the 
Coulombic energy between the electron and proton of H given by Eq. (1.243). For each of 
the other halogens, CI, Br, and /, ^ is the hybridization factor of the participating carbon 
atom since the halogen atom is energy matched to the carbon atom. C 2 of the terminal-atom 
bond matches that used to determine the energies of the corresponding C-X-bond MO. 
Then, C 2 is the hybridization factor of Eq. (15.52) given by Eq. (15.70) with c 2 (l) for the 

halogen atom given by Eq. (15.68) that matches the valence energy of X ( E^alence)) to 
that of the C2sp 3 HO (E 2 (valence)^ -14.63489 eV , Eq. (15.25)) and to the hybridization of 
C2sp 3 HO (c 2 (2)- 0.91771, Eq. (13.430)). E T (atom - atomjnsp\A<j) of Eq. (15.90) is 

the maximum for the participating atoms which is -1.85836 eV , -0.92918 eV , 
-0.72457 eV , and -0.33582 eV for F, CI, Br , and /, respectively. 

Consider the case that the terminal bond is H - X corresponding to the angle of the 
atoms HCX where C is a carbon atom and X is a halogen atom. The factors c x and C 1 of 
Eq. (15.90) are 0.75 for all halogen atoms. For X = F , c' 2 is given by Eq. (15.69) with c 2 
of the participating carbon and F atoms given by Eq. (15.62) and Eq. (15.65), respectively. 
The factor C 2 of fluorine is one since it is the only halogen wherein the ionization energy is 
greater than that 13.605804 eV , the magnitude of the Coulombic energy between the electron 
and proton of H given by Eq. (1.243). For each of the other halogens, CI, Br , and /, c r 2 is 
also given by Eq. (15.69) with c 2 of the participating carbon given by Eq. (15.62) and c 2 of 
the participating X atom given by c 2 = 0.91771 (Eq. (13.430)) since the X atom is energy 
matched to the C2sp 3 HO. In these cases, C 2 is given by Eq. (15.65) for the corresponding 
atom X where C 2 matches the energy of the atom X to that of H . 
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(15.98) 



Using the distance between the two atoms A and B of the general molecular 
group ACB when the total energy of the corresponding A- B MO is zero, the corresponding 
bond angle can be determined from the law of cosines: 

s x +*y 2 2 ~ 2s x s 2 cosme9 = s 3 2 (15.96) 

With s t = 2c' c _ A9 the internuclear distance of the C- A bond, s 2 = 2c' c _ B , the internuclear 

distance of each C- B bond, and s s = 2c\_ B 9 the internuclear distance of the two terminal 

atoms, the bond angle Q^ CB between the C - A and C-B bonds is given by 

(2c' c _J +(2^^) ^"2(2c^X 2c V 5 >sine^=(2 C vJ (15.97) 

[ 2(2 C ._X2^) . 

Consider the exemplary structure C b C a {O a )O b wherein C a is bound to C b9 O , and 0 5 . In 
the general case that the three bonds are coplanar and two of the angles are known, say 0 X and 
9 %9 then the third 0 3 can be determined geometrically: 

'•0 3 =-360-4 (15.99) 
In the general case that two of the three coplanar bonds are equivalent and one of the angles is 
known, say G x , then the second and third can be determined geometrically: 

(360- ft) 

0 2=Gi = K 2 J (15.100) 

ANGLES AND DISTANCES FOR AN MO THAT FORMS AN ISOSCELES 
TRIANGLE ( 

In the general case where the group comprises three A- B bonds having B as the central 
atom at the apex of a pyramidal structure formed by the three bonds with the A atoms at the 
base in the xy-plane. The C 3v axis centered on B is defined as the vertical or z-axis ? and any 
two A- B bonds form an isosceles triangle. Then, the angle of the bonds and the distances 
from and along the z-axis are determined from the geometrical relationships given by Eqs. 
(13.412-13.416): 

the distance d originB from the origin to the nucleus of a terminal B atom is given by 
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2C ' B - B (15.101) 



origin- B 2 sin 60° 

the height along the z-axis from the origin to the A nucleus d height is given by 

the angle O v of each A - B bond from the z-axis is given by 
9 v =tan- 1 f ^-^ (15.103) 

V ^height ) 

Consider the case where the central atom B is further bound to a fourth atom C and 
the B-C bond is along the z-axis. Then, the bond #^ 5C given by Eq. (14.206) is 

^c = 180 ™^ (15.104) 

DIHEDRAL ANGLE 

Consider the plane defined by a general ACA MO comprising a linear combination of two 
C - A -bond MOs where C is the central atom. The dihedral angle ^ ZBCMC4 between the 
ACA -plane and a line defined by a third bond with C , specifically that corresponding to a 
C - B -bond MO, is calculated from the bond angle and the distances between the A 9 

5, and C atoms. The distance ^ along the bisector of 9^ 4CA from C to the internuclear- 
distance line between A and A , 2c , is given by 

^-2c' c ^cos^^ . (15.105) 

where 2c y C _ A is the internuclear distance between A and C. The atoms A, A, and 5 
define the base of a pyramid. Then, the pyramidal angle can be solved .from the 

internuclear distances between A and A , 2c'^_^ 5 and between .4 and 5, 2^ ^ ? using the 
law of cosines (Eq. (15.98)): 

r (^vJ + (2c',,J^(2cyJ^ 



Q ZAEA = COS" 1 



(15.106) 



2(2c'^X 2c ^) 

Then, the distance <i 2 along the bisector of &^ BA from B to the internuclear-distance line 
2c* ,is given by 
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A-B 



cos- 



ts* 



ZABA 



(15.107) 



The lengths d { , d 2 , and 2c } C B define a triangle wherein the angle between d { and the 



internuclear distance between B and C , 2c' c is the dihedral angle 9 



'ZBCIACA 



solved using the law of cosines (Eq. (15.98)): 



0 ZBC/ACA= COS ~ 



that can be 



(15.108) 



SOLUTION OF GEOMETRICAL AND ENERGY PARAMETERS OF 
MAJOR FUNCTIONAL GROUPS AND CORRESPONDING ORGANIC 
MOLECULES 

The exemplary molecules given in the following sections were solved using the solutions of 
organic chemical functional groups as basis elements wherein the structures and energies 
where linearly added to achieve the molecular solutions. Each functional group can be treated 
as a building block to form any desired molecular solution from the corresponding linear 
combination. Each functional group element was solved using the atomic orbital and hybrid 
orbital spherical orbitsphere solutions bridged by molecular orbitals comprised of the H 2 - 

type prolate spheroidal solution given in the Nature of the Chemical Bond of Hydrogen-Type 
Molecules section. The energy of each MO was matched at the HO or AO by matching the 
hybridization and total energy of the MO to the AOs and HOs. The energy E (e.g. given 

by Eq. (15.58)) for a C2sp 3 HO and Eq.(15.59) for an Olp AO) was subtracted for each set 
of unpaired electrons created by bond breakage. 

The bond energy is not equal to the component energy of each bond as it exists in the 
molecule; although, they are close. The total energy of each group is its contribution to the 
total energy of the molecule as a whole. The determination of the bond energies for the 
creation of the separate parts must take into account the energy of the formation of any 
radicals and any redistribution of charge density within the pieces and the corresponding 
energy change with bond cleavage. Also, the vibrational energy in the transition state is 
dependent on the other groups that are bound to a given functional group. This will effect the 
functional-group energy. But, because the variations in the energy based on the balance of the 
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molecular composition are typically of the order of a few hundreds of electron volts at most, 
they were neglected. 

The energy of each functional-group MO bonding to a given carbon HO is 
independently matched to the HO by subtracting the contribution to the change in the energy 
of the HO from the total MO energy given by the sum of the MO contributions and 



(15,71-15.78) using the final radius of the HO of each atom. The final carbon-atom radius is 
determined using Eqs. (15.32) wherein the sum of the energy contributions of each atom to all 
the MOs in which it participates in bonding is determined. This final radius is used in Eqs. 
(15.19) and (15.20) to calculate the final valence energy of the HO of each atom at the 
corresponding final radius. The radius of any bonding heteroatom that contributes to a MO is 
calculated in the same manner, and the energy of its outermost shell is matched to that of the 
MO by the hybridization factor between the carbon-HO energy and the energy of the 
heteroatomic shell. The donation of electron density to the AOs and HOs reduces the energy. 
The donation of the electron density to the MO's at each AO or HO is that which causes the 
resulting energy to be divided equally between the participating AOs or HOs to achieve 
energy matching. 

The molecular solutions can be used to design synthetic pathways and predict product 
yields based on equilibrium constants calculated from the heats of formation. New stable 
compositions of matter can be predicted as well as the structures of combinatorial chemistry 
reactions. Further important pharmaceutical applications include the ability to graphically or 
computationally render the structures of drugs that permit the identification of the biologically 
active parts of the molecules to be identified from the common spatial charge-density 
functions of a series of active molecules. Drugs can be designed according to geometrical 
parameters and bonding interactions with the data of the structure of the active site of the 
drug. 

To calculate conformations, folding, and physical properties, the exact solutions of the 
charge distributions in any given molecule are used to calculate the fields, and from the fields, 
the interactions between groups of the same molecule or between groups on different 
molecules are calculated wherein the interactions are distance and relative orientation 
dependent. The fields and interactions can be determined using a finite-element-analysis 
approach of Maxwell's equations. 




The intercept angles are determined from Eqs. 
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AROMATIC AND HETEROCYCLIC COMPOUNDS 

Aromatic and heterocyclic molecules comprise at least one of an aromatic or a cyclic 
conjugated alkene functional group. The latter was described in the Cyclic and Conjugated 
Alkenes section. The aromatic bond is uniquely stable and requires the sharing of the electrons 
of multiple H 2 -type MOs. The results of the derivation of the parameters of the benzene 

molecule given in the Benzene Molecule (C 6 H 6 ) section can be generalized to any aromatic 
function group(s) of aromatic and heterocyclic compounds. 

C 6 H 6 can be considered a linear combination of three ethylene molecules wherein a 
C — H bond of each CH 2 group of H 2 C = CH 2 is replaced by a C-C bond to form a six- 
member ring of carbon atoms. The solution of the ethylene molecule is given in the Ethylene 
Molecule (CH 2 CH 2 ) section. The radius (0.85252a Q ) of the C2sp 3 shell of ethylene 

calculated from the Coulombic energy is given by Eq. (14.244). The Coulombic energy 
-^coulomb (£ethyiene^ s P 3 ^) (-15.95955 eV ) of the outer electron of the C2sp 3 shell is given by Eq. 
(14.245). The energy E(C ethylem ,2$p 3 ^ (-15.76868 eV ) of the outer electron of the C2sp 3 shell 
is given by Eq. (14.246). E T (c = C,2^p 3 ) (-1.13380 eV) (Eq. (14.247), the energy change of 
each C2$p 3 shell with the formation of the C = C-bond MO is given by the difference between 

C 6 H 6 can be solved using the same principles as those used 

to solve ethylene wherein the 2s and 2p shells of each C hybridize to form a single 2sp 3 shell 
as an energy minimum, and the sharing of electrons between two C2sp 3 hybridized orbitals 
(HOs) to form a molecular orbital (MO) permits each participating hybridized orbital to 
decrease in radius and energy. Each 2sp 3 HO of each carbon atom initially has four impaired 
electrons. Thus, the 6 H atomic orbitals (AOs) of benzene contribute six electrons and the six 
sp 3 -hybridized carbon atoms contribute twenty-four electrons to form six C - H bonds and six 
C = C bonds. Each C - H bond has two paired electrons with one donated from the H AO 
and the other from the C2sp 3 HO. Each C-C bond comprises a linear combination of a factor 

of 0.75 of four paired electrons (three electrons) from two sets of two C2sp 3 HOs of the 
participating carbon atoms. Each C - H and each C = C bond comprises a linear combination 
of one and two diatomic H 2 -type MOs developed in the Nature of the Chemical Bond of 
Hydrogen-Type Molecules and Molecular Ions section, respectively. 

Consider the case where three sets of C = C-bond MOs form bonds between the two 
carbon atoms of each molecule to form a six-member ring such that the six resulting bonds 
comprise eighteen paired electrons. Each bond comprises a linear combination of two MOs 
wherein each comprises two C2sp 3 HOs and 75% of a H 2 -type ellipsoidal MO divided between 
the C2sp 3 HOs: 

3(2 C2sp 3 + 0.75 H z MO}-* s(c= cX ethylene - type - bond MO 
( 3e } 

\-> 6 \C = C/— bond MO of benzene 



(15.142) 
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The linear combination of each # 2 -type ellipsoidal MO with each C2sp* HO further comprises 
an excess 25% charge-density contribution per bond from each C2sp 2 HO to the C = C -bond 
MO to achieve an energy minimum. Thus, the dimensional parameters of each bond C = C - 
bond are determined using Eqs. (15.42) and (15.1-15.5) in a form that are the same equations as 
those used to determine the same parameters of the C = C -bond MO of ethylene (Eqs. (14.242- 
14.268)) while matching the boundary conditions of the structure of benzene. 

Hybridization with 25% electron donation to each C-C -bond gives rise to the 

C w~2*P 3 HO-shell Coulombic energy E Coulomb (q_,25/? 3 ) given by Eq. (14.245). To meet 
the equipotential condition of the union of the six C2sp 3 HOs, c 2 and C 2 of Eq. (15.42) for the 

3e 

aromatic C-C -bond MO is given by Eq. (15.62) as the ratio of 15.95955 eV 9 the magnitude of 
E couio t „t,(Ci>eKene> 2s P 3 ) ( Ec l- (14.245)), and 13.605804 eV 9 the magnitude of the Coulombic 
energy between the electron and proton of H (Eq. (1 .243)): 

C 2 (henzeneC2sp 3 HO)= c 2 (benzeneClsp 3 HO)= 13 - 605804 eV = q.85252 (15.143) 

15.95955 

3e ' 

The energies of each C-C bond of benzene are also determined using the same 

equations as those of ethylene (Eqs. (14.251-14.253) and (14.319-14.333) with the parameters of 

■ 3e ■ 

benzene. Ethylene serves as a basis element for the C~C bonding of benzene wherein each of 

3e 

the six C=C bonds of benzene comprises (0.75)(4)= 3 electrons according to Eq. (15.142). 



The total energy of the bonds of the eighteen electrons of the C=C bonds of benzene, 
E t (c6 H 6> c=c )> is given by (6)(o.75> times E T+ao (jC = c) (Eq. (14.492)), the total energy of 
the C-C -bond MO of benzene including the Doppler term, minus eighteen times i?(c,2sp 3 ) 



3e 



(Eq. (14.146)), the initial energy of each C2sp 3 HO of each C that forms the C=C bonds of 



3e 



bond order two. Thus, the total energy of the six C=C bonds of benzene with three electron 
per aromatic bond given by Eq. (14.493) is 

E T [c 6 H 6 ,C=cj = (6X0 J5)E T+o , c (c « c)~ (6X3>(c ? 2^ 3 ) 

= (6)(0.75)(-66.05796eF)-18(-14.63489eF) (15.144) 
= -297.26081 eV - (-263.42798 eV) 
= -33.83284 eV 

The results of benzene can be generalized to the class of aromatic and heterocyclic compounds. 
E hv of an aromatic bond is given by E T (# 2 ) (Eqs. (1 1.212) and (14.486)), the maximum total 
energy of each H 2 -type MO such that 



E — n, \ E„ + E r ., )= n 

arc 1 \ D Kvib / 1 



- 31 .63536 8 31 C ^^ + I S jij (15.145, 

The factor of 0.75 corresponding to the three electrons per aromatic bond of bond order two 
given in the Benzene Molecule (C 6 H 6 ) section modifies Eqs. (15.52-15.56). Multiplication of 
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the total energy given by Eq. (15.55) by /, = 0,75 with the substitution of Eq. (15.145) gives the 
total energy of the aromatic bond: 



f E(basis energies) + E T (atom - atom,rnsp* .AO) 



E T (group)— f x 



-31.63536831 eV* 



2h\ 



Kc 2 y 

4ns R 3 

o 

m 



mc 



(15.146) 



The total bond energy of the aromatic group E D (g^) is the negative difference of the total 
energy of the group (Eq. (15.146)) and the total energy of the starting species given by the sum 
of o,E inml (c 4 AO I HO) and c^E mm (c 5 AO 1 HO); 

f ( E(hasis energies) + E T (atom- atom,msp 3 .AO) ^ 



E D <GraMp)= — 



-31.63536831 eV ] 



2h 



C, C, e 2 



m 



\ 



m c 



8Ws 



(15.147) 



" (f Aw fa I HO) + c 5 E tnillal (c 5 AO I HO)) 
Since there are three electrons per aromatic bond, c A is three times the number of aromatic 
bonds. 

Benzene can also be considered as comprising chemical bonds between six CH radicals 
wherein each radical; comprises a chemical bond between carbon and hydrogen atoms. The 
solution of the parameters of CH is given in the Hydrogen Carbide ( CH) section. Those of the 
benzene are given in the Benzene Molecule (C 6 H 6 ) section. The energy components of V e , V p9 

T, V m , and E T are the same as those of the hydrogen carbide radical, except that 
£ r (c==C,2^> 3 )= -1.13379 eV (Eq. (14.247)) is subtracted from E T (CH) of Eq. (13.495) to 
match the energy of each C-H-bond MO to the decrease in the energy of the corresponding 
C2sp 3 HO. In the corresponding generalization of the aromatic CH group, the geometrical 
parameters are determined using Eq. (15.42) and Eqs. (15.1-15.5) with 
E T (atom- atom 9 msp\AO)^ -1.13379 eV . 

The total energy of the benzene C - H -bond MO, E T [C - H\ 9 given by Eq. (14.467) 
is the sum of 0,5E T (c = C.lsp 3 ), the energy change of each C2$p* shell per single bond due to 

3e 

the decrease in radius with the formation of the corresponding C~C -bond MO (Eq. (14.247)), 
and E T (CH), the a MO contribution given by Eq. (14.441). In the corresponding 

benzene \ ✓ 

generalization of the aromatic CH group, the energy parameters are determined using Eqs. 

s \ -1 13379 eV 
(15.146-15.147) with j\ = 1 and E T [atom- atom,m$p 3 .AO — : — . Thus, the energy 
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3e 

contribution to the single aromatic CH bond is one half that of the C=C double bond 

3e 

contribution. This matches the energies of the CH and C=C aromatic groups, conserves the 
electron number with the equivalent charge density as that of s = lin Eqs. (15.18-15.21), and 

3e 

further gives a minimum energy for the molecule. Breakage of the aromatic C—C bonds to 
give CH groups creates unpaired electrons in these fragments that corresponds to c 3 = 1 in Eq. 
(15.56) with E mag given by Eq. (15.58). 

Each of the C-H bonds of benzene comprises two electrons according to Eq. (14.439). 
From the energy of each C-H bond, ( 12 CH) (Eq. (14.477)), the total energy of the 

twelve electrons of the six C — H bonds of benzene, E T (C 6 H 6 ,C - H^, given by Eq. (14.494) 
is 

E T (C 6 H 6 ,C-H)= (6)^E D ^ ( 12 C/f))- 6(-3.90454 eV)= -23.42724 eV (15.148) 

The total bond dissociation energy of benzene, E D (C 6 H 6 given by Eq. (14.495) is the negative 

( ' 3t O / x 

sum of E T ^C 6 H 6 ,C = CJ (Eq. (14.493)) and E T {C 6 H 6 ,C-H) (Eq. (14.494)): 

E D (C 6 H e y-^ 

= -((-33.83284 eV)+ (-23.42724 eV^j) (15.149) 
= 57.2601 eV 

Using the parameters given in Tables 15.214 and 15.216 in the general equations (Eqs. (15.42), 
(15.1-15.5), and (15.146-15.147)) reproduces the results for benzene given in the Benzene 
Molecule (C 6 H 6 ) section as shown in Tables 15.214 and 15.216. 

The symbols of the functional groups of aromatics and hertocyclics are given in Table 
15.213. The geometrical (Eqs. (15.1-15.5) and (15.42)), intercept (Eqs. (15.7l-15.78)), and 
energy (Eqs. (15.6-15,11), (15.17-15.56), and (15.146-15.147)) parameters of aromatics and 
hertocyclics are given in Tables 15.214, 15.215, and 15.216, respectively. The total energy of 
benzene given in Table 15.217 was calculated as the sum over the integer multiple of each 
E D (group) of Table 15.216 corresponding to functional-group composition of the molecule. The 

bond angle parameters of benzene determined using Eqs. (15.79-15.108) are given in Table 
15.218. 



Table 15.213. The symbols of functional groups of aromatics and hertocyclics. 



Functional Group 


Group Symbol 


CC (aromatic bond) 


3e 

C=C 


CH (aromatic) 


CH(i) 
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Table 15.214.. The geometrical bond parameters of aromatics and hertocyclics and 



experimental values [1], 



Parameter 


3e 

— V^- 

Group 


CH 
Group 


a ( a o) 


1.47348 


1.60061 




1.31468 


1.03299 


Bond Length 2c' (□) 


1.39140 


1.09327 


Exp. Bond Length 

(□) 


1.399 
(benzene) 


1.101 

(benzene) 


b,c (a 0 ) 


0.66540 


1.22265 


e 


0.89223 


0.64537 



Table 15.216. The energy parameters (eV) of functional groups of aromatics and 
heterocyclics. 



Parameters 


3e 

C=C 
Group 


CH 
Group 


A. 


0.75 


1 


n l 


2 


1 




0 


0 




0 


0 




0.5 


0.75 


c 2 


0.85252 


1 


c l 


1 


1 


C 2 


0.85252 


0.91771 


C l 


0 


1 


C 4 


3 


1 


C 5 


0 


1 




0.5 


0.75 




0.85252 


1 


K (eV) 


-101.12679 


-37.10024 


V p (eV) 


20.69825 


13.17125 


T(eV) 


34.31559 


11.58941 




-17.15779 


-5.79470 


E(aoiho) (eV) 


0 


-14.63489 




0 


-1.13379 



58 



WO 2007/053486 



PCT/US2006/042140 



TP ( Am rrr\ ( ^T/\ 

Et T [ao/ ho) \ev ) 


u 


-I J.5U1 1U 








E T [atom — atom, msp .AO \ (eV) 


-2.26759 


-0.56690 


E T ( ¥ o){eV) 


-65.53833 


-32.20226 


a> (lO 15 radl 


49.7272 


26 4826 


E (eV) 


32.73133 


17 43132 


E (eV) 


-0 35806 


-0 261 30 


77* / T7"\ 

E K Vi t ( eV ) 


0.19649 
T49] 


0.35532 
Eq. (13.458) 


E (eV) 

osc ^ / 


-0.25982 


-0.08364 


E (eV) 

mag v s 


0.14803 


0.14803 


E T (group) (eV) 


-49.54347 


-32.28590 




-14.63489 


-14.63489 




0 


-13.59844 


E D (group) (eV) 


5.63881 


3.90454 
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■J 3 


0.79597 


0.79597 


J3 


r- 
o\ 
o 


T-H 

OS 

o 


Final Total 
Energy 

C2sp 3 
(eV) 1 


CO 

oo 
oo 

CO 

in 
*— < 
i 


cn 

CO 
CO 

oo 

CO 

in 
■ 


CP 


o 


o 


PQ 


o 

SO 

so 
tn 


o 

OS 
SO 

so 
m 

9 


CN 

PQ 


m 
ro 
O 
in 
OO 

9 


*n 

o 
in 
oo 

9- 


' PQ 


in 
ro 
O 
in 
oo 

9 


m 

,s 
</■> 

oo 

9 


Atom 


O 




Bond 


l 

U 


O 

<3 II 

<S II 

O 



^"3 


21379 


.55533 




o 


o 


^"3 


oo 


0.75935 




GO 
CO 

cn 


OO 
OS 
OO 

m 




oo 
in 

in 
© 


so 
f": 
in 


CD O 




Si 

CO 
1— 1 




s 

as 

SD 


s 

OS 

SO ' 




i 


> — 1 
1 










■17.09334 


•17.09334 


1 












Bond 


£ 

ft! 
1 

O 


£ It 
3 II 



CO 

* i— i 

§ 



Pi 

:1 

CO 

O 

& 

O 

o 

§, 

•I 







Relative Errc 


SO 

o 
o 


© 
o 


S 




i § 

J « 


.26340 




s* ° 
« ffl 


in 


I 








Calculated 
tal Bond Energy 
(eV) 


57.26008 


Toi 




CH 


so 


O 

<S 11 

u 










S 

3 








1 
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P « a 

^ 2 

< 


0.79232 




< 


0.79232 




Atom 2 
Hybridization 
Designation 

(Table 
15.3.A) 


m 




•i 1 


-17.17218 




Atoml 
Hybridization 
Designation 

j (Table 
1 15.3.A) 






E 

Couiombtc 
Atom 1 


-17.17218 




- "is.— 


4.5585 




c« , s 

CQ s — ' 


2.62936 




CQ s — ' 


2.62936 




Atoms of 
Angle 


ZCCC 
(aromatic) 


ZCCH 
(aromatic) 



<3s 


120 [50-52] 
(benzene) 


120 [50-52] 
(benzene) 




120.19 


119.91 












120.19 










-1.85836 




V 


0.79232 










Atoms of Angle 


ZCCC 
(aromatic) 


ZCCH 
(aromatic) 



61 



WO 2007/053486 



PCT/US2006/042140 



REFERENCES FOR DETAILED DESCRIPTION OF INVENTION 

1 . D. R. Lide, CRC Handbook of Chemistry and Physics, 86th Edition, CRC Press, Taylor & 
Francis, Boca Raton, (2005-6), pp. 9-19 to 9-45. 

2. G. Herzberg, Molecular Spectra and Molecular Structure II Infrared and Raman Spectra 
of Polyatomic Molecules, Van Nostrand Reinhold Company, New York, New York, 
(1945), p. 344. 

3. D. R. Lide, CRC Handbook of Chemistry and Physics, 86th Edition, CRC Press, Taylor & 
Francis, Boca Raton, (2005-6), pp. 9-63; 9-64 to 9-76. 

4. R. J. Fessenden, J. S. Fessenden, Organic Chemistry, Willard Grant Press. Boston, 
Massachusetts, (1979), p. 320. 

5 . cyclohexane at http://webbook.nist.gov/ . 

6. G. Herzberg, Molecular Spectra and Molecular Structure II Infrared and Raman Spectra 
of Polyatomic Molecules, Van Nostrand Reinhold Company, New York, New York, 
(1945), p. 326. 

7 . b 1 2 -methyl- 1 -propene at http ://webbook.nist. gov/ . 

8 . al 2-methyl- 1 -propene at http ://webbook.nist . gov/ . 

9. 2-butyne at http ://webbook.nist gov/ . 

1 0. trifluoromethane at http://webbook.nist.gov/ . 

1 1 . fluoroethane at http://webbook.nist.gov/ . 

12. G. Herzberg, Molecular Spectra and Molecular Structure II Infrared and Raman Spectra 
of Polyatomic Molecules, Krieger Publishing Company, Malabar, FL 5 (1991), p. 195. 

13. chloromethane at http://webbook.nist.gov/ . 

14. methyl bromide at http://webbook.nist.gov/ . 

15. tetrabromomethane at http ://webbook.nist gov/ . 

16. methyl iodide at http://webbook.nist.gov/ . 

17. K. P. Huber, G. Herzberg, Molecular Spectra and Molecular Structure, IK Constants of 
Diatomic Molecules, Van Nostrand Reinhold Company, New York, (1979). 

18. J. Crovisier, Molecular Database — Constants for molecules of astrophysical interest in the 
gas phase: photodissociation, microwave and infrared spectra, Ver. 4.2, Observatoire de 
Paris, Section de Meudon, Meudon, France, May 2002, pp. 34-37, available at 
http://wwwusr.obspm.fr/~crovisie/. 

19. methanol at http://webbook.nist. gov/ . 



62 



WO 2007/053486 



PCT/US2006/042140 



20. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 46. 

21 ( dimethyl ether at http ://webbook,nist. gov/ . 

>2. T. Amano, P. F. Bernath, R. W. McKellar, "Direct observation of the v Y and v 3 
fundamental bands of NH 2 by difference frequency laser spectroscopy", J. MoL 
Spectrosc, Vol. 94, (1982), pp. 100-113. 

13. methylamine at http ://webbook.nist gov/ . 

24. D. R. Lide, CRC Handbook of Chemistry and Physics, 79th Edition, CRC Press, Boca 
Raton, Florida, (1998-9), pp. 9-80 to 9-85. 

25. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 482. 

26. acetaldehyde at http ://webbook,nist gov/ , 

27. acetone at http://webbook.nist.gov/ . 

28. 2-butanone at http ://webbook.nist. gov/ . 

29. acetic acid at http://webbook.nist. gov/ . 

30. formic acid at http://webbook.nist. gov/ . 

31. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 138. 

32. methyl formate at http://webbook.nist.gov/ . 

33. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, Th& Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 144. 

34. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 147. 

35. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 143. 



63 



WO 2007/053486 



PCT/US2006/042140 



36. H. J. Vledder, F. C. Mijlhoff, J. C. Leyte, C. Romers, "An electron diffraction 
investigation of the molecular structure of gaseous acetic anhydride", Journal of 
Molecular Structure, Vol. 7, Issues 3-4, (1971), pp. 421-429. 

37. acetonitrile at http://webbook.nist. gov/ . 

$8. D. R. Lide, CRC Handbook of Chemistry and Physics, 86th Edition, CRC Press, Taylor & 
Francis, Boca Raton, (2005-6), pp. 10-202 to 10-204. 

39. D. R. Lide, CRC Handbook ofChemisPy and Physics, 86th Edition, CRC Press, Taylor & 
Francis, Boca Raton, (2005-6), p. 9-79. 

40. D. R. Lide, CRC Handbook of Chemistry and Physics, 86th Edition, CRC Press, Taylor & 
Francis, Boca Raton, (2005-6), pp. 9-82 to 9-86. 

41 . thiirane at htt p://webbook.nist.gov/ . 

42. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 242. 

43. D. Lin-Vien. N. B, Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 241. 

44. R. M. Ibberson, M. T. F. Telling, S. Parsons, "Structural determination and phase 
transition behavior of dimethyl sulfate", Acta. Cryst, Vol. B62, (2006), pp. 280-286. 

45. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 181. 

46. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 480. 

47. D. Lin-Vien. N B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 187. 

48. 1.3-butadiene at http://webbook.nist.gov/ . 

49. G. Herzberg, Molecular Spectra and Molecular Structure II Infrared and Raman Spectra 
of Polyatomic Molecules, Van Nostrand Reinhold Company, New York, New York, 
(1945), pp. 362-369. 



64 



WO 2007/053486 



PCT/US2006/042140 



50, W. L F. David, R. M. Ibberson, G. A. Jeffrey, J. R. Ruble, "The structure analysis of 
deuterated benzene and deuterated nitromethane by pulsed-neutron powder diffraction: a 
comparison with single crystal neutron analysis", Physica B (1992), 180 & 181, pp. 597- 
600.. 

>1.G. A. Jeffrey, J. R. Ruble, R. K. McMullan, J. A. Pople, "The crystal structure of 
deuterated benzene," Proceedings of the Royal Society of London. Series A, 
Mathematical and Physical Sciences, Vol. 414, No. 1846, (Nov. 9, 1987), pp. 47-57. 

52. H. B. Burgi, S. C. Capelli, "Getting more out of crystal-structure analyses," Helvetica 
Chimica Acta, Vol. 86, (2003), pp. 1625-1640. 

53. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 481. 

54. R. Boese, D. B laser, M. Nussbaumer, T. M. Krygowski, "Low temperature crystal and 
molecular structure of nitrobenzene", Structural Chemistry, Vol. 3, No. 5 5 (1992), pp. 
363-368. 

55. G. A. Sim, J. M. Robertson, T. H. Goodwin, "The crystal and molecular structure of 
benzoic acid", Acta Cryst., Vol 8, (1955), pp. 1 57- 1 64 

56. D. Lin-Vien. N. B. Colthup, W, G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 301. 

57. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
jovanovich, Boston, (1991), p. 303. 

58. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), p. 207. 

59. S. Martinez-Carrera, "The crystal structure of Imidazole at -150°C", Acta. Cryst., Vol. 20, 
(1966), pp. 783-789. 

60. D. Lin-Vien. N. B. Colthup, W. G. Fateley, J. G. Grasselli, The Handbook of Infrared and 
Raman Frequencies of Organic Molecules, Academic Press, Inc., Harcourt Brace 
Jovanovich, Boston, (1991), pp. 296-297. 



65 



WO 2007/053486 



PCT/US2006/042140 



61. H. Sternglanz, C. E. Bugg, "Conformations of N 6 -monosubstituted adenine derivatives 
crystal structure of N 6 ~methyladenine" ? Biochimica et Biophysica Acta, Vol. 308, (1973), 
pp. 1-8. 



66 



WO 2007/053486 



PCT/US2006/042140 



Software Program 

The present invention relates to a system of computing and rendering the nature of at 
least one specie selected from a group of diatomic molecules having at least one atom that is 
other than hydrogen, polyatomic molecules, molecular ions, polyatomic molecular ions, or 
molecular radicals, or any functional group therein, comprising physical, Maxwellian solutions 
of charge, mass, and current density functions of said specie, said system comprising: 
processing means for processing physical, Maxwellian equations representing charge, mass, 
and current density functions of said specie; and an output device in communication with the 
processing means for displaying said physical, Maxwellian solutions of charge, mass, and 
current density functions of said specie. 

In one embodiment, for example, the system comprises five components: (1) the 
graphical user interface (GUI); (2) the routine for parsing between an input chemical structure 
or name and taking the input and activating a routine to call up the parts of the molecule 
(functional groups), which are used for determining the energies and structure to be rendered; 
(3) the functional-group data base that has an organization of the theoretical solutions; (4) the 
rendering engine, which calculates and enables manipulations of the image, such as a three- 
dimensional model in response to commands, as well as responds to commands for data 
parameters corresponding to the image such as bond energies and charge distribution and 
geometrical parameters; and (5) data transfer system for inputting numerical data into or out of 
the computational components and storage components of the main system. The system 
further comprises spreadsheets with solutions of the bond parameters with output in any 
standard spreadsheet format. The system also comprises a data-handling program to transfer 
data from the spreadsheets into the main program. 

The output may be, for example, at least one of graphical, simulation, text, and 
numerical data. The output may be the calculation of at least one of: (1) a bond distance 
between two atoms; (2) a bond angle between three of the atoms; (3) a bond energy between 
two atoms; (4) orbital intercept distances and angles; and (5) charge-density functions of 
atomic, hybridized, and molecular orbitals, wherein the bond distance, bond angle, and bond 
energy are calculated from physical solutions of the charge, mass, and current density 
functions of atoms and atomic ions, which solutions are derived from Maxwell's equations 
using a constraint that a bound electron(s) does not radiate under acceleration. 

In other embodiments, the charge, current, energy, and geometrical parameters are 
output to be inputs to other programs that may be used in further applications. For example, 
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the data of heats of formation may be input to another program to be used to predict stability 
(existence of compounds) equilibrium constants and to predict synthetic pathways. That is, a 
novel composition of matter may be discovered by calculating at least one of a bond distance 
between two of the atoms, a bond angle between three of the atoms, and a bond energy 
between two of the atoms, orbital intercept distances and angles, charge-density functions of 
atomic, hybridized, and molecular orbitals, the bond distance, bond angle, and bond energy 
being calculated from physical solutions of the charge, mass, and current density functions of 
atoms and atomic ions, which solutions are derived from Maxwell's equations using a 
constraint that a bound electron(s) does not radiate under acceleration. The charge and current 
density functions may be used to predict the electric and magnetic fields of the species to 
determine other properties due to the interaction of the fields between species. These fields 
and the predictions of field interactions may be computed using Maxwell's equations. In one 
embodiment, finite-element analysis is used to predict or calculate the interaction and resulting 
properties, such as the freezing point, boiling point, density, viscosity, and refractive index. 
Furthermore, the output data can be used to give thermodynamic, spectroscopic, and other 
properties, aid in drug design and other applications with or without direct visualization. 
Furthermore, the data can be input into other programs of the system, which calculate 
thermodynamic and other properties, or performs a simulation, such as a chemical reaction or 
molecular dynamics. 

The output data may be used to predict a composition of matter comprising a plurality 
of atoms, the improvement comprising a novel property or use discovered by calculation of at 
least one of a bond distance between two of the atoms, a bond angle between three of the 
atoms, and a bond energy between two of the atoms, orbital intercept distances and angles, 
charge-density functions of atomic, hybridized, and molecular orbitals, the bond distance, bond 
angle, and bond energy being calculated from physical solutions of the charge, mass, and 
current density functions of atoms and atomic ions, which solutions are derived from 
Maxwell's equations using a constraint that a bound electron(s) does not radiate under 
acceleration. The novel property, for example, may be a new pharmaceutical use, or stability 
at room temperature of a novel arrangement of atoms or ions. 

In one embodiment, the output device of the system is a display that displays at least 
one of visual or graphical media. The display may be at least one of static or dynamic. At 
least one of vibration, rotation, and translation may be displayed. The displayed information 
may be used for at least one of modeling reactivity, predicting physical properties, and aiding 
in drug and material design. The output device may be a monitor, video projector, printer, or 
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one-, two- or three-dimensional rendering device. The displayed information may be used to 
model other molecules and provides utility to anticipate their reactivity and physical properties. 
Additionally, data may be output and used in the same and additional applications as the 
rendered models and representations of the calculated physical solutions. The processing 
means of the system may be a general-purpose computer. The general-purpose computer may 
comprise a central processing unit (CPU), one or more specialized processors, system memory, 
a mass storage device such as a magnetic disk, an optical disk, or other storage device, an input 
means. The input means may comprise a serial port, USB port, microphone input, camera 
input, keyboard or mouse. The processing means comprises a special purpose computer or 
other hardware system. The system may comprise computer program products such as 
computer readable medium having embodied therein program code means. The computer 
readable media may be any available media which can be accessed by a general purpose or 
special purpose computer. The computer readable media may comprise, for example, at least 
one of RAM, ROM, EPROM, CD ROM, DVD or other optical disk storage, magnetic disk 
storage or other magnetic storage devices, or any other medium that can embody the desired 
program code means and which can be accessed by a general purpose or special purpose 
computer. The program code means may comprise executable instructions and data, which 
cause a general purpose computer or special purpose computer to perform a certain function of 
a group of functions. Commercial examples of suitable program language includes, for 
example, C++, C, JAVA, FORTRAN, Python and Assembly Languages, programmed with an 
algorithm based on the physical solutions, and the computer may be a PC, mainframe, 
supercomputer, or cluster of computers. Commercial examples of suitable programs include, 
for example, APIs like OpenGL, DirectX, FOX GUI toolkit, and Qt. This program may be 
developed to run on at least one of operating systems like Windows XP, Windows 2000, 
Windows Vista, MAC OS, MAC OS X, Linux, Unix, Irix and other Unix-type operating 
systems. 

Millsian software is designed to render 3-D models of molecules, molecular ions, 
molecular radicals, functional groups thereof, and related structure and property information 
and produce useful data output and application of the parameters of these species, wherein the 
nature of their bound electrons and chemical bonds are solved using Dr. Randell L. Mills' 
Classical Quantum Mechanics theory described in the reference: R. L. Mills, "The Grand 
Unified Theory of Classical Quantum Mechanics", June 2006 Edition, Cadmus Professional 
Communications-Science Press Division, Ephrata, PA, ISBN 0963517171, Library of 
Congress Control Number 2005936834; posted at 
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http://wvw.blacldightpow which is incorporated by this reference 

in its entirety. Figure 76 provides a flow chart diagram, which is an example of a software 
system that can be utilized for this purpose, which example is not intended to limit the scope of 
the disclosed inventions. The main parts of this exemplary software system illustrated in 
Figure 1 will now be further explained: 

Start: A user can start the program by running an executable program file. That might 
be done, for example, by double clicking the program icon on a Windows-based operating 
system, or typing the name of the executable file on the command line and pressing the 'Enter 5 
key on a Linux or Unix operating system. The program initially starts by reading data files 
located in specific directories. The names of those directories and locations are fixed 
according to the type and format of the data files. 

Data files: There are two types of data files used in the software system: functional 
group data files and molecule data files. Functional group data files contain information about 
various functional groups. Functional groups are the basic bonding elements or units that each 
typically comprise an atom, or at least two atoms bound together as found within a molecule 
(e.g. -CI , C=C, C=(X CH3). Functional groups typically dictate or define properties and 
structure of the molecule. Similar functional groups in different molecules typically react in 
similar ways when subjected to a particular set of reaction conditions. Molecule data files 
contain information about molecules, molecular ions, and molecular radicals. These data files 
are processed according to their file formats. 

File formats: Millsian software employs two kinds of file formats for storing 
information about the structure, energies and names of molecules and functional groups: raw- 
data format and hierarchical format. In a raw-data file, all information is stored as is, below 
the header describing the type of information. As shown in the Table below, for example, the 
names of the molecule or functional group are listed below the #NAMES header. The names 
and positions of the atoms are listed below the #ATOMS header, and so on. 
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#NAMES 




CH3 




Alkane CH3 




#SMILES 




C- 




#ATOMS 




1 C 0 0 


0 


2 H 1.9775 -0.698 


0 


3 H -0.9888 


-0.698 1.712 


4 H -0.9888 


-0.698 -1.712 


#AO 




1 1 Is 0.171 2 




2 1 2sp3 0.864 -1 




#BONDS 




1 12 1 1.649 




2 13 1 1.649 




3 14 1 1.649 




#BONDAXES 




11 1 90 


0 0 


#DATA 




RCHbond angle: 109.44° 




CH bond length: 2.097 a 




HFORM: 12.492 





Most of the functional groups files and some simple molecule files are stored in raw 
data format. The other file format used with the Millsian software system, the hierarchical file 
format, represents the information in a graph style, in which nodes are connected to other 

i 

nodes through links. This format is designed to construct molecules by attaching different 
functional groups in a desired manner. Using this file format, the user can construct complex 
molecules built from solved functional groups. 

For example, as shown in the Table below, the pentane data file includes a #GROUPS 
header below which is list all of the functional groups that form a part of the pentane molecule. 
Under the #GROUP_LINKS header is information about how these functional groups are 
connected to each other to construct pentane. 



71 



WO 2007/053486 



PCT/US2006/042140 



#NAMES 
Pentane 

#SMILES 

ccccc 

#GROUPS 

1 C- 

2 -C- 

3 -C- 

4 -C- 

5 C- 



#GROUP LINKS 



1 


1 1 


2 1 


180 


2 


22 


3 1 


180 


3 


32 


4 1 


180 


4 


42 


5 1 


180 



Processing data files: As further shown in the flow diagram of Figure 1 for the 
exemplary Millsian software system, the program first processes the functional group data file 
and constructs the functional-group objects, which are complete 3-D representations of the 
functional groups and their related information supplied in the corresponding file. These 
objects are then ready to be visually displayed through use of a molecule viewer. Next, the 
program processes the molecule data file. If the molecule date file is in raw-data format, then 
the program makes a molecule object directly from it. If the file is in hierarchical format, then 
the program calculates geometric parameters from listed functional groups. The program 
stores all functional-group objects and molecule objects using internal data structures. 

Visualization/ User Interactions: As shown in Figures 2 and 3, the molecule viewer 
displays the functional- group objects and molecule objects and provides basic interaction 
capabilities with the displayed objects, such as rotating, scaling, and moving the objects. The 
molecule viewer also provides other visualization options, such as viewing molecules in wire 
frame mode, viewing coordinate axes, and changing of the transparency and lighting. The 
user, for example, can also select parts of a molecule for visualization, like a nucleus, atomic 
orbital, molecular orbital, or bond axis. The viewer also includes a drop down information 
window, which provides, for example, related information about molecules, such as bond 
angles, component functional groups, and total heat of formation. In one embodiment, a user 
can create new molecules by joining functional groups. For this purpose, the user can select a 
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first functional group. Next, the user can select an open bond from the functional group where 
the user desires to attach another group. Next, the user can select another or the same 
functional group, followed by selecting an open bond from the other group. The user can join 
the two selected functional groups at the selected open bonds by clicking on 'Join Groups'. 
This method of joining functional groups at open bonds can be repeated to form the desired 
molecule. 

SMILES input: The Simplified Molecular Input Line Entry Specification or SMILES is 
a specification for unambiguously describing the structure of chemical molecules using short 
ASCII strings. Through the user interface, a user can enter SMILES to construct molecules, 
provided the new molecule comprises functional groups that are in the database. Once a 
SMILES is entered, a parser reads in and breaks it down into component functional groups. 
The software system then attaches the component functional groups to create the new molecule 
object, which can then be viewed using the molecule viewer. 

While the claimed invention has been described in detail and with reference to specific 
embodiments thereof, it will be apparent to one of ordinary skill in the art that various changes 
and modifications can be made to the claimed invention without departing from the spirit and 
scope thereof. 
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CLAIM: 

1. A system of computing and rendering the nature of at least one specie selected from a 
group of diatomic molecules having at least one atom that is other than hydrogen, 
polyatomic molecules, molecular ions, polyatomic molecular ions, or molecular radicals, 
or any functional group therein, comprising physical, Maxwellian solutions of charge, 
mass, and current density functions of said specie, said system comprising: 

processing means for processing physical, Maxwellian equations representing charge, 
mass, and current density functions of said specie; and 

an output device in communication with the processing means for displaying said 
physical, Maxwellian solutions of charge, mass, and current density functions of said specie. 

2. The system of claim 1 wherein the output device is a display that displays at least one of 
visual or graphical media. 

3. The system of claim 2 wherein the display is at least one of static or dynamic. 

4. The system of claim 3 wherein at least one of vibration and rotation is be displayed. 

5. The system of claim 4 wherein displayed information is used to model reactivity and 
physical properties. 

6. The system of claim 5, wherein the output device is a monitor, video projector, printer, or 
three-dimensional rendering device. 

7. The system of claim 6 wherein displayed information is used to model other species and 
provides utility to anticipate their reactivity and physical properties. 

8. The system of claim 7 wherein the processing means is a general purpose computer. 

9. The system of claim 8 wherein the general purpose computer comprises a central 
processing unit (CPU), one or more specialized processors, system memory, a mass 
storage device such as a magnetic disk, an optical disk, or other storage device, an input 
means. 
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10. The system of claim 9, wherein the input means comprises a serial port, usb port, 
microphone input, camera input, keyboard or mouse. 

1 1 . The system of claim 10 wherein the processing means comprises a special purpose 
computer or other hardware system. 

12. The system of claim 1 1 further comprising computer program products. 

13. The system of claim 12 further comprising computer readable medium having embodied 
therein program code means. 

14. The system of claim 13 wherein the computer readable media is any available media 
which can be accessed by a general purpose or special purpose computer. 

15. The system of claim 14 wherein the computer readable media comprises at least one of 
RAM, ROM, EPROM 5 CD ROM 5 DVD or other optical disk storage, magnetic disk 
storage or other magnetic storage devices, or any other medium which can embody the 
desired program code means and which can be accessed by a general purpose or special 
purpose computer. 

16. The system of claim 15 wherein the program code means comprises executable 
instructions and data which cause a general purpose computer or special purpose 
computer to perform a certain function of a group of functions. 

17. The system of claim 16 wherein the program code is Mathematica programmed with an 
algorithm based on the physical solutions, and the computer is a PC. 

18. The system of claim 17 wherein the algorithm is ParametricPlot3D[{2 >H Sqrt[l- 
z*z]*Cos[u],Sqrt[(l-z*z)]*Sin[u],z} 9 {u,0,2* Pi},{z,-1,.9999}], and the rendering is 
viewed from different perspectives. 
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19. The system of claim 18 wherein the algorithms for viewing from different perspectives 
comprises Show[Out[l], ViewPoint->{x,y 5 z}] where x, y, and z are Cartesian 
coordinates, 

20. The system of claim 19 wherein the physical, Maxwellian solutions of the charge, mass, 
and current density functions of said specie comprises a solution of the classical wave 



equation i V 2 -—— w(r 9 0,^,0 = 0 which is the equation of motion of the charge. 



21 . The system of claim 20 wherein the boundary constraint of the wave equation solution is 
nonradiation according to Maxwell's equations. 

22. The system according to claim 1, wherein the force generalized constant k 1 of a H 2 -type 



where Cj is the fraction of the H 2 -type ellipsoidal molecular orbital basis function of a 
chemical bond of the specie and C 2 is the factor that results in an equipotential energy match 
of the participating at least two molecular or atomic orbitals of the chemical bond. 

23. The system according to claim 1 ? wherein the distance from the origin of the MO to each 
focus c ' is given by: 




ellipsoidal MO due to the equivalent of two point charges at the foci is given by: 



CC 2e 2 



(15.1) 



4tts 0 




(15.2) 



the internuclear distance is 




(15.3) 



the length of the semiminor axis of the prolate spheroidal MO b = c is given by 



b = «Ja 2 -c' 



(15.4) 
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and, the eccentricity, e , is 



£- (15.5). 

a 



24. The system according to claim 23, wherein a potential energy of the electrons in the 
central field of the nuclei at the foci is 



T _ -2e 2 . a + 'sja 2 -b 2 n<&\ 

Sxe^cP-b 2 a-\la 2 -b 2 

where 7^ is the number of equivalent bonds of the MO for functional groups and in the case 

of independent MOs not in contact with the bonding atoms, the terms based on charge are 

multiplied by c BO , the bond-order factor, which is 1 for a single bond, 4 for an independent 

double bond and 9 for an independent triplet bond. 

25. The system according to claim 1, wherein the potential energy of the two nuclei is 

V p =fh < 15 - 7 ^ 

%7teja 2 -b 2 

26. The system according to claim 1 3 wherein the kinetic energy of the electrons is 

T = I 2 ,2^ 1 2 ,2 (15 * 8) 

2m e ayja —b a-^a -b 
where n x is the number of equivalent bonds of the MO for functional groups and in the case 
of independent MOs not in contact with the bonding atoms, the terms based on charge are 
multiplied by c BQ , the bond-order factor, which is 1 for a single bond, 4 for an independent 
double bond and 9 for an independent triplet bond. 

27. The system according to claim 1 5 wherein the energy, V m 5 of the magnetic force between 
the electrons is 

V m = n x c x c 2 , In . (15.9) 
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where is the number of equivalent bonds of the MO for functional groups and in the case 
of independent MOs not in contact with the bonding atoms, the terms based on charge are 
multiplied by c BO , the bond-order factor, which is 1 for a single bond, 4 for an independent 
double bond and 9 for an independent triplet bond. 

28. The system according to claim 1, wherein total energy of the H 2 -type prolate spheroidal 
MO, E T (h 2 mo) , is given by the sum of the energy terms: 

E T M = V e +T + V m +V p (15.10) 



87tsJa 2 -b 2 



c x c 2 \ 2 — 2- 



a 0 \ a + ^)a 2 -b 2 1 
In , - 1 



o 

a j 



a-4a 2 -h 2 



2 r 



87V£ Q C ' |_ 



c x c 2 \ 2 



a J a-c' 



(15.11) 



where is the number of equivalent bonds of the MO for functional groups and in the case 
of independent MOs not in contact with the bonding atoms, the terms based on charge are 
multiplied by c BO 9 the bond-order factor, which is 1 for a single bond, 4 for an independent 
double bond and 9 for an independent triplet bond. 



29- The system according to claim 1, wherein the total energy E T [atom, msp 3 j (m is the 

integer of the valence shell) of the AO electrons and the hybridized shell is given by the 
sum of energies of successive ions of the atom over the n electrons comprising total 
electrons of the at least one AO shell according to the formula: 

E T (atom,msp 3 ) = ~^IP m (15.12 
where IP m is the m th ionization energy (positive) of the atom. 



30. The system according to claim 29, wherein the radius r^ 3 of the hybridized shell is given 
by: 

r >B § " ( f" g>2 3 , (15.13). 

msp q =z-n 8rcs^E T iatom, msp J 
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3 1 . The system according to claim 29, wherein the Coulombic energy E Coulomh [atom, msp 3 ) 
of the outer electron of the atom msp 3, shell is given by 



-e 2 



Ecou^b (atom,msp 3 ) = (15-14). 



u msp' 



32. The system of claim 31 ? wherein in the case that during hybridization at least one of the 
spin-paired AO electrons is unpaired in the hybridized orbital (HO), the energy change for 
the promotion to the unpaired state is the magnetic energy E (magnetic) at the initial 
radius r of the AO electron: 



E(magnetic) = 



= 2xju 0 e 2 h 2 = *w 0 nl (15.15) 



m 2 r 3 r 3 



then, the energy E [atom, msp 3 ) of the outer electron of the atom msp 3 shell is given by the 
sum of E Coulomb [atom, msp 3 ) and E (magnetic) : 

E (atom, msp 3 ) = — — 1- (15.16). 

33. The system according to claim 29, wherein the total energy E T [mol.atom,msp 3 ) (m is 

the integer of the valence shell) of the HO electrons is given by the sum of energies of 
successive ions of the atom over the n electrons comprising total electrons of the at least 
one initial AO shell and the hybridization energy: 

» 

E T [moLatom.msp 3 ) = E [atom, msp 3 ) - ^ IP m (15.17) 

where IP m is the m th ionization energy (positive) of the atom and the sum of -IP X plus the 
hybridization energy is E[atom 9 msp 3 ) . 

34. The system of claim 33, wherein the radius r 3 of the hybridized shell is given by: 
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r , = 



\q=Z~n 



%7te Q E T {mohatom, msp 3 ) 
where s = 1,2,3 for a single, double, and triple bond, respectively. 



(15.18) 



35. The system of claim 34, wherein the Coulombic energy E Coulomh {moLatom,msp 3 ) of the 
outer electron of the atom msp 3 shell is given by 



— e 2 

E couiom b (mol.atom,msp 3 ) = ~ (15.19). 



u msp* 



36. The system of claim 35, wherein in the case that during hybridization at least one of the 
spin-paired AO electrons is unpaired in the hybridized orbital (HO), the energy change for 
the promotion to the unpaired state is the magnetic energy E(magnetic) at the initial 

radius r of the AO electron given by Eq. (15.15). Then, the energy E{moLatom,msp 3 ) of 

the outer electron of the atom msp 3 shell is given by the sum of E Coulomb [mol. atom, msp 3 ) 
and EQnagnetic) : 

E { moLatom, msp 3 ) = + J 1 (15.20) 

v , ' %7vs Q r 3 m;r 

, u msp 3 e 

E T [atom - atom, msp 3 ) , the energy change of each atom msp 3 shell with the formation of 
the atom-atom-bond MO is given by the difference between E^mol.atom, msp 3 ) and 
E[atom,msp 3 ) : 

E T {atom - atom, msp 3 ) = E(mol.atom, msp 3 )~E [atom, msp 3 ) (1 5.21). 

37. The system of claim 31, wherein E Coulomb (mol atom, msp 3 ) is one of: 

E Coulomb ( C ethylene > ^Sp 3 ) s E C oulomb (C ethanQ , 2sp 3 ) 9 E Cou}omb [C acetyhm , 2sp 3 ) , and 
^Coulomb {palkam ' 2>Sp 3 ) ; 

E couio,nb(atom,msp 3 ) is one of E Couhmb (c,2sp 3 ) and E Coulomb (Cl,3sp 3 ); 

E{moLatom,msp 3 ) is one of E(C ethylene ,2sp 3 ), E(C etham ,2sp 3 ) , 



80 



WO 2007/053486 



PCT/US2006/042140 



E(atom,msp 3 ) is one of and E(C, 2 sp 3 ) and E(Cl,3sp 3 ); 
E T (atom - atom, msp 3 ) is one of E (C - C, 2sp 3 ), E(C = C, 2sp 3 ) , and. 

E[C = C,2sp 3 ); 

atom msp 3 is one of C2sp 3 , CISsp 3 
E T (atom -atom(s 1 ), msp 3 ) is E T (c - C, 2sp 3 ) and E T (atom- atom (s 2 ), msp 3 ) is 

E T (C = C,2sp 3 ), and 

In the case of the C2sp 3 HO, the initial parameters (Eqs. (14.142-14.146)) are 

r ,=y (Z ~ w)g2 - = lOf! - = 0.91771a, (15.22) 

2sp %£toB 0 (el4&J25751eV) 8xe 0 (el48.25751 eV) 

1 2 

E r ' Ac,2sp 3 ) = ^- = — = -14.82575 eV (15.23) 

co»,o mb \ » P ) sxe^ 8^ 0 0.91771o 0 

Eimagnetic) = ^wf§l = = 0.19086 ,F (15.24) 

(0.84317a 0 ) 3 



£(C,2jrp 3 ) = 



-e 2 2n/j 0 e 2 h 2 



- + ■ 



8 ^oV ™ 2 (r 3 ) 

= -14.82575 eF + 0.19086 (15.25) 
= -14.63489 eV 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 

38. The system of claim 37, wherein 

]T(Z-?) = 10 (15.26) 

q=Z-n 

Equations (1 4. 1 47) and (1 5 . 1 7) give 

E T (moLatom,msp 3 )^ E T (C etham ,2sp 3 ) = -151.61569 eV (15.27) 
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and using Eqs. (15.18-15.28), the final values of r C2sp , , E Coulomb (C2sp 3 ), and E(C2sp 3 ),and 
the resulting E T \ C ^ C, C2sp 3 ] of the MO due to charge donation from the HO to the MO 



BO 



where C - C refers to the bond order of the carbon-carbon bond for different values of the 
parameter s are given in Table 15.1. 



Table 15.1: 



MO 
Bond 
Order 
(BO) 


s 1 


s 2 


Final 


-Ec— (C2*y)(eV) 
Final 


E(C2sp 3 ) 

(eV) 
Final 


E T [c-C,C2sp 3 ^j 
(eV) 


I 


1 


0 


0.87495 


-15.55033 


-15.35946 


-0.72457 


n 


2 


0 


0.85252 


-15.95955 


-15.76868 


-1.13379 


in 


3 


0 


0.83008 


-16.39089 


-16.20002 


-1.56513 


IV 


4 


0 


0.80765 


-16.84619 


-16.65532 


-2.02043 



wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 

39. The system of claim 33, wherein a minimum-energy bond with the constraint that it must 
meet the energy matching condition for all MOs at all HOs or AOs, the energy 
E{rnoLatom, msp 3 ) of the outer electron of the atom msp 3 shell of each bonding atom is 

the average of E(inolatom, msp 2 ) for two different values of s : 

, „ x EimoLatom (s } ) , msp* ) + Eimol.atom (s 7 ) 9 msp 3 ) 
E(mol.atom,msp 2 ) = -± K 1,9 F} K V 2) 1 (15.28). 



40. The system of claim 39, wherein in the case, E T {atom - atom, msp 3 ) , the energy change 

of each atom msp 3 shell with the formation of each atom-atom-bond MO, is average for 
two different values of s : 



82 



WO 2007/053486 



PCT/US2006/042140 



„ / -x E T ( atom - atom ( s*) , msp 3 ) + E T (atom- atom (s 0 ), msp 2 ) 

E T [atom - atom, msp 3 ) = — ^ — — ; ^ J ^ v 2J — t—l (15.29). 

41. The system of claim 39, wherein a first MO and its HOs comprising a linear combination 
of bond orders and a second MO that shares a HO with the first. In addition to the mutual 
HO, the second MO comprises another AO or HO having a single bond order or a mixed 
bond order, in order for the two MOs to be energy matched, the bond order of the second 
MO and its HOs or its HO and AO is a linear combination of the terms corresponding to 
the bond order of the mutual HO and the bond order of the independent HO or AO, and in 
general, E T [atom - atom, msp 3 ) , the energy change of each atom msp 3 shell with the 

formation of each atom-atom-bond MO, is a weighted linear sum for different values of s 
that matches the energy of the bonded MOs, HOs, and AOs: 

N 

E T [atom - atom, msp 3 ) = ^ c Sn E T [atom - atom [s n ) , msp 3 ) (15.30) 
where c s is the multiple of the BO of s n . The radius r 3 of the atom msp 3 shell of each 
bonding atom is given by the Coulombic energy using the initial energy E Coulomb [atom, msp 3 ) 
and E T [atom - atom, msp 3 ) , the energy change of each atom msp 3 shell with the formation 



of each atom-atom-bond MO: 

f* 3 — ~^ _ (15 31) 

msp 8ft£ Q a 0 [E Coulonb [atom, msp 3 ) + E T [atom - atom, msp 3 )) 

where E Coulomb [C2sp 3 ) = -14.825751 eV . 



42. The system of claim 41, wherein the Coulombic energy E Couhmb [molatom.msp 3 ^ of the 
outer electron of the atom msp 3 shell is given by Eq. (15.19). 

43. The system of claim 42, wherein in the case that during hybridization, at least one of the 
spin-paired AO electrons is unpaired in the hybridized orbital (HO), the energy change for 
the promotion to the unpaired state is the magnetic energy E {magnetic) (Eq. (15.15)) at 

the initial radius r of the AO electron. Then, the energy E[molatom,msp 3 } of the outer 
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electron of the atom msp 3 shell is given by the sum of E Coulomb {molatom, msp 3 ) and 
E (magnetic) (Eq. ( 1 5 .20)) . 

44. The system of claim 42, wherein, E T [atom - atom, msp 3 ) , the energy change of each 

atom msp 3 shell with the formation of the atom-atom-bond MO is given by the difference 
between E{mol.atom,msp 3 } and E{atom,msp 3 ) given by Eq. (15.21). 



45. The system of claim 44, wherein using the equation (15.23) for E Coulomb (C,2sp 3 ) in 

equation (15.31), the single bond order energies given by Eqs. (15.1 8-15.27) and shown in 
Table 15.1, and the linear combination energies (Eqs. (15.28-15.30)), the parameters of 
iinear combinations of bond orders and linear combinations of mixed bond orders are 
given in Table 15.2, wherein the calculated and measured values and constants recited in 
the equations herein can be adjusted, for example, up to ± 10%, if desired. 

46. The system of claim 29, wherein the radius r moJ2 ^ of the C2sp 3 HO of a carbon atom of a 
given specie is calculated using Eq. (14.5 14) by considering ]>] E^ [MO, 2sp 3 ) , the total 
energy donation to each bond with which it participates in bonding. 



47. The system of claim 1, wherein equation for the radius is given by 

2 



—e 

r 

mol2sp- 



8^o {E Coulomb (C^p 3 ) + 2X o( (MO,2sp 3 )) 



e 2 



(15.32). 



8^ 0 (el4.825751 eV + ^E^ (MO,2sp 3 )\) 



48. The system of 29, wherein the Coulombic energy E Cmlomb [mol.atom,msp^ of the outer 
electron of the atom msp 3 shell is given by Eq. (15.19). 



84 



WO 2007/053486 



PCT/US2006/042140 



49. The system of claim 29 5 in the case that during hybridization, at least one of the spin- 
paired AO electrons is unpaired in the hybridized orbital (HO), the energy change for the 
promotion to the unpaired state is the magnetic energy E (magnetic) (Eq. (15.15)) at the 

initial radius r of the AO electron. Then, the energy E {mohatom, msp 3 ) of the outer 

electron of the atom msp 3 shell is given by the sum of E Coulomb (mol.atom^msp 3 ^ and 
E(magnetic) (Eq. (15.20)). 

50. The system of claim 49, wherein for the Clsp 2 HO of each methyl group of an alkane 
contributes -0.92918 eV (Eq. (14.513)) to the corresponding single C-C bond; the 
corresponding C2sp 3 HO radius is given by Eq. (14.514). The C2sp 3 HO of each 
methylene group of C n H 2}1 + 2 contributes -0.92918 eV to each of the two corresponding 
G — C bond MOs, wherein the calculated and measured values and constants recited in 
the equations herein can be adjusted, for example, up to ± 10%, if desired. 

51. The system of claim 50, wherein the radius (Eq. (15.32)), the Coulombic energy (Eq. 
(15.19)), an d the energy (Eq. (15.20)) of each alkane methylene group are 

r 

alkaneC melhy!e , ie 2sp 



8^o {E Cou!omb {C, 2sp 3 ) + Y i E Tamme {methylene C-C, 2sp 3 )) 



e 2 



(15.33) 



&re 0 (el 4.825751 eF+e0.92918 eF + e0.92918 eV) 
= 0.81549a 0 



^4^^)- 8 , go(o l 2 549 fl0 ) = - 16 - 68412 eV (1534) 

E{C methyhm 2sp>)= ~f + 2 **£* , 3 =-16.49325 eV (15.35) 

v ' ' 8;zz? o (O.81549tf 0 ) m] (0.843 17a 0 ) 

wherein the calculated and measured values and constants recited in the equations 
herein can be adjusted, for example, up to + 10%, if desired. 

52. The system of claim 51, wherein in the determination of the parameters of functional 
groups, heteroatoms bonding to C2sp 3 HOs to form MOs are energy matched to the 
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C2sp 3 HOs, the radius and the energy parameters of a bonding heteroatom are given by 
the same equations as those for C2sp 3 HOs. 

53. The system of claim 52, wherein using Eqs. (15.15), (15.19-15.20), (15.24), and (15.32) 
in a generalized fashion, the final values of the radius of the HO or AO, r AtomJJOAO , 

E c ^ emb (molMtom 9 m^) 9 and E(C mol 2sp 3 ) are calculated using £^(A»,2jp 3 ), 
the total energy donation to each bond with which an atom participates in bonding 
corresponding to the values of E T - C, C2$p 3 J of the MO due to charge donation 

from the AO or HO to the MO given in Tables 15.1 and 15.2 and the final values of 
r Atom . H0 .Ao> E couio,n b (molMom,msp 3 ), and E(C mol 2sp 3 ) calculated using the values of 

E r ^C B °C,C2sp 3 ^j given in Tables 15.1 and 15.2 are shown in Tables 15.3Aand 15.3B 

in the specification, wherein the calculated and measured values and constants recited in 
the equations herein can be adjusted, for example, up to ± 10%, if desired. 

54. The system of claim 53, wherein the energy of the MO is matched to each of the 
participating outermost atomic or hybridized orbitals of the bonding atoms wherein the 
energy match includes the energy contribution due to the AO or HO's donation of charge 
to the MO. 

55. The system of claim 54, wherein the force constant k % (Eq. (15.1)) is used to determine 
the ellipsoidal parameter c' (Eq. (15.2)) of the each H 2 -type-ellipsoidal-MO in terms of 
the central force of the foci, c' is substituted into the energy equation (from Eq. (15.11))) 
which is set equal to ^ times the total energy of H 2 where ^ is the number of equivalent 
bonds of the MO and the energy of H 2 , -31.63536831 eV , Eq. (11.212) is the minimum 
energy possible for a prolate spheroidal MO, wherein the calculated and measured values 
and constants recited in the equations herein can be adjusted, for example, up to ± 10%, if 
desired. 
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56. The system of claim 55, wherein the energy equation and the relationship between the 
axes, the dimensions of the MO are solved, the energy equation has the semimajor axis a 
as it only parameter, the solution of the semimajor axis a allows for the solution of the 
other axes of each prolate spheroid and eccentricity of each MO (Eqs. (15.3-15.5)), and 
the parameter solutions then allow for the component and total energies of the MO to be 
determined. 



57. The system of claim 1, wherein the total energy, E T {h 2 mo) , is given by the sum of the 
energy terms (Eqs. (15.6-15.1 1)) plus E T (AO I HO) : 
E T {h.mo) = V e +T + V m +V p +E T (AO I HO) (15.36) 



E T [h 2 mo) = ■ 



87rsJa 2 -b 2 



c x c 2 2-^ 



In , -1 



a ) a-^fa 





( 

2- 




In 


a + c' j 




a j 




a—c' 



2 b 2 

+ E T (AO/HO) 



+ E T (AO/HO) 



(15.37) 

where i\ is the number of equivalent bonds of the MO, c x is the fraction of the H 2 -type 
ellipsoidal MO basis function of a chemical bond of the group, c 2 is the factor that results in 
an equipotential energy match of the participating at least two atomic orbitals of each 
chemical bond, and E T (AO I HO) is the total energy comprising the difference of the energy 

E(AOI HO) of at least one atomic or hybrid orbital to which the MO is energy matched and 
any energy component AE H2M0 (AO/HO) due to the AO or HQ's charge donation to the MO. 



58. The system of claim 57, wherein 

E T (AO/HO) = E(AO/HO)-AE H2MO (AO/HO) (15.38) 



59. The system of claim 58, wherein as specific examples, E T (AO I HO) is one from the 
group of 

E T (AO /HO) = E(02p shell) = -E (ionization; O) = -13.6181 eV ; 
E T (AO I HO) = E(N2p shell) = ~E(ionization; N) = -14.53414 eV ; 
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E T {AOIHO 
E T (AO/HO 
E T (AO/HO 
E T (AO/HO 
E T (AO/HO 
E T (AO/HO 
E T (AOIHO 
E T (AOIHO 
E T {AO I HO 
E T (AOJHO 
E T (AO/HO 



= E(C,2sp 3 ) = -14.63489 eV ; 



: E Coulomb (Cl,3sp 3 ) = -14.60295 eV ; 



■ E(ionization; C) + EQonization; C + ) ; 
■E(C eiham ,2sp 3 ) = -15.35946 eF; 

: E(C, 2sp 3 ) - 2E T (C = C, 2sp 3 ) = -14.63489 eV -(-2.26758 eF) ; 
: ^(<^^ 

■ E(C,2sp 3 )-2E T (c = C,2sp 3 ) = -14.63489 eF -(-3.13026 eF); 
; ^(C^,2jp 3 )-£(c^ 9 2sp 3 ) ; 

£(C 3 2^ 3 )-£ r (C - C,2^ 3 ) = -14.63489 eV - (-1.13379 eF) , and 
E [C alkam , 25p 3 ) =- 1 5.56407 e V 9 wherein the calculated and measured 

values and constants recited in the equations herein can be adjusted, for example, up to + 
10%, if desired. 



> 60. The system of claim 1, wherein to solve the bond parameters and energies, 



c' = a. 



Tl 2 47V€ ( 



0 __ 



m e e 2 2C x C 2 a \2C X C 2 



(Eq. (15.2)) is substituted into E T {h 2 mo) to give 



Jb T \H z MO) = 1= 



*) a-4a 2 -b 2 



■E T (AOIHO) 



%7CS Q C ' 



c \ c 2 \ 2 ~ In :"1 



)2C l C 2 



a-c* 



+ E T {AOIHO) 



c,c,\ 2-^-|ln- 
a 



a+ i 




2C,C 2 






2C t C 2 



-1 



+ 



E T (AOIHO) 



(15.39) 

wherein the total energy is set equal to E(basis energies) which in the most general case is 
given by the sum of a first integer ^ times the total energy of H 2 minus a second integer n 2 
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times the total energy of H , minus a third integer n z times the valence energy of E(AO) 
(e.g. E{N) = -14.53414 eV) where the first integer can be 1,2,3,.., and each of the second 
and third integers can be 0 5 1, 2, 3 . . . . 

61. The system of claim 60, wherein 

E(basis energies) = ^ (-31.63536831 eV)-n 2 (-13.605804 eV)-n 3 E(AO) 
(15.40) 

in the case that the MO bonds two atoms other than hydrogen, E(basis energies) is y\ times 
the total energy of H 2 where 7\ is the number of equivalent bonds of the MO and the energy 
of H 2 , -31 .6353683 1 eV , Eq. (1 1 .212) is the minimum energy possible for a prolate 
spheroidal MO: 

E(basis energies) = n, (-31.63536831 eV) (15.41) 

E T {h 2 mo) , is set equal to E(basis energies) , and the semimajor axis a is solved. 

62. The system of claims 23 or 60, wherein the semimajor axis a is solved from the equation 
of the form: 



r^e 2 



%7te, 1 - 



* xl 2C x C 2 



aa 



'0 



^12-^lln V C ^-1 



] 2C X C 2 



+ E T (AO/ HO) = E (basis energies) 



(15.42) 

The distance from the origin of the H 2 -type-ellipsoidal-MO to each focus c ' , the internuclear 
distance 2c ' ? and the length of the semiminor axis of the prolate spheroidal H 2 -type MO 
b = c are solved from the semimajor axis a using Eqs. (15.2-15.4). 

63. The system of claim 62, wherein the component energies are given by Eqs. (15.6-15.9) 
and (15.39). 

64. The system of claim 57, wherein the total energy of the MO of the functional group, 
E T {mo) , is the sum of the total energy of the components comprising the energy 
contribution of the MO formed between the participating atoms and 
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E T {atom - atom, msp 2 .AO^ , the change m the energy of the AOs or HOs upon forming 
the bond. 

65. The system of claim 64, wherein from Eqs. (15.39-15.40), E T {mo) is 

E T (mo) = E(basis energies) + E T {atom - atom, msp* AO} (15.43) 

66. The system of claim 65, wherein during bond formation, the electrons undergo a reentrant 
oscillatory orbit with vibration of the nuclei, and the corresponding energy E osc is the sum 

of the Doppler, E D , and average vibrational kinetic energies, E Kvib : 



^osc=n l {E D ^E Kvib )^n i 



f 



v \m e c 2 yjLi; 



(15.44) 



where j\ is the number of equivalent bonds of the MO , k is the spring constant of the 
equivalent harmonic oscillator, and ju is the reduced mass. 

67. The system of claim 66, wherein the angular frequency of the reentrant oscillation in the 
transition state corresponding to E D is determined by the force between the central field 
and the electrons in the transition state. 



68. The system of claim 67,wherein the force and its derivative are given by 

/ w-- e -*^ (is - 45) 

and 

r r <, 2 

(15.46) 



f'(a) = 2c B Cl ° C *° e2 



"BO 



such that the angular frequency of the oscillation in the transition state is given by 



CD = — 



-3 
a 



/(*)-/■(") 



C C e 2 



*™« K - (15.47) 



where R is the semimajor axis a or the semiminor axis b depending on the eccentricity of 
the bond that is most representative of the oscillation in the transition state, c B0 is the bond- 
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order factor which is 1 for a single bond and when the MO comprises 7\ equivalent single 
bonds as in the case of functional groups, c BO is 4 for an independent double bond as in the 
case of the C0 2 and N0 2 molecules and 9 for an independent triplet bond, C 1o is the 
fraction of the H 2 -type ellipsoidal MO basis function of the oscillatory transition state of a 
chemical bond of the group, and C 2o is the factor that results in an equipotential energy match 
of the participating at least two atomic orbitals of the transition state of the chemical bond. 
Typically, C lo = C x and C 2o = C 2 , the kinetic energy, E K , corresponding to E D is given by 
Planck's equation for functional groups: 



69. The system of claim 68, wherein the Doppler energy of the electrons of the reentrant orbit 




4xs 0 R 3 



(15.48) 



is 




(15.49) 



E osc given by the sum of E D and E Kvib is 



E, 



osc 




(15.50) 



J 



E hv of a group having i\ bonds is given by E T (mo)/^ such that 




Et+osc ( g ™p) * s B iven by the sum of E T (mo) (Eq. (15.42)) and E osc (Eq. (15.51)): 
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E T+osc ( G ™P) = E T ( M °) + E o* 

rr 



\2C,C 2 



c,c 2 2—2- In 



a 



a + <-^- 



a- 



aa 



0 



2C,C 2 



(.4 O I HO) + E T [atom - atom, ww/? 3 .,40) 



2h\ 



\ C lo C 2o e 2 
4ns„R 3 



= [E (basis energies) 4- i? r (afo/w - atom, msp^.AO^ 




(15.52). 



70. The system of claim 69, wherein the total energy of the functional group E T (group) is the 
sum of the total energy of the components comprising the energy contribution of the MO 
formed between the participating atoms, E(basis energies) 9 the change in the energy of 

the AOs or HOs upon forming the bond ( E T [atom - atom, msp 3 .AO^ ), the energy of 
oscillation in the transition state, and the change in magnetic energy with bond formation, 



mag * 



71. The system of claim 70, wherein from Eq. (15.52), the total energy of the group E T (Group) 
is 



E*p {Group) = 



(15.53). 



{E(basis energies) + E T [atom - atom, msp^.Acty 
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72. The system of claim 71, wherein the change in magnetic energy E mag which arises due to 
the formation of impaired electrons in the corresponding fragments relative to the bonded 
group is given by 



E mag ~ C 3 2 3 ~ C 3 

mzr 



(15.54) 



where r 3 is the radius of the atom that reacts to form the bond and c 3 is the number of 

electron pairs. 

r 



Ej, (Group) — 



(15.55). 



(E(basis energies) + E T [atom — atom, msp 3 



1 + 



i 



2h\ 



C C e 



mc 



73 . The system of claim 72, wherein the total bond energy of the group E D (Group) is the 
negative difference of the total energy of the group (Eq. (15.55)) and the total energy of 
the starting species given by the sum of c A E initial (c 4 AO I HO) and c 5 E mml (c 5 AO I HO) : 



E D (Group) — — 



{^E(basis energies) + E T {atom - atom, msp^.AO^ 



1 + 



2h \ 


C lo C,„e 2 


4tt£ o R 3 




m e c 2 



+n x E Kvib + c 3 ^^-(c A E mtial {AO I HO) + c s E miial (c 5 AO I HO)) 



(15.56). 



74. The system of claim 70, wherein in the case of organic molecules, the atoms of the 
functional groups are energy matched to the C2sp 3 HO such that 
E(AO/HO) = -14.63489 eV (15.57) 
^0 For examples of E mag from previous sections: 
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E mag (C2sp>) = c 3 ^ = c 3 (Q 8 9 ^ )3 =C3 0.14803 eV (15.58) 
E mog {02p) = c 3 = c 3 ^£ = C3 0.1 1441 eV (15.59) 



r a n 

2 



E ^ {N2p) = Ci ^, c -e,o. mK .y (15.60). 

""^ ^ 3 r 3 3 (0.93084a 0 ) 3 



75. The system of claim 74, wherein in the general case of the solution of an organic 

functional group, the geometric bond parameters are solved from the semimajor axis and 
the relationships between the parameters by first using Eq. (15.42) to arrive at a , the 
remaining parameters are determined using Eqs. (15.1-15.5), the energies are given by 
Eqs. (15.52-15.59), and to meet the equipotential condition for the union of the H 2 -type- 
ellipsoidal -MO and the HO or AO of the atom of a functional group, the factor c 2 of a 
H 2 -type ellipsoidal MO in principal Eqs. (15.42) and (15.52) may given by 

(i) one: 

c 2 =l (15.61) 

(ii) the ratio that is less than one of 13.605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1.243), and the magnitude of the 
Coulombic energy of the participating AO or HO of the atom, E Coulomb {MO.atom^msp 2 ) 

given by Eqs. (15.19) and (15.31-15.32). For \E Coulomb (MO.a^ 5 w^ 3 )|>13.605804 eV : 



e 2 



87V£ 0 a 0 _ 13.605804 eV 



e 2 



For 



^Coulomb 

(MO.atom s msp 3 )\ 

Eco Ulon ,t (MO.atom,msp*)\ <13.605804 eV : 
e 2 

%7ES * r A-BAor Bsp > _ \ E cou lo ,r, b {MO.atom,m S p 3 )\ 



(15.62) 



„2 



13.605804 eV 



(15.63) 



8jre 0 a 0 
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(iii) the ratio that is less than one of 13.605804 eV , the magnitude of the Coulombic 
energy between the electron and proton of H given by Eq. (1 .243), and the magnitude of the 
valence energy, Eiyalence) , of the participating AO or HO of the atom where Eiyalence) is 

the ionization energy or E^MO.atom, msp 3 ) given by Eqs. (15.20) and (15.31-15.32). For 
\E(valence)\ >13.605804 eV : 



e 1 



e 



%K S 0 r A -BAorBzi? 



\E(yalence)\ 



For \E(yalence)\< 13. 605804 eV : 



e 2 



J-^ ^\Eiyalenaei 
2 e 2 13.605804 eV 



(iv) the ratio that is less than one of the magnitude of the Coulombic energy of the 
participating AO or HO of a first atom, E Coulomb (MO. atom, msp 3 } given by Eqs. (15.19) and 
(15.31-15.32), and the magnitude of the valence energy, Eiyalence) 9 of the participating AO 
or HO of a second atom to which the first is energy matched where E{valence) is the 
ionization energy or E [MO.atom, msp 3 ) given by Eqs. (15.20) and (15.31-15.32). For 

\ E cauh»b {MO.atom, msp 3 )| > E{valence) : 
\E(yalence)\ 



(15.66) 



For 



Ecouiomb (MO.atom, msp 3 )| 
E couhmb {MO,atom,msp 3 ^ <E(yalence) : 

_ \Ecouio,nb (MO.atom,msp 3 )\ (156?) 
\E(valence)\ 



(v) the ratio that is less than one of the magnitude of the valence-level energies, 
E n (yalence) , of the AO or HO of the nth participating atom of two that are energy matched 
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where E(yalence) is the ionization energy or E^MO.atom^msp 3 ^ given by Eqs. (15.20) and 

(15.31-15.32): 

^(valence) 

E 2 (yalence) 

(vi) the factor that is the ratio of the hybridization factor c 2 (l) of the valence AO or 
HO of a first atom and the hybridization factor c 2 (2) of the valence AO or HO of a second 
atom to which the first is energy matched where c 2 (n) is given by Eqs. (15.62-15.68); 
alternatively c 2 is the hybridization factor c 2 (l) of the valence AOs or HOs a first pair of 
atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 
pair to which the first two are energy matched: 

05,9) 

(vii) the factor that is the product of the hybridization factor c 2 (l) of the valence AO 
or HO of a first atom and the hybridization factor c 2 (2) of the valence AO or HO of a second 
atom to which the first is energy matched where c 2 (n) is given by Eqs. (15.62-15.69); 
alternatively c 2 is the hybridization factor c 2 (l) of the valence AOs or HOs a first pair of 
atoms and the hybridization factor c 2 (2) of the valence AO or HO a third atom or second 
pair to which the first two are energy matched: 

c 2 =c 2 (l)c 2 (2) (15.70) 

The hybridization factor c 2 corresponds to the force constant k (Eqs. (1 1.65) and (13.58)). 

In the case that the valence or Coulombic energy of the AO or HO is less than 13.605804 eV 5 
the magnitude of the Coulombic energy between the electron and proton of H given by Eq. 
(1.243), then C 2 corresponding to k 1 (Eq. (15.1)) is given by Eqs. (15.62-15.70). 

76. The system of claim 75, wherein specific examples of the factors c 2 and C 2 of a H 2 -type 
ellipsoidal MO of Eq. (15.51) given in following sections are 
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c 2 (C2sp*HO to F) = E(C J^K (C2sp*HO) = ~ 14 - 634 * 9 gF (Q.9177l) = 0.77087 ; 
x ; E(F) 2V J -17.42282 eV K ' 

C 2 (C2sp 3 HO to Cl)= J} Cl \. c 2 (C2 S p 3 HO) = - l2 - 96764eV (0.91771) = 0.81317 ; 
v ' E[C,2sp 3 ) 2 \ ' -14.63489 eV K } 

c,(C2*HO ,o 0 )__^^(cV«>)-Ej^£ (C.9 17 71 M . 8 5395; 

= *(jfl- - - 14 - 534 ' 4 ^ . 0.94627 ; 
' £(C,2sp 3 ) -15.35946 eV 

v J i'fc^sp*) -15.56407«f 

C 2 (a, » *) = *i£M . -'"^W ^ , 0.76144 ; 
v ' -13.60580 sV 

C 2 (cisp'HO to S) = -^fL^ (C2jp"/fO) = Zj|j^57^(0-91771) = 0.64965 ; 

c, (O to S3sp' to Clsp'HO) = c 2 (C2sp s flO) 

-13.61806 ,r (0 91771) . 



-10.36001 
= 1.20632 



/«-, 3\ £ c«fc«A (-S^^ 3 ) -11.57099 
^ )" J( g ) = -13.60580 e F = °- 85Q45; 

C 2 (C2sp?HO to S3 S p> ) ~ (S3sp> U - 11 ' 52126 gF (0.85045) 

V 7 £(c,2*p 3 ) n ' -14.63489 eF V ; 
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C 2 (S3sp 3 to O to C2sp 3 HO) = ^^^c 2 (C2sp 3 HO) 

-11.52126 eV 



-13.61806 eV 
= 0.77641 



(0.91771) ; 



c 2 (O to N2p to C2sp 3 HO) = c 2 (C2sp 3 HO) 



13.61806 eV {0 91J7l) 



-14.53414 eV 
= 0.85987 



,„„ . c 2 (C2sp 3 HO to N) 0 91140 

c 2 (N2p to 02p) = -±± - ) uvvvw 06727- 

2V V} c 2 (C2sp 3 HOtoO) 0.85395 > 

C 2 (benzeneC2sp 3 HO) = c 2 ( W e «eC2sp 3 M?) = 13 - 605804 eV = 0>85252 

15*95955 ^ 1^ 

c 2 (arylC2sp 3 HO to 0)= f^°\. c 2 (arylC2sp 3 HO) 

EiC,2sp J 7 

= -13.61806 eF 

-14.63489 eF V ; 
= 0.79329 

, 2 (H to anline JV) = ?W = -14.53414 

V ' £(c,2sy) -15.76868 ^ Z1/i > 

(arylC2sp 3 HO to N)= (arylC2sp 3 HO) 

_ -14.53414 «K , v 
~ -14.63489 e F ( °- 85252) 
= 0.84665 

/ej-, v »r\ -14.53414 
<? 2 (i7 to anlme iV r )= , v = 0.92171 : 

7 £(C,2.sp 3 ) -15.76868 eV 



c 2 



c 2 (arylC2sp 3 HO to N)= c 2 (arylC2sp 3 HO) 

-14.53414 _ . 

= 0.84665 

C 2 (Bp * axyMype C2*»flo)- ^4 - - 10 ' 36001 eK = 0.65700 . 
v 7 £(C,2sp 3 ) -15.76868 
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wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10%, if desired. 



77. The system of claim 1, wherein the parameters of the point of intersection of each H 2 - 
type ellipsoidal MO and the A -atom AO are determined from the polar equation of the 
ellipse: 

1 + e 



r = r r 



o l + ecos0' 



(15.71). 



78. The system of claim 77, wherein the radius of the A shell is r A , and the polar radial 
coordinate of the ellipse and the radius of the A shell are equal at the point of intersection 
such that 



r A ={a-c*)- 



1 + - 



1 + — cos#' 
a 



(15.72). 



79. The system of claim 78 5 wherein the polar angle 6 ' at the intersection point is given by 



f f 



<9' = cos" 1 



a 



(15.73). 



J J 



80. The system of claim 79, wherein the angle 0 AAO the radial vector of the A AO makes 
with the internuclear axis is 



(15.74). 



81. The system of claim 80, wherein the distance from the point of intersection of the orbitals 
to the internuclear axis is the same for both component orbitals such that the angle 
cot = d„ ,,n between the internuclear axis and the point of intersection of each H 2 -type 
ellipsoidal MO with the A radial vector obeys the following relationship: 



(15.75) 



such that 
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Ow='to* r '* m ° AM (15-76). 

82. The system of claim 81, wherein the distance d H ^ f0 along the internuclear axis from the 
origin of H 2 -type ellipsoidal MO to the point of intersection of the orbitals is given by 

= aoosO HiMO (15.77). 

83. The system of claim 82, wherein the distance d AAO along the internuclear axis from the 
origin of the A atom to the point of intersection of the orbitals is given by 

d AA0 =c^d HiM0 (15.78). 

84. The system of claim 1, wherein in ACB MO comprising a linear combination of C - A - 
bond and C - B -bond MOs where C is the general central atom and a bond is possible 
between the A and B atoms of the C-A and C-B bonds, the ZACB bond angle is 
determined by the condition that the total energy of the H 2 -type ellipsoidal MO between 
the terminal A and B atoms is zero. 

85. The system of claim 84, wherein the force constant k ' of a H 2 -type ellipsoidal MO due 
to the equivalent of two point charges of at the foci is given by: 

C C 2e 2 

v ^iw (1579) 

Atvs 0 

where C x is the fraction of the H 2 -type ellipsoidal MO basis function of a chemical bond of 
the specie which is 0.75 (Eq. (13.59)) for a terminal A-H {A is H or other atom) and 1 
otherwise and C 2 is the factor that results in an equipotential energy match of the 
participating at least two atomic orbitals of the chemical bond and is equal to the 
corresponding factor of Eqs. (15.42) and (15.52). 

86. The system of claim 85, wherein the distance from the origin of the MO to each focus c ' 
of the A-B ellipsoidal MO is given by: 

e , / h 2 4„ So = r^r 



100 



WO 2007/053486 



PCT/US2006/042140 



87. The system of claim 85, wherein the intemuclear distance is 



ci5 - 8,) - 



88. The system of claim 85, wherein the length of the semiminor axis of the prolate 
spheroidal A-B MO b = c is given by Eq. (15.4). 

89. The system of claim 85, wherein the component energies and the total energy, E T (h 2 mo) , 
of the A-B bond are given by the energy equations (Eqs. (1 1 .207-1 1 .212), (1 1 .213- 
11.217), and (11.239)) of H 2 except that the terms based on charge are multiplied by c BO , 
the bond-order factor which is 1 for a single bond and when the MO comprises ^ 
equivalent single bonds as in the case of functional groups. c BO is 4 for an independent 
double bond as in the case of the C0 2 and N0 2 molecules. 

90. The system of claim 85, wherein the kinetic energy term is multiplied by c\ Q which is 1 
for a single bond, 2 for a double bond, and 9/2 for a triple bond, the electron energy terms 
are multiplied by c x , the fraction of the H 2 -type ellipsoidal MO basis function of a 
terminal chemical bond which is 0.75 (Eq. (13.233)) for a terminal A-H {A is H or 
other atom) and 1 otherwise. 

91 . The system of claim 85, wherein the electron energy terms are further multiplied by c' 2 , 
the hybridization or energy-matching factor that results in an equipotential energy match 
of the participating at least two atomic orbitals of each terminal bond. 

92. The system of claim 84, wherein when A-B comprises atoms other than H , 

E T [atom - atom,msp 3 .AO) , the energy component due to the AO or HO's charge 
donation to the terminal MO, is added to the other energy terms to give E T [h 2 mo) : 



E T {hmo) = 



87C£ 0 C ' 




(15.82). 
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93. The system of claim 84, the radiation reaction force in the case of the vibration of A - B 
in the transition state corresponds to the Doppler energy, E D 9 given by Eq. (11.181) that 
is dependent on the motion of the electrons and the nuclei. 

94. The system of claim 93 , wherein the total energy that includes the radiation reaction of the 
A-B MO is given by the sum of E T (h 2 mo) (Eq. (15.82)) and E osc given Eqs. (1 1 .213- 

1 1.220), (1 1.23 1-1 L236), and (1 1.239-1 1.240). 



95. The system of claim 94, wherein the total energy E T (A-B) of the A-B MO including 
the Doppler term is 



E T (A-B) = 



a J a — c 



+ E T {atom — atom, msp 3 .AO^j 



2%\ 



'BO 



C la C 2o e 2 



1+' 



m a 



m n c 



cA e2 



'BO 



(15.83) 

where C 1o is the fraction of the H 2 -type ellipsoidal MO basis function of the oscillatory 
transition state of the A - B bond which is 0.75 (Eq. (13.233)) in the case of H bonding to a 
central atom and 1 otherwise, C 2o is the factor that results in an equipotential energy match of 
the participating at least two atomic orbitals of the transition state of the chemical bond, and 



_ m x ?n 2 
m 1 + m 2 



is the reduced mass of the nuclei given by Eq. (11.1 54). 



96. The system of claim 95, wherein to match the boundary condition that the total energy of 
the A-B ellipsoidal MO is zero, E T (A - B) given by Eq. (15.83) is set equal to zero and 
substitution of Eq. (15.81) into Eq. (15.83) gives 



102 



WO 2007/053486 



PCT/US2006/042140 



0 = 



8** 



'2QC 2 



2c fl0 -c' B0 -*-|ln- 



aa 0 
'2C,C 2 



-1 



f 2C,C 2 



■ i? r {atom - atom, msp 3 .AO} 



1 + 
(15.84). 



2h 



i 



m. 



mc 



+ 



j 2 



C BO e 



'BO 



%7ts n a 



\2C X C 2 



97. The system of claim 96, wherein the vibrational energy-term of Eq. (15.84) is determined 
by the forces between the central field and the electrons and those between the nuclei 
(Eqs. (11.141-11.145)). 



98. The system of claim 96, wherein the electron-central-field force and its derivative are 
given by 



J 2 



-£^7 (15-85) 



and 



4/re 0 a 



(15.86). 



99. The system of claim 97, wherein the nuclear repulsion force and its derivative are given 
by 

2 



/(a + c') = 



8/zs 0 (a + c') 



and 



f'(a + c') = ~ 



(15.87) 



(15.88) 



%7ue 0 (a + c') 

such that the angular frequency of the oscillation is given hy 



co = 



f /(«)-/>) 



'BO 



4xe 0 a 3 Sxe 0 (a + c') 



(15.89). 
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1 00. The system of claim 99, wherein since both terms of E osc ~E D + E Kvib are small due to 

the large values of a and c\ an approximation of Eq. (15.84) which is evaluated to 
determine the bond angles of functional groups is given by 



0 = 




+ E T ( atom - atom, msp 2 .A 



(15.90). 



101. The system of claim 1 00, wherein from the energy relationship given by Eq. (15 .90) 
and the relationship between the axes given by Eqs. (1 5.2-1 5.5), the dimensions of the 
A-B MO can be solved. 

1 02. The system of claim 101, wherein Eq. (1 5.90) is solved by the reiterative technique 
using a computer. 

103. The system of claim 100, wherein a factor c 2 of a given atom in the determination of 
c' 2 for calculating the zero of the total A-B bond energy is given by Eqs. (15.62-15.65). 

1 04 . The system of claim 1 03 , wherein in the case of a H - H terminal bond of an alkyl 
or alkenyl group, c' 2 is the ratio of c 2 of Eq. (15.62) for the H- H bond which is one and 
c 2 of the carbon of the corresponding C — H bond: 

c <- 1 13.605804 eV 

2 c % (C2sp^ E Coulomb {C-HC2sp>) 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 
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105. The system of claim 100, wherein in the case of the determination of the bond angle 
of the ACH MO comprising a linear combination of C - A -bond and C-H -bond MOs 
where A and C are general, C is the central atom, and c 2 for .an atom is given by Eqs. 
(15.62-15.70), c' 2 of the A- H terminal bond is the ratio of c 2 of the A atom for the 
A — H terminal bond and c 2 of the C atom of the corresponding C — H bond: 



106. The system of claim 100, wherein in the case of the determination of the bond angle 
of the COH MO of an alcohol comprising a linear combination of C - O -bond and 
O—H -bond MOs where C 5 O, and H are carbon, oxygen, and hydrogen, respectively, 
c 2 of the C-H terminal bond is 0.91771 since the oxygen and hydrogen atoms are at 
the Coulomb potential of a proton and an electron (Eqs. (1.236) and (10.162), 
respectively) that is energy matched to the C2sp 3 HO, wherein the calculated and 
measured values and constants recited in the equations herein can be adjusted, for 
example, up to ± 10%, if desired. 

107. The system of claim 100, wherein in the determination of the hybridization factor c 2 
of Eq. (15.90) from Eqs. (15.62-15.70), the Coulombic energy, E Coulomb {MO.atom,msp 3 ) , 
or the energy, E [MO.atom, msp 3 ) , the radius r A _ B AorBsp , of the A or B AO or HO of the 
heteroatom of the A~B terminal bond MO such as the C2sp 3 HO of a terminal C-C 
bond is calculated using Eq. (1 5.32) by considering ]JjT E Tmot {M0 9 2sp 3 ) , the total energy 
donation to each bond with which it participates in bonding as it forms the terminal bond. 

108. The system of claim 107, wherein the Coulombic energy E Coulomb (^MO.atom.msp 3 ^ of 
the outer electron of the atom msp 3 shell is given by Eq. (15.19). 

109. The system of claim 108, wherein in the case that during hybridization, at least one of 
the spin-paired AO electrons is unpaired in the hybridized orbital (HO), the energy change 
for the promotion to the unpaired state is the magnetic energy E(magnetic) (Eq. (15.15)) 



c 2 {A{A~H)msp 3 ) 
c 2 (C(C-H)(msp 3 ) 



(15.92). 
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at the initial radius r of the AO electron, and the energy E^MO.atom, msp 3 ) of the outer 
electron of the atom msp* shell is given by the sum of E CouJomb [MO.atom,msp 3 ^ and 
EQnagnetic) (Eq. (15.20)). 

1 10. The system of claim 99, wherein in the specific case of the terminal bonding of two 
carbon atoms, the c 2 factor of each carbon given by Eq. (15.62) is determined using the 

Coulombic energy E Coulomh {C-C C2sp 3 ) of the outer electron of the C2sp 3 shell given 

byEq. (15.19) with the radius r ccc2sp2 of each C2sp 3 HO of the terminal C-C bond 

calculated using Eq. (1 5 .32) by considering ]T E^ [MO, 2sp 3 ) , the total energy donation 

to each bond with which it participates in bonding as it forms the terminal bond including 
the contribution of the methylene energy, 0.92918 eV (Eq. (14.513)), corresponding to 
the terminal C-C bond. The corresponding E T (atom -atom, msp 3 .AO) in Eq. (15.90) is 
E T (C-C C2sp 3 ) = -1 .85836 eV , wherein the calculated and measured values and 

constants recited in the equations herein can be adjusted, for example, up to + 10%, if 
desired. 

111. The system of claim 1 00, wherein in the case that the terminal atoms are carbon or 
other heteroatoms, the terminal bond comprises a linear combination of the HOs or AOs; 
thus, c 2 is the average of the hybridization factors of the participating atoms 
corresponding to the normalized linear sum: 

°2 ^{c f 2 {atom l) + c' 2 (atom 2)) (15.93) 

In the exemplary cases of C-C , 0-0,and N-N where C is carbon: 



f e 2 ^ 



87re 0 a 0 $7re 0 a Q 



e 2 e 2 



^^cFa-aa^aoiho ^^a-a k^aoiho J 
r 13.605804 eV 13.605804 eV 



(15.94) 



I Eco Ulomb (A - AA x AO I HO) E Coulomb {A - A^AOIHO) 
In the exemplary cases of C-N , C-O, and C-S, 
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r 




2{E Coulami (C-BC2sp>) 



1 13.605804 eV 



+ c 2 (C toB) 



(15.95) 



where C is carbon and c 2 (C to B) is the hybridization factor of Eqs. (15.52) and (15.84) that 
matches the energy of the atom B to that of the atom C in the group, 

112. The system of claim 97, wherein the corresponding E T {atom - atom, msp* .AO} term 

in Eq. (15.90) depends on the hybridization and bond order of the terminal atoms in the 
molecule, but typical values matching those used in the determination of the bond 
energies (Eq. (15.56)) are 

E T (C-0 C2sp\02p) = -1.44915 eV ; E T (C- O C2sp\02p) - -1.65376 eV ; 
E T (C-N C2sp\N2p) = -1.44915 eV ; E T (c~8 C2sp\S2p) = -0.72457 eV ; 
E T (0-0 02 pD2 p) = -l .449 1 5 eV ; E T (O - O 02 p.02p) = -1.65376 eV; 
E T (N~ N N2p.N2p) = -1.44915 eV ; E T (N-ON2p.02p) = -1.44915 eV; 
E T (F-F F2p.F2p) - -1 .44915 eV ; E T (CI - CI C13 p.CB p) = -0.9291 8 eV ; 
E T (Br - Br Br 4 p3r4 p) = -0.9291 8 e V ; E T (l - J 15 p.I5 p) = -0.36229 eV ; 
E T (C-F C2sp\F2p) = -1.85836 eV ; E T (c -CI C2sp\Cl3p) = -0.92918 eV; 
E T (C - Br C2sp\Br4p) « -0.72457 eV ; E T (C - / C2sp 3 .I5p) = -0.36228 eV , and 



constants recited in the equations herein can be adjusted, for example, up to + 10% 5 if desired. 

113. The system of claim 1 00, wherein in the case that the terminal bond is X - X where 
X is a halogen atom, c x is one, and c\ is the average (Eq. (15.93)) of the hybridization 
factors of the participating halogen atoms given by Eqs. (15.62-15.63) where 



E Coulomb (MO.atom > msp*) -13.605804 eV for X = 7. 
1 14. The system of claim 113, wherein the factor C x of Eq. (15.90) is one for all halogen 



E T (0-CI 02p.Cl3p) = -0.9291 8 eV 



y wherein the calculated and measured values and 



^coulomb {MO.atom,msp % } is determined using Eq. (15.32) and 



atoms. 
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115. The system of claim 113, wherein the factor C 2 of fluorine is one since it is the only 
halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude 
of the Coulombic energy between the electron and proton of H given by Eq. (1 .243), 
wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10%, if desired. 

116. The system of claim 113, wherein for each of the halogens, CI 9 Br , and / , C 2 is 

the hybridization factor of Eq. (15.52) given by Eq. (15.70) with c 2 (l) being that of the 

halogen given by Eq. (1 5.68) that matches the valence energy of X (E t (valence) ) to that 

of the C2sp 3 HO (E 2 (valence) = -14.63489 eV, Eq. (15.25)) and to the hybridization of 

C2sp 3 HO (c 2 (2) = 0.91771 , Eq. (13.430)), wherein the calculated and measured values 

and constants recited in the equations herein can be adjusted, for example, up to + 10%, if 
desired. 

117. The system of claim 1 16, wherein E T [atom - atom, msp 3 .AO) of Eq. (15.90) is the 

maximum for the participating atoms which is -1.44915 eV , -0.92918 eV , 
-0.92918 eV , and -0.33582 eV for F , CI , Br , and / , respectively, wherein the 
calculated and measured values and constants recited in the equations herein can be 
adjusted, for example, up to + 10%, if desired. 

118. The system of claim 100, wherein in the case that the terminal bond is C - X where 
C is a carbon atom and X is a halogen atom, the factors c 1 and C x of Eq. (15.90) are 
one for all halogen atoms. 

119. The system of claim 118, wherein for X = F , c' 2 is the average (Eq. (15.95)) of the 
hybridization factors of the participating carbon and F atoms where c 2 for carbon is 
given by Eq. (15.62) and c 2 for fluorine matched to carbon is given by Eq. (15.70) with 
c 2 (l) for the fluorine atom given by Eq. (15.68) that matches the valence energy of F 
{E i (valence) = -17.42282 eV) to that of the C2sp 3 HO (E 2 (valence) = -14.63489 eV 9 
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Eq. (15.25)) and to the hybridization of C2sp* HO (c 2 (2) = 0.91771, Eq. (13.430)), 

wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10%, if desired. 

120. The system of claim 1 19, wherein the factor C 2 of fluorine is one since it is the only 
halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude 
of the Coulombic energy between the electron and proton of H given by Eq. (1.243), 
wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10% 5 if desired. 

121 . The system of claim 119, wherein for each of the other halogens, CI , Br , and / , c' 2 

is the hybridization factor of the participating carbon atom since the halogen atom is 
energy matched to the carbon atom. 

122. The system of claim 121, wherein the C 2 of the terminal-atom bond matches that 
used to determine the energies of the corresponding C - X -bond MO. 

123. The system of claim 122, wherein C 2 is the hybridization factor of Eq. (15.52) given 
by Eq. (15.70) with c 2 (l) for the halogen atom given by Eq. (15.68) that matches the 
valence energy of X (E x (valence) ) to that of the C2sp 3 HO 

(E 2 (valence) = -14.63489 eV , Eq. (15.25)) and to the hybridization of C2sp* HO 
(c 2 (2) = 0.91771 , Eq. (13.430)), wherein the calculated and measured values and 

constants recited in the equations herein can be adjusted, for example, up to ± 10%, if 
desired. 

124. The system of claim 123, wherein E T (atom ~ atonumsp* .AO) of Eq. (15.90) is the 

maximum for the participating atoms which is -1.85836 eV , -0.92918 eV , 
-0.72457 eV , and -0.33582 eV for F , CI , Br , and / , respectively, wherein the 
calculated and measured values and constants recited in the equations herein can be 
adjusted, for example, up to ± 10%, if desired. 
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125. The system of claim 119, wherein in the case that the terminal bond is H - X 
corresponding to the angle of the atoms HCX where C is a carbon atom and X is a 
halogen atom, the factors c x and C x of Eq. (15.90) are 0.75 for all halogen atoms. 

126. The system of claim 1 19, wherein for X = F , c' 2 is given by Eq. (15.69) with c 2 of 
the participating carbon and F atoms given by Eq. (15.62) and Eq. (15.65), respectively. 

127. The system of claim 126, wherein the factor C 2 of fluorine is one since it is the only 
halogen wherein the ionization energy is greater than that 13.605804 eV , the magnitude 
of the Coulombic energy between the electron and proton of H given by Eq. (1 .243), 
wherein the calculated and measured values and constants recited in the equations herein 
can be adjusted, for example, up to ± 10%, if desired. 

128. The system of claim 127, wherein for each of the other halogens, CI , Br , and / , c' 2 
is also given by Eq. (15.69) with c 2 of the participating carbon given by Eq. (15.62) and 
c 2 of the participating X atom given by c 2 = 0.91771 (Eq. (13.430)) since the X atom is 
energy matched to the C2sp 3 HO. 

129. The system of claim 128, wherein C 2 is given by Eq. (15.65) for the corresponding 
atom X where C 2 matches the energy of the atom X to that of H . 

130. The system of claim 100, wherein the distance between the two atoms A and B of 
the general molecular group ACS when the total energy of the corresponding A-B MO 
is zero, the corresponding bond angle can be determined from the law of cosines: 

+ s 2 2 - 2^»y 2 cosine^ = s 3 2 (15.96). 

131. The system of claim 130, wherein with s x = 2c' c _ A , the internuclear distance of the 
C-A bond, s 2 = 2c' c _ B ? the internuclear distance of each C-B bond, and s 3 = 2c \ 4 _ B , 
the internuclear distance of the two terminal atoms, the bond angle #^ Cfi between the 
C-A and C-B bonds is given by 
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(2c< c _ A f + {2c< c _ B ) 2 -2(2c^)(2c' c _ B )cosine^ = (2c^ B ) 2 (15.97) 

(15.98). 



&ZACB = C0S 1 



(2c' c _,) 2 + (2c^) 2 -(2c^) 2 ^ 



2(2c' c _J(2c' c „ 5 ) 



132. The system of claim 131 ? wherein the structure C b C a (O a )O b wherein C a is bound to 
C b 9 O a , and 0 6 , the three bonds are coplanar and two of the angles are known, say 6 X 
and 0 2 , then the third 0 3 can be determined geometrically: 

d 3 =36O-0 l -0 2 (15.99) 

133. The system of claim 131, wherein in the general case that two of the three coplanar 
bonds are equivalent and one of the angles is known, say 0 1 , then the second and third can 
be determined geometrically: 

02=&3 = y 2 lJ (15.100). 

134. The system of claim 1, wherein in the general case where the group comprises three 
A-B bonds having B as the central atom at the apex of a pyramidal structure formed by 
the three bonds with the A atoms at the base in the xy-plane. 

135. The system of claim 134, wherein the C 3v axis centered on B is defined as the 
vertical or z-axis, and any two A-B bonds form an isosceles triangle, and the angle of 
the bonds and the distances from and along the z-axis are determined from the 
geometrical relationships given by Eqs. (13.412-13.416): 

the distance d origin _ B from the origin to the nucleus of a terminal B atom is given by 

dorizinB^ ^ *"* (15.101) 
on S ,n-B 2sin6Q o ^ ' 

the height along the z-axis from the origin to the A nucleus d height is given by 
the angle 9 V of each A-B bond from the z-axis is given by 
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tan 



origin— B 
\ ^height J 



(15.103). 



136. The system of claim 135, wherein in the case where the central atom B is further 
bound to a fourth atom C and the B - C bond is along the z-axis. Then, the bond 9^^ 
given by Eq. (14.206) is 

^*c= 180-0,. (15.104). 



137. The system of claim 134, wherein in the plane defined by a general AC A MO 
comprising a linear combination of two C - A -bond MOs where C is the central atom, 
the dihedral angle 9 ZBCfACA between the ACA -plane and a line defined by a third bond 
with C 9 specifically that corresponding to a C -B -bond MO, is calculated from the bond 
angle 9^^ and the distances between the A 9 B , and C atoms. 

138. The system of claim 137, wherein the distance d x along the bisector of 6 ZACA from 
C to the internuclear-distance line between A and A , 2c\_ A , is given by 

d x =2c' c „ A cos-^^- (15.105) 

where 2c } C _ A is the internuclear distance between A and C . 

139. The system of claim 138, wherein the atoms A , A , and B define the base of a 
pyramid. 

140. The system of claim 139, wherein the pyramidal angle 9^ BA can be solved from the 
internuclear distances between A and A 5 2c * A _ A , and between A and B , 2c \_ B , using 
the law of cosines (Eq. (15.98)): 



^=cos" 1 



(2c^) 2 + (2c^) 2 ^(2c^) 2 
2(2c' A _ B )(2c* A _ B ) 



(15.106). 
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141. The system of claim 140, wherein the distance d 2 along the bisector of from B 
to the internuclear-distance line 2c\_ A , is given by 



</ 2 = 2^_ 5 cos%^ (15.107). 



142. The system of claim 141 , wherein the lengths d l9 d 2 , and 2c ! c _ 5 define a triangle 
wherein the angle between d x and the internuclear distance between B and C , 2c , c _ 5 , is 
the dihedral angle 0 ABCIACA that can be solved using the law of cosines (Eq. (15.98)): 



&ZBC/ACA ™ COS 



-1 



d?+(2c< c _ B ) 2 ~d : 
2^(2c' s _ c ) 



(15.108). 



143. The system of claim 1, wherein the specie are solved using the solutions of organic 
chemical functional groups as basis elements wherein the structures and energies where 
linearly added to achieve the molecular solutions, each functional group can be treated as 
a building block to form any desired molecular solution from the corresponding linear 
combination, each functional group element was solved using the atomic orbital and 
hybrid orbital spherical orbitsphere solutions bridged by molecular orbitals comprised of 
the H 2 -type prolate spheroidal solution given in the Nature of the Chemical Bond of 
Hydrogen-Type Molecules section, the energy of each MO was matched at the HO or AO 
by matching the hybridization and total energy of the MO to the AOs and HOs, the energy 
E mag (e.g. given by Eq. (15.58)) for a C2sp 3 HO and Eq.(15.59) for an Olp AO) was 
subtracted for each set of unpaired electrons created by bond breakage. 

144. They system of claim 143, wherein the bond energy is not equal to the component 
energy of each bond as it exists in the specie, although, they are close. 

145. The system of claim 143, wherein the total energy of each group is its contribution to 
the total energy of the specie as a whole. 

146. The system of claim 143, wherein the determination of the bond energies for the 
creation of the separate parts must take into account the energy of the formation of any 
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radicals and any redistribution of charge density within the pieces and the corresponding 
energy change with bond cleavage. 

147. The system of claim 143, wherein the vibrational energy in the transition state is 
dependent on the other groups that are bound to a given functional group, which will 
effect the functional-group energy, however because the variations in the energy based on 
the balance of the molecular composition are typically of the order of a few hundreds of 
electron volts at most, they are neglected. 

148. The system of claim 143, wherein the energy of each functional-group MO bonding to 
a given carbon HO is independently matched to the HO by subtracting the contribution to 
the change in the energy of the HO from the total MO energy given by the sum of the MO 
contributions and E(C,2sp 3 ) = -14.63489 eV (Eq. (13.428)). 

149. The system of claim 143, wherein the intercept angles are determined from Eqs. 
(15.71-15.78) using the final radius of the HO of each atom. 

1 50. The system of claim 143, wherein a final carbon-atom radius is determined using Eqs. 
(15.32) wherein the sum of the energy contributions of each atom to all the MOs in which 
it participates in bonding is determined. 

151. The system of claim 1 50, wherein the final radius is used in Eqs. (15.19) and (1 5 .20) 
to calculate the final valence energy of the HO of each atom at the corresponding final 
radius. 

1 52. The system of claim 151, wherein the radius of any bonding heteroatom that 
contributes to a MO is calculated in the same manner, and the energy of its outermost 
shell is matched to that of the MO by the hybridization factor between the carbon-HO 
energy and the energy of the heteroatomic shell. 

153. The system of claim 150, wherein the donation of electron density to the AOs and 
HOs reduces the energy. 
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1 54. The system of claim 143, wherein the donation of the electron density to the MO's at 
each AO or HO is that which causes the resulting energy to be divided equally between 
the participating AOs or HOs to achieve energy matching. 

1 55. The system of claim 1, wherein the molecular solutions are used to design synthetic 
pathways and predict product yields based on equilibrium constants calculated from the 
heats of formation. 

156. The system of claim 1 , wherein the new stable compositions of matter are predicted as 
well as the structures of combinatorial chemistry reactions. 

157. The system of claim 1 5 wherein pharmaceutical applications include the ability to 
graphically or computationally render the structures of drugs that permit the identification 
of the biologically active parts of the specie to be identified from the common spatial 
charge-density functions of a series of active species. 

1 58. The system of claim 1, wherein novel drugs are designed according to geometrical 
parameters and bonding interactions with the data of the structure of the active site of the 
drug. 

1 59. The system of claim 1, wherein to calculate conformations, folding, and physical 
properties, the exact solutions of the charge distributions in any given specie are used to 
calculate the fields, and from the fields, the interactions between groups of the same 
specie or between groups on different species are calculated wherein the interactions are 
distance and relative orientation dependent. 

160. The system of claim 159, the fields and interactions can be determined using a finite- 
element-analysis approach of Maxwell's equations. 

161. The system of claim 1 , wherein in the case where three sets of C — C -bond MOs 
form bonds between the two carbon atoms of each molecule to form a six-member ring 
such that the six resulting bonds comprise eighteen paired electrons, and each bond 
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comprises a linear combination of two MOs wherein each comprises two C2sp 3 HOs and 
75% of a i^-type ellipsoidal MO divided between the C2sp 3 HOs: 

3^2 C2$p 3 + 0.75 H 2 MOy~* 3 G = c)- ethylene - type - bond MO 

( 3e \ 

V-> 6\C=C/- bond MO of benzene 



(15.142). 



162. The system of claim 162, wherein the linear combination of each H 2 -type ellipsoidal 
MO with each Clsp 3 HO further comprises an excess 25% charge-density contribution 
per bond from each Clsp 2 HO to the C = C -bond MO to achieve an energy minimum. 

1 63 . The system of claim 1 63, wherein the dimensional parameters of each bond C = C- 
bond are determined using Eqs. (15.42) and (15.1-15.5) in a form that are the same 
equations as those used to determine the same parameters of the C = C -bond MO of 
ethylene (Eqs. (14.242-14.268)) while matching the boundary conditions of the structure 
of benzene. 

1 64. The system of claim 1 6 1 , wherein hybridization with 25% electron donation to each 
C = C-bond gives rise to the C hemen 2sp 3 HO-shell Coulombic energy 

E CoM (C be ^2s P >) given by Eq. (14.245). 

1 65. The system of claim 1 61 , wherein to meet the equipotential condition of the union of 

3e 

the six C2sp 3 HOs 5 c 2 and C 2 of Eq. (15.42) for the aromatic C= C-bond MO is given 
by Eq. (15.62) as the ratio of 15.95955 eV , the magnitude of E Cadomb (C betjzem ,2sp 3 ) (Eq. 

(14.245)), and 13.605804 eV , the magnitude of the Coulombic energy between the 
electron and proton of H (Eq. (1.243)): 

C 2 (benzeneC2sp 3 HOy c 2 (benzeneC2sp 3 s HO)^ ^f^^V " 0,85252 ( 15 * 143 ) 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 
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3e 

166. The system of claim 161, wherein the energies of each C=C bond of benzene are 
determined using the same equations as those of ethylene (Eqs. (14.251-14.253) and 
(14.319-14.333) with the parameters of benzene. 

3e 

1 67. The system of claim 161, wherein ethylene serves as a basis element for the C = C 

bonding of benzene wherein each of the six C=C bonds of benzene comprises 
(0.75)(4)= 3 electrons according to Eq. (15.142). 

1 68 . The system of claim 1 6 1 , wherein the total energy of the bonds of the eighteen 
electrons of the C = C bonds of benzene, E T (g 6 H 6 ,C=c), is given by (6)(0.75) times 

3e 

E T+osc (c = C) (Eq. (14.492)), the total energy of the C=C -bond MO of benzene 
including the Doppler term, minus eighteen times E(<J 9 2sp 3 ) (Eq. (14.146)), the initial 

3e 

energy of each C2sp 3 HO of each C that forms the C = C bonds of bond order two. 

3e 

169. The system of claim 161, wherein the total energy of the.six C=C bonds of benzene 
with three electron per aromatic bond given by Eq. (14.493) is 

E T (c 6 tf 6 ,C = c) = (6X0-75^ (c = dy (^)e(c^) 

= (6X0.75X"66.05796 eV)- 18(-14.63489 eV) (15.144) 
- -297.26081 eV - (-263.42798 eF) 
= -33.83284 eV 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to + 10%, if desired. 

1 70. The system of claim 161, wherein the results of benzene can be generalized to the 
class of aromatic and heterocyclic compounds. E hv of an aromatic bond is given by 
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E T (tf 2 ) (Eqs. (1 L212) and (14.486)), the maximum total energy of each H 2 -type MO 
such that 



E =n(E n + Er ih }=n, 

osc 1 \ D Kvib J 1 



-31.63536831 eV< 



l Mc 2 



(15.145), 



wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 

171 . The system of claim 170, wherein the factor of 0.75 corresponding to the three 
electrons per aromatic bond of bond order two given in the Benzene Molecule (C 6 H 6 ) 
section modifies Eqs. (15.52-15.56). 

172. The system of claim 171, wherein the multiplication of the total energy given by Eq. 
(15.55) by f x = 0.75 with the substitution of Eq. (15.145) gives the total energy of the 
aromatic bond: 



E T (Group)— f x 



Eibasis energies) + E T (atom - atom jnsp 3 .AO^ 



-31.63536831 eV 



2W 



Atvs R 3 



m 



m c 



- %KU It 

■ + n,E r ..+c, " 

1 Kvib 3 3 

r 



(15.146) 



wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 



1 73. The system of claim 161, wherein the total bond energy of the aromatic group 

E D (group) is the negative difference of the total energy of the group (Eq. (15.146)) and the 

total energy of the starting species given by the sum of c A E initial (c 4 AO I HO^ and 
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( AO i HO Y cArtu. *o / HO)) 



(15.147), wherein the calculated and measured values and constants recited in the 
equations herein can be adjusted, for example, up to + 10%, if desired. 

174. The system of claim 161, wherein benzene is considered as comprising chemical 
bonds between six CH radicals wherein each radical comprises a chemical bond between 
carbon and hydrogen atoms, energy components of V g9 V p9 T , V m , and E T are the same 

as those of the hydrogen carbide radical, except that E T (c = C,2sp 3 )= -1.13379 eV (Eq. 

(14.247)) is subtracted from E T (CH) of Eq. (13.495) to match the energy of each 

C - H -bond MO to the decrease in the energy of the corresponding C2sp 3 HO, wherein 
the calculated and measured values and constants recited in the equations herein can be 
adjusted, for example, up to ± 10%, if desired. 

175. The system of claim 174, wherein in the corresponding generalization of the aromatic 
CH group, the geometrical parameters are determined using Eq. (15.42) and Eqs. (15.1- 

15.5) with E r (atom - atom, msp\AO)= -1.13379 eV wherein the calculated and 

measured values and constants recited in the equations herein can be adjusted, for 
example, up to ± 10%, if desired. 

176. The system of claim 174, wherein the total energy of the benzene C-H -bond MO, 
E Tbe ^ {C~H) 9 given by Eq. (14.467) is the sum of 0.5E T (c = C,2sp 3 ), the energy 

change of each C2sp 3 shell per single bond due to the decrease in radius with the 



3e 



formation of the corresponding C=C -bond MO (Eq. (14.247)), and E (CH 



benzene ^ 




MO contribution given by Eq. (14.441). In the corresponding generalization of the 
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aromatic CH group, the energy parameters are determined using Eqs. (15.146-15.147) 
with f x ~l and E T (atom - atom 7 msp 3 .A O ^- l'-^379 e ^ ^ w i iere } n the calculated and 

measured values and constants recited in the equations herein can be adjusted, for 
example, up to + 10%, if desired. 

177. The system of claim 176, wherein the energy contribution to the single aromatic CH 

3s 

bond is one half that of the C= C double bond contribution, which matches the energies 

3e 

of the CH and C=C aromatic groups, conserves the electron number with the equivalent 
charge density as that of s = 1 in Eqs. (15.18-15.21), and further gives a minimum energy 
for the molecule. 

3e 

178. The system of claim 1 77, wherein the breakage of the aromatic C = C bonds to give 
CH groups creates unpaired electrons in these fragments that corresponds to c 3 = 1 in Eq. 
(15.56) with E mag given by Eq. (15.58). 

179. The system of claim 178, wherein each of the C - II bonds of benzene comprises two 
electrons according to Eq. (14.439). 

180. The system of claim 179, wherein from the energy of each C-H bond, 

-E D ( n CH) (Eq. (14.477)), the total energy of the twelve electrons of the six C - H 

benzene \ S 

bonds of benzene, E T (C 6 H 6 ,C - given by Eq. (14.494) is 

E T (C 6 H 6 ,C - H)= (^E D ^ ne (^CH^y 6(-3.90454 eV)= -23.42724 eV 

(15.148), wherein the calculated and measured values and constants recited in the 
equations herein can be adjusted, for example, up to ± 10%, if desired. 

181. The system of claim 1 74, wherein the total bond dissociation energy of benzene, 
E D (C 6 H 6 ), given by Eq. (14.495) is the negative sum of E T \C 6 H 69 C=Cj (Eq. 



(14.493)) and E T (f 6 H e ,C - h) (Eq. (14.494)): 
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(15.149) 



= 57.2601 eV 

wherein the calculated and measured values and constants recited in the equations herein can 
be adjusted, for example, up to ± 10%, if desired. 

182. The system of claim 181, wherein using the parameters given in Tables 15.214 and 
15.216 in the general equations (Eqs. (15.42), (15.1-15.5), and (15.146-15.147)) 
reproduces the results for benzene given in the Benzene Molecule (C 6 H 6 ) section as 
shown in Tables 15.214 and 15.216. 

183. A system of computing and rendering the nature of at least one specie selected from a 
group of diatomic molecules having at least one atom that is other than hydrogen, 
polyatomic molecules, molecular ions, polyatomic molecular ions, or molecular radicals, 
or any functional group therein, comprising physical, Maxwellian solutions of charge, 
mass, and current density functions of said specie, said system comprising: 

processing means for processing physical, Maxwellian equations representing charge, 
mass, and current density functions of said specie to produce at least one Maxwellian 
solution; and 

an output means for outputting the Maxwellian solution. 

1 84. The system of claim 1 83, further comprising a data transfer system for inputting 
numerical data into or out of a computational components and storage components of the 
main system. 

185. The system of claim 183, further comprising a spreadsheet containing solutions of the 
bond parameters with output in a standard spreadsheet format. 

186. The system of claim 185, further comprising a data-handling program to transfer data 
from the spreadsheets into the main program. 
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1 87. The system of claim 1 83, wherein output may be at least one of graphical, simulation, 
text, and numerical data. 

188. The system of claim 1 87, wherein the output may be the calculation of at least one of: 
(1) a bond distance between two of the atoms; (2) a bond angle between three of the 
atoms; (3) a bond energy between two of the atoms; (4) orbital intercept distances and 
angles; and (5) charge-density functions of atomic, hybridized, and molecular orbitals, 
wherein the bond distance, bond angle, and bond energy are calculated from physical 
solutions of the charge, mass, and current density functions of atoms and atomic ions, 
which solutions are derived from Maxwell's equations using a constraint that a bound 
electron(s) does not radiate under acceleration 

189. The system of claim 183, wherein the charge, current, energy, and geometrical 
parameters are output to be inputs to other programs that can be used in further 
applications. 

190. The system of claim 1 89, wherein the data of heats of formation can be input to 
another program to be used to predict stability (existence of compounds) equilibrium 
constants and to predict synthetic pathways. 

191 . The system of claim 190, wherein novel composition of matters can be discovered by 
calculating at least one of a bond distance between two of the atoms, a bond angle 
between three of the atoms, and a bond energy between two of the atoms, orbital intercept 
distances and angles, charge-density functions of atomic, hybridized, and molecular 
orbitals, the bond distance, bond angle, and bond energy being calculated from physical 
solutions of the charge, mass, and current density functions of atoms and atomic ions, 
which solutions are derived from Maxwell's equations using a constraint that a bound 
electron(s) does not radiate under acceleration. 

192. The system of claim 190, wherein the charge and current density functions can be 
used to predict the electric and magnetic fields of the species to determine other properties 
due to the interaction of the fields between species. 



122 



WO 2007/053486 



PCT/US2006/042140 



193. The system of claim 192 5 wherein finite-element analysis is used to predict or 
calculate the interaction and resulting properties, such as the freezing point, boiling point, 
density, viscosity, and refractive index. 

194. The system of claim 1 83, wherein the output data can be used to give thermodynamic, 
spectroscopic, and other properties, aid in drug design and other applications with or 
without direct visualization. 

195. The system of claim 1 83, wherein the data can be input into other programs of the 
system, which calculate thermodynamic and other properties, or performs a simulation, 
such as a chemical reaction or molecular dynamics. 

196. A method of using any of the systems of claims 1 - 1 95 . 

1 97. A use of any system of claims 1 - 1 95 . 

198. A composition of matter comprising a plurality of atoms, the improvement comprising 
a novel property or use discovered by calculation of at least one of 

a bond distance between two of the atoms, 
a bond angle between three of the atoms, and 
a bond energy between two of the atoms, 
orbital intercept distances and angles, 

charge-density functions of atomic, hybridized, and molecular orbitals, 
the bond distance, bond angle, and bond energy being calculated from physical solutions of 
the charge, mass, and current density functions of atoms and atomic ions, which solutions are 
derived from Maxwell's equations using a constraint that a bound electron(s) does not radiate 
under acceleration. 

199. A composition of matter according to claim 1 98, wherein the novel property is a new 
pharmaceutical use. 

200. A composition of matter according to claim 198 that is novel, wherein the novel 
property is stability at room temperature of a new arrangement of atoms or ions. 
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201 . A novel composition of matter discovered by calculating at least one of 
a bond distance between two of the atoms, 

a bond angle between three of the atoms, and 
a bond energy between two of the atoms, 
orbital intercept distances and angles, 

charge-density functions of atomic, hybridized, and molecular orbitals, 
the bond distance, bond angle, and bond energy being calculated from physical solutions of 
the charge, mass, and current density functions of atoms and atomic ions, which solutions are 
derived from Maxwell's equations using a constraint that a bound electron(s) does not radiate 
under acceleration. 

202. A system of determining at least one of 
a bond distance between two of the atoms, 

a bond angle between three of the atoms, and 
a bond energy between two of the atoms, 
orbital intercept distances and angles, 

charge-density functions of atomic, hybridized, and molecular orbitals, 
wherein the bond distance, bond angle, and bond energy being calculated from physical 
solutions of the charge, mass, and current density functions of atoms and atomic ions, which 
solutions are derived from Maxwell's equations using a constraint that the bound electron(s) 
does not radiate under acceleration. 

203. The system of claim 202, further comprising discovering a new composition of 
matter. 

204. A system of claim 1-195, wherein the calculated and measured values and constants 
recited therein can be adjusted up to ± 10% as desired. 



124 




1/3 



WO 2007/053486 



PCT/US2006/042140 



Fig. 2 




2/3 



WO 2007/053486 



PCT/US2006/042140 



Fig. 3 



9& Tetrachtorotnethane 
N'rtrogen Dioxide 
Hydrogen Nitride 
#£> Formic Acid 
m> Nitric Oxide Radical 

Triatomic Hydrogen MolectJar 
rg£3 Aldehydes 
3 23 Aromalics 
■SC3 Carboxylate_Esiers 
SiD Carboxylic_Aads 
Adds 

#e* Carbolic Add 
#a Adipic Arid 

SalleycllcAdd 
9& Lactic Add 
•e* Carbonic Arid 
ss-£3 Alcohols 
fid Amtna^Acids 
S-Ll Ketones 
sO Ethers 
•*! 23 Amines 
Si £3 Alkanes 
| ""-:'2J Functional Groups 
>j iJ Aromatlcs 
*0 Alcohols 
*CJ Ethers 
iil3 Nitno 
'*iZJ Amines 
*£2J Alkanes 
3d Amides 
0S> Alkyt Chloride 

Ketone 
m> Acid Haiide 
9g> Acetate Ester 
V> Acetic Anhydride 

i Formic Acetic Anhydride 



Wit 



3/3 



